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PREFACE 


The electric field which exists in the clouds is mainly responsible 
for enhanced ice crystal production within the supercooled fogs. Evans 
(1973) experimentally demonstrated that charged water molecules move 
to the crystal tips, thereby increasing the nucleation rates. 


Singh et al. (1986) and Sharma et al. (1992) discussed the effect of 
external electric field on the relaxation time in the nucleation process of 
water vapour condensation and ice glaciations. The effect of electric field 
on condensation and further polari 2 ability of water molecules have been 
discussed. (Singh, N., et. al. 2000, 2001). 


Effect of orographic features on atmospheric electrical parameters 
of different cities of India was discussed by Kumar et al. (1998). Shaw 
(1998) has shown the effect of electrical discharge in the nucleation of 
water droplets. 


Pinsky et al. (1999) showed that inertia of cloud droplets play a 
crucial role in the evolution of spectral during drop growth by diffusion. 
The change in electric field due to lightening in thunderstorm has been 
observed by Beasley and Eack (2000). 


OBJECTIVE: 


We have studied the effect of external electric field on the 


condensation of water and glaciations of ice in clouds. We have 
calculated the values of relaxation time by the help of Collins formula. 



ie: 


li 


0 ) 
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The equivalence between electric field with supersaturation ratio 
and Temperature has also been calculated for the present study in which 
spherical liquid drop model has been considered. 

According to this model, an embryo of water molecule is formed. 
The embryo acquires a critical size (called nucleus). At a certain 
supersaturation ratio and tempetature below this size, nucleus of water 
vapour molecules is formed. In presence of electric field, the Gibb’s free 
energy is modified due to the polarization of water vapour molecules and 
hence the nucleus of smaller size is formed earlier and hence the 
nucleation rate enhanced. 

THEORETICftL CONSIDERATION: 

In chapter-1, we have reviewed the nucleation phenomenon in 
presence and absence of electric field due to workers. The nucleation rate 
and the relaxation time in the electric field induced cases have not been 
compared for temperature and supersaturation with the electric field. The 
models for the nucleation and electrification of clouds have also been 
studied. 

Possible influence of cosmic rays on climate through thunderstorm 
clouds was studied by Lev. I. Dormana et al. (2005). Cosmic ray-induced 
Ionization in the atmosphere has been studied by LG. Usoskin (2004). 
Gondot et al. (1988) developed a new type of sensor to measure ionic 
electrical conductivity and, convection conductivity both around and 
inside convection clouds. Vonnegut (1988) proposed the mechanism to 
explain thunderstorm electrification differ greatly in how they might 
affect or be affected by the cloud. Electric field changes of lightning 
observed in thunderstorms had been studied by William H. Beasley 


(2000). During water vapour condensation the equivalence between 
supersaturation ratio and the external electric field has been established 
recently by Dwivedi, D. et al. (2008). 


In chapter-2, we have estimated the effect of external electric field 
on Homogeneous nucleation phenomenon of water vapour condensation 
and ice glaciations under the variation of density and surface free energy 
with temperature at different values of electric field. 

The frequency of growth of embryo beyond a critical size in the 
nucleation process corresponds to a free energy maximum, this energy is 

AGjy = (/Tj, fli)Tl^ + ^w/v^n (0 

Where is the chemical potential of vapour and /ij; chemical 
potential of liquid. (T^/v, the surface free energy of water vapour 
interface; and surface area of embryo; n^,, is the number of water 
molecules in the embryo. 


Corresponding to the free energy maximum for a critical size 
nucleus, we have, 

aAGw _ Q 

anw 


and thus the eq. (1) becomes 


• _ 20vv/j^ 4. n- n' 

T ^w/vH-w'hv 

3tIv\/ 


( 2 ) 


Where = to (3) 

m^, being the molecular mass and the density of water. 

{iii) 



In the absence- of an electric field Collins (1955) evaluated the 
relaxation time. 

To = (4) 

The number of water molecules in a critical nucleus is given by: 


n* = V 

\1.958xl0~sy 


Where is the radius of the critical nucleus with 


•w 


RT £tl Sp y| 


(5) 


(6) 


In this expression Sp.vi; is the supersaturation ratio of vapour over 
the plane water surface. 


In eq. (4) putting the values from eq. (5), eq. (6) and than 
evaluating the values of Tq, putting the values of m^= 3x10 gm.; 
=18; N = 6.023x10^^; k = 1.38xlO'’Srg K*’; R=8.317xl07erg K'* mole'*; 
n'^ = 2.m2xl0-^\ We have 


To = 3782.026 x 10'^ x 


<Tvv/v 




(7) 


Variation of density and surface free energies of water and ice with 
temperature has taken into account (Edward et al. 2003; Gupta et al. 
2009). 

In presence of an electric field we have 

2/3 


r; = [3p„(9 « XEym„y'Hn/2p^] 


(S) 



In chapter-3, we have established equivalence between 
supersaturation ratio and the external electric field during water vapour 
condensation. Nucleus of a particular radius is obtained fr om eg. (6), 
whereas in presence of electric field the radius is given by eg. ( 8^ If the 
size of the nucleus is obtained in two cases, than R.H.S. of tiiese two 

t . — — — 

equations may be equated to get. 


2A^wO'w/v 

Py^RT -STI S-iy yy 


3py,TE 

2pyy \myyj J 


( 9 ) 


Which reduces to 


4 

X 



9p„T 

\py^CC2,J 

\RT in Sv.w^ 


From eq. (10) we may write 



( 10 ) 


(II) 


Where E = Eg^, is equivalent electric field corresponding to 
supersaturation ratio 5^ vv- 


From eq. (1 1) we can write 


£■2/3 33 X X f ) (12) 

\9pvrJ \pwaV \RT in ^ ^ 


Which reduces to 


X X (13) 

from eq. (11) and eq. (13) we have calculated the values of 
equivalent electric field and equivalent supersaturation ratios. 


(V) 



From the equations we came to know that the supersaturation ratio 
varies exponentially with equivalent electric field. 

Also the nucleation rate J and Gibb’s free energy for the formation 
of water and ice phase have been calculated. 

In chapter-4, we have established an equivalence between 
temperature and the electric field during water vapour condensation and 

ice glaciations. ' 

< 

In presence of an external electric field, the rate of increase of 
radius with respect to time is given by Singh, N. (1986). 

In absence of electric field the radius of the critical nucleus is given 

by 

But in presence of external electric field, the same critical radius is 
given by 



3PyTjV 

2p^ 



(15) 


In electric field induced nucleation, for a given size of the critical 
nucleus, we have described a relation between temperature and electric 
field, using eq. (14) and eq. (15) as 

(16) 

or £'73/2=Kfor a given supersaturation ratio and given relaxation 

time. 


r 2My^cr-^fv 

3/2 

^2pu/^ 

1/2 1 

LPvv/jRT ■^72 


' TjvT3/2 


In chapter-5, we have described the effect of ions on the electric 
field nucleation phenomenon as well as the study of thunderstorms has 
been dealt with in atmospheric clouds. We have also discussed the 
electrical conductivity of air in thunderstorms, variation of diffusion 
coefficient of ions, attachment coefficient and studied the use of 


(Vi) 



atmospheric electrical conductivity as an air, pollution indicator through 
correlation analysis based on two year of daily averaged observations of 
the atmospheric electrical conductivity. 

The chief ionizing agents in the lower atmosphere are the 
penetrating radiation from the radioactive substances in the earth’s crust, 
radiation from the radioactive substances in the atmosphere and the 
cosmic rays. 

In the Retezat Mountains concentrations of O3, NO2 and SO2 in 
summer season 2000-2002 were calculated and were low and below 
toxicity levels for forest trees. Ambient O3 levels found in the Retezat did 
not affect crown condition of Norway spruce or European beech (Akula 
Venkatram, 2006). 

Condensation of water vapour of ice germ formation in 
supercooled liquid, when ionizing radiations are passed (Varshney, 1969) 
is a milestone in the understanding of ions as centers of condensation. 
This basis has been confirmed theoretically (Vashneya, 1971)., The 
pioneer worker (Castleman, 1979; Mohanen, 1970, 1974; Chan, 1980a, 
1980b) has placed the belief on a firm footing. 

IMPORTANCE OF PROPOSED STUDY: 

The work proposed for investigation is very important as it may 
lead to air explanation of atmospheric problems related with nucleation. 
Thunderclouds have an electric field developed inside which causes 
nucleation. The several charging machines have been proposed. The 
present research work proposed would set up the effective role of external 
electric field in development of clouds and also set up a tool for artificial 
rain making. 


(vii) 
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CHAPTER- 1 

INTRODUCTION 


In atmosphere there are a number of gases and water vapour are 
present in different proportion at different temperature and pressure. The 
amount of water vapour in the air determines how fast each molecule will 
return back to the surface. The amount of water vapour in an atmosphere 
is constrained by the restrictions of partial pressures and temperature. The 
water vapour exist in all three forms i.e. Gas (vapour), Liquid (water) and 
Solid (ice). The phase change of water vapour takes place the central role 
in cloud microphysics in the following manners- 

Condensation 

1 . Water vapour ^ . — ^ Liquid water. 

Evapouration 

Condensation Crystalization 

2. Water vapour > Liquid water > Solid water (ice). 

Glaciation 

3. Water vapour > Solid water (ice) 

Approx. 99.99% of water vapour is contained in the troposphere. 

1.1 CHEMICAL POTENTIAL: 

To every phase certain energy in the form of latent heat is 
associated with it. This energy is called the chemical potential energy. 
When the phase change takes place there is a certain amount of energy 
given out or absorbed. According to “Willard Gibbs” an American 
physicist chemical potential '\l' is defined as follows - “If to any 
homogeneous mass in a state of hydrostatic stress we suppose an 
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infinitesimal quantity of any substance to be added, the mass remaining 
homogeneous and its entropy and volume remaining unchanged, the 
increase of the energy of the mass divided by the quantity of the 
substance added is the potential for that substance in the mass 
considered.” 

Chemical potential is also referred to as partial molar Gibbs 
energy. 

1.1 (i) Chemical Potential of water vapour in Humid Air: Let us 

derive chemical potential of water vapour and through the equilibrium the 
chemical potential of water in aqueous solutions. We have equation 
(Maxwell relations) 


^ ^^"k/T,p,nj9iK 


( 1 . 1 ) 


From eq. (1 . 1) we have for an ideal gas K in mixture of ideal gases 



So that, upon integration 

l^K = hK,o + RT^npk (1.3) 

Where the integration constant depends only on the 

temperature. For such a mixture the potential pressure is = X^p, so 
that also 


M-k = Rk,o + RT^np + RT^nXk 


(1.4) 
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Therefore if we assume pure water vapour at pressure e is an ideal 
gas, its chemical potential is 

Hv,o(e,T) = <o(T) + RT^ne (1.5) 

Where gy.oC^) is the chemical potential at a standard state of unit 
pressure. Similarly the chemical potential py.CP'T) of water vapour in 
humid air at total pressure p is 

gv(p,T,Xv) = <o(T) + RT^np + RT^nX^ (1.6) 

From eq. (1.5) with e = p, we see that this last result may also be 
expressed as 


^v(p,T,Xv) = gv,o(P»T) + RTfnXv (1.7) 

Which shows that Py ^ Pv,o since Xy < i 


In contrast to pure gases whose chemical potentials vary 
logarithmically with pressure, the chemical potential of a pure liquid is 
proportional to pressure, to an excellent approximation. This is obvious 
from. 


■ ^T,nj=iK 


Vk 


( 1 . 8 ) 


On realizing that liquids are nearly incompressible thus for water 
we have 



Constant 


(1.9) 


(thesis of devesh dwivedi) 


From which the chemical potential is found to be 

Uw.oCP-T) « Pw,o(0»T) + Vw,oP (1-10) 

As we have seen, if a liquid and gas are in equilibrium the chemical 
potential of a given component will be the same in both phases, 
consequently, from eq. (1.3) the chemical potential of component K in a 
liquid solution, which is in equilibrium with its vapour at partial pressure 
Pk, is 

Pw = Pk,o+ RT-^npk (1.11) 

In addition, experiments show that for so-called ‘ideal’ solutions, 
for which there are no interactions between the solvent and solute 
molecules, the equilibrium vapour pressure of any component is 
proportional to its mole fraction in the solution. (This is known as 
Raoult’s law). 

Assuming Raoult’s law then we have Pk = Xk,^Pk,o, where Xk,^ is 
the mole fraction of component K in the solution, and Pk,o, is the partial 
pressure of component K in equilibrium with the pure liquid phase of K 
at the same temperature. Then as a fraction of Xk,^, the chemical potential 
becomes 


Pk,^ = Pw + RT^nXk,, (1.12) 

Where Pk,^ is a function of both temperature and total pressure, but 
is independent of the composition of the solution. 
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1.1 (ii) Chemical Potential of water in aqueous solutions: In clouds the 
liquid phase is rarely present in the form of pure water, but rather is 
generally a dilute aqueous salt solution. 

Therefore, + RT-^nXk,^ is especially relevant to us, 

and we may use it to write the chemical potential for water in an ideal 
aqueous salt solution in the following form: 

liw(p.T,X^) = p+(p,T) + RT^nX^ (1.13) 

Where X^^ = = 1 — Xg is the mole fraction of water, n^ 

(Hw + ^^s) 

and ng being respectively the number of moles of water and salt in the 
solution. By analogy one would expect that the chemical potential of the 
salt component could be expressed in the same way, viz; 

Hs(P.T,Xw) = nJ-Cp.T) + RT^nX (1.14) 

In passing, we may note that for 

Xw = l,Pw(PtT) = Pw,o(P«T) = Pw(PtT) 

The chemical potential of pure water at (p,T), there is no analogous 
simple physical interpretation for the quantity Pg (p,T). 

Experiments shows that most dilute solutions of non-electrolytes 
are in conformity with eq. (1.13) & eq. (1.14). In general, however, real 
aqueous solution depart from such ideal behaviour. It is customary to 
account for non-ideal solutions through the replacement of the mole 
fraction X by the activity a = fx, where f is called the rational activity 
coefficient. Thus for real aqueous salt solutions we write 

PwCPr'f'^w) ~ Pw,o (P' "^3 "b RT■^na■^y; 3-^^ (1.15) 


Ps(p-T,as) = Ps (PrT) + RTfnag; ag = fgXg 


( 1 . 16 ) 
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The importance of the activity to us is that it provides a direct 
measure of the equilibrium water vapour pressure over a real salt 
solution, or in other words, the generalization of Raoult’s law to real 
solution we now turn to a demonstration of this property. 

1 .2 CHANGE m CHEMICAL POTENTIAL: 

Change in chemical potential of a system is given by 

dH=-(^v-M (1-17) 

Where |iv is the chemical potential of the vapour phase and the 
chemical potential of liquid phase. This change in chemical potential 
appears as the energy of the formation of the liquid phase. 

1 .3 GIBBS PHASE RULE: 

This rule enables us to determine the variance of a system i.e. the 
number of intensive variables, which may be jfreely specified without 
cause the system to depart from equilibrium. In equilibrium the system is 
characterized by a common T&p, and by a number of mole fractions in 
the various phases. 

Let us denote the mole fraction of the k* component in the j* phase by 

= n®/2kn^’ 

Thus, for every phase we have the simple mass conservation 
constraint that = 1 for a total of (j) constraints. In addition we 


have the condition 


THESIS OF DBVESH DWTVEDI 


Uk = = •”= 1,2, 3, 4,- , c. On the 

chemical potentials, which constitute another (cj)-!) constraints for every 
k, for a total of c( cj) -1) constraints. Therefore the number of 
thermodynamic degrees of freedom, is 

w = 2 + (|)C - (|) - c((|) - 1) =2 + c — (|) (1-18) 

Which is the Gibb’s phase rule for the bulk phases. As a simple 
application of above equation consider a homogeneous fluid, we have- 

4> = 1, c = 1 and so w = 2. 

For a mixture of two gases, c=2, 4» = 1 and hence w=3. For water 
in equilibrium with a vapour and ice, c=l, cj) = 3 and w=0, equilibrium is 
possible only for a single choice of temperature and pressure, which 
defines the triple point temp. T|r of this system. If this system is now 
exposed to the atmosphere, c=2 (water substance and air), 4> = 3 and so 
w=l. 

1 .4 STRUCTURE OF AN ISOLATED WAT^ MOLECULE: 

Heat capacity of water vapour at constant volume near room 
temperature yield a value of approximately 3k per molecule, where k is 
the Boltzmann Constant. We must interpret the heat capacity 
measurements in terms of a contribution of (V 2 ) k from each of the three 
translational degrees of freedom. This interpretation implies that the 
water molecule cannot have its three atoms arranged in a linear fashion, 
which is the same conclusion by investigating electrical properties of the 
water molecule. In water each hydrogen nucleus is bound to the central 
oxygen atom by a pair of electrons that are shared between them, these 
shared electron pair are called covalent chemical bond. Since such 

CD 


measurements reveal a large electric dipole moment of v= 1.83x10'^* 
e.s.u. cm. (Hobb’s, P.V., 1974; Eisenberg et al., 1969). 

On the basis of such measurements, it has been concluded that the 
three atoms are situated at the vertices of a triangle, the geometry of 

which is given in fig. (1.1). Recent experiments show that the equilibrium 

0 

0-H bond length is 0.95718 A and that the equilibrium H-O-H bond 
angle is 104.523° (Fletcher, N.H., 1970a; Hobb’s, 1974). 

In the structure of water molecule in its ground state, an oxygen 
atom has two electrons in the spherical Is orbital, where they are bound 
tightly to the atomic nucleus, and two electrons, less tightly bound in the 
spherical 2s orbital (fig. 1 .2). In addition, two electrons can be considered 
to occupy the 2Px orbital, one electron in 2Py orbital and one electron in 
2Pz orbital. Since the 2Py and 2Pz orbitals may contain two electrons 
each, these are incomplete. Therefore the electrons in these orbitals are 
free to couple with the Is orbital of two hydrogen atoms, allowing them 
to form two 0-H bonds. If these describe them to form two 0-H bond of 
a water molecule one would expect water to have a bond angle of 90°. 
However experimentally, one find that the bond angle is some 15° larger. 
Since oxygen is more electronegative than hydrogen, oxygen exerts a 
greater force on the shared electron pair. Electrostatic repulsion between 
the two hydrogen atoms must, consequently, lead to an increase of the 
bond angle. 

However, Heath and Linnet (1948) showed that this repulsion is 
insufficient to account for the experimentally, found bond angle. They 
suggested that the ‘hybridization’ of the 2s orbital of the oxygen atom; 
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fig. 1.2 : Hybridization of orbitals of oxygen atom. 
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with its 2py and 2pz orbitals, resulting in the formation of four Sp^-hybrid 
atomic orbitals. Two of these overlap with the hydrogen orbitals, while 
the two remaining orbitals form two lobes on the side of the Oxygen atom 
away from the Hydrogen atoms fig. (1.2). These lobes called lone-pair 
hybrids, are symmetrically located above and below the molecular plane 
and form roughly tetrahedral angels with the bond hybrids. (Exact 
109.467°). It is this tetrahedral character of the water molecule which 
gives rise to a tetrahedral coordination of water molecules in water and 
ice. 


Duncan and Pople (1953)& Bader and Jones (1963) have carried 
out quantum mechanical (molecular orbital theory) calculations of the 
electron density distribution shown in fig. (1.2), and show that there are 
four locations in the water molecules with high electron density; close to 
the Oxygen atom, close to each Hydrogen atom, and at the location of the 
lone pair orbitals, which appear as an electron density bulge behind the 
oxygen atom. 

The charged distribution around a water molecule may also be 
approximated by ‘electrostatic’ or point charge models. In these models 
charges are assigned whose sign, magnitude, and location are such that 
the molecules as a whole is electrically neutral, and the electric dipole 
moment is equal to that experimentally measured. Such models have been 
worked out (Rowlinson, 1951; Bjerrum, 1951; Pople, 1951; Cambell, 
1952; Cohen, 1962). Although such models are convenient in same cases, 
they generally do not correctly predict the higher electric moments. The 
current views of water structure is influenced greatly by molecular 
modeling simulations beginning in the 1980’s, water is more like a ‘gel’ 
consisting of a single, huge hydrogen-bonded cluster. On a 10' "-10' sec 
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time scale, rotations and other thermal motions cause individual hydrogen 
bonds to break and re-form in new configurations, inducing ever- 
changing local discontinuities whose extent and influence depends on the 
temperature and pressure. 

1.5 STRUCTURE OF Water VAPOUR: 

In contrast to ideal gas behaviour experiments show that water 
molecules in water vapour tend to interact & form clusters (Kell, 1972a). 
Though only in small concentration dimmers as well as higher order 
polymers are considered to be present in water vapour. Recent 
experiments involving molecular beam techniques (Lin, 1973; Searcey, 
1974) suggest that in highly supersaturated water vapour, clusters of up to 
180 water molecules may be present. Clusters of 21 water molecules 
seemed to exhibit particularly large stability. It is interesting to note that 
21 water molecules can be arranged in the form of a pentagonal 
dodecahedron with a molecules at each comer and a single molecule in 
the centre of the ‘cage’. 

At present however, no conclusive evidence of the actual geometric 
arrangement, if any, of water molecules in such clusters in vapour is 
available. Studies on the possible and more likely cluster types have been 
reviewed (Rao, 1972; Kell, 1972a).Recent theoretical studies of the 
formation of water clusters have been carried out (Kistenmacher et al., 
1974a, 1974b; Abraham, 1974). Kistenmacher et al. found two possible 
stable configurations for the dimmers, a cyclic form and an open form, 
which was more stable. For the large clusters the authors suggested not a 
single stmcture, but a statistical distribution of different configuration. 
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Between a pair of water molecules the potential energy of 
interaction U, has the general character of being strongly repulsive at very 
close separations and weakly attractive at longer range, so a simple 
expression for it is due to Stockmayer (1941). 


H^f C , Cor^^ 

r3 1-6 ‘ y24 


(1.19) 


Whei;e r is the separation of the molecules; p, is the dipole moment 
of an isolated water molecules; o is the collision diameter (the molecular 
separation at which U==0 if p=0), c is an adjustable constant, and f is 
known function of the mutual orientation of the two molecules. 


The first term on the right of eq. (1.19), is just the dipole-dipole 
contribution to the interaction energy, and may be attractive or repulsive. 
The second term represents contribution from- 

1. The interaction energy between a permanent dipole of one 
molecule and the dipole it induces in the other (dipole polarization 
or induction interaction). 

2. The net energy arising from momentary, fluctuating dipoles 
interacting with the corresponding induced dipoles (polarization-2 
or dispersion interaction). Even though the time average of these 
dipole fluctuations may be zero, the energy contribution is 
proportional to their mean square, which is finite and positive. 

Both (1) & (2) are usually referred to as Vander Wall’s interaction, 
which by its nature can be seen to bring about an attractive force between 
the molecules. The third term in eq. (1.19), represents the short range 
repulsive forces, which loosely may be ascribed to the overlap of 
electronic orbitals, which are incompatible according to the Pauli 
exclusion principle. 
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There is little doubt that the Stockmayer potential or similar ones 
such as (Rowlinson, 1949, 1951), potential portray with fair accuracy the 
interaction between pairs of water molecules at large separations in dilute 
water vapour. This is evidenced by the fact that values for the second 
virial coefficient computed via. eq. (1.19), can be made to fit 
experimental values on the other hand, the same potential functions yields 
values for the third virial coefficient of water vapour, which disagree 
substantially with experiment. Partly, this is due to the approximate 
nature of eq. (1.19), and partly because three-body interactions should be 
included also, since other molecules in the system can significantly 
modify the interaction of a given pair. In particular the Stockmayer 
potential is insufficiently ‘directional’ in character to account for the 
geometry of cluster formation in water vapour. 

1.6 STOUCTURE OF ICE: 

At temperature between about -100°C & 0°C and at atmospheric 
pressure, water substance crystallizes from its gaseous or its liquid state 
to form a six fold symmetric or hexagonal solid called ice-Ih. At different 
pressure and temperatures ice assumes other crystalline modifications, 
which are discussed, for example in Fletcher (1970a) and Hobbs (1974), 
we shall concern ourselves here only ice-Ih, hence fourth referred to 
simply as ‘Ice’. 

X-ray diffraction studies demonstrate that in ice each oxygen atom 
is surrounded by four nearest neighbour oxygen atoms at a distance of 
about 2.76 x 10'® cm. These four atoms form an almost regular 
tetrahedron. These tetrahedrons are joined together to form a hexagonal 
lattice fig. (1.3).These hexagonal space group is denoted by or P 63 
/mmc and is characterized by 1 six-fold axis of rotation. Perpendicular 
to 1 mirror plane, (3+3) mirror planes, and a centre of symmetry. 





fig. 1.3 : Position of Oxygen atoms in ice-Ih: (a) view along c-axis, (b) view 
perpendicular to c-axis, (c) Hexagonal Ice-Structure. 




Near O'^C, any given atom in ice also has 12 second nearest 

0 

neighbours at a distance of about 4.52 A,i third nearest neighbour at 4.59 

0 0 

A, 6 fourth nearest neighbour at 5. 26 A, 3 fifth nearest neighbour at 5.31 

0 0 

A, 6 sixth nearest neighbour at 6.36 A^ 6 seventh nearest neighbour at 

0 0 

6.46A, 9 eighth nearest neighbours at 6.69A, 2 ninth nearest neighbours 

0 0 

at 7.36 Aand 18 tenth nearest neighbours at 7.81 A . 

Each water molecule is hydrogen bonded to its four nearest 
neighbours. The hydrogen atoms in a water molecule are about 104 
degrees apart, which is just about the right angle to build tetrahedral 
structures. The positively charged portion of one water molecule is 
attracted to the negatively, charged portion of a neighbour. This 
attraction, termed hydrogen bonding, gives water many of its imique 
properties. This is brought about through the formation of two hydrogen 
(O-H-0) bonds by each water molecule, each bond being directed 
towards a lone electron pair of a neighbouring water molecule. This 
manner of bonding leads to an open lattice structure, as illustrated in fig. 
(1.3). Perpendicular to the c-axis, the ice lattice consists of open-puckered 
hexagonal rings with oxygen atoms. Comparisons show that the 
arrangement of oxygen atoms in ice is isomorphous with the wurtzite 
structure of ZnS and the tridymite structure of SiOa- 

Each unit cell of ice, a four sided prism set on a rhombic base, 
contains four water molecules and is characterized by the lattice constants 
ao and Cq fig. (1.4). X-Ray data for ao and Cq (Blackman,1957; 
Lonsdale,1958; Laplaca,1960; Brill, 1967 ; Kumai,1968) are summarized 
in fig. ( 1.5), as a function of temperature. 




fig. 1.4; A typical disordered arrangement of protons in the ice-Ih structure; 

oxygens (1) and (2) contribute 12/12 each oxygen (3), (4), (5), 
(6) contribute 1/12 each, oxygens (7) to (12) contribute 2/12 
each, and oxygens (13) and (14) contribute 4/2 each, for a total of 
48/12 =4 oxygens. 






These measurements show that ao and Co decrease with decreasing 
temperature- the rate of decrease smaller, the lower the temperature-such 
that (co/ao) = 1.629 for all temperatures. Using the values for ao and Co 
given in fig.(1.5), the volume of a unit cell of ice Vuj. = 2(ao-73/4) cq, 
varies from 1.305 x lO'^^cm^ (0"C) to 1.281 x 10'^^ cm^ (-180°C). Thus 
the number of water molecules cm" varies from 3.06x 10 (0 C) to 

3.12x10^^ (-180®C). According to Eisenberg and Kauzmann (1969), three 
points of view may be taken to define the hydrogen bond energy Eh in 
ice. First, one may assume that Eh is given by lattice energy El of one 
mole of ice (the difference in energy between one mole of isolated water 
molecules and one mole of ice, both at 0®K and with motionless atoms) 
divided by the number of hydrogen bonds in a mole. 

Since both hydrogen atoms of a water molecule participate in one 
H-bond (excluding the molecules in the surface of ice), one may estimate 
Eh as Eh = El/ 2 = 7.04 K cal mol’*, based on the value for El found by 
Whalley et al. (1973). Second one may define Eh, possible more 
appropriately, in terms of sublimation enthalpy ( AH )s by writing 
Eh=(AH)s/2 Since (AH)s is temperature dependent, so is Eh, and we find 
Eh (0°K) = 5.66 K cal mole* and Eh (0“C) = 6.15 K cal mole*. In both of 
these definitions we ascribe the entire intermolecular energy in ice to 
hydrogen bonding we therefore include in Eh the effects of dispersion and 
short range repulsive forces which are present not only in the ice but also 
in crystals of non hydrogen bonded substance. The third definition of Eh 
is based on premise that the contribution to (AH)s from the hydrogen 
bonds is distinct from that of other forces. One may therefore set 
Eh = [(AH)s — Eother]/2. where Eother represents the intermolecular 
energy associated with the other forces. This definition suffers from the 
fact that Eother is not an observed quantity and cannot presently be 
accurately calculated. 




An accurate theoretical calculation of a single hydrogen bond in ice 
should include at least the effects of nearest neighbours. Generally the 
approach taken has been to assume the total hydrogen bond energy is 
given by the sum of the four component energies (dipole-dipole, dipole- 
polarization, polarization-polarization and short range inert-actions) and 
to evaluate each of these by approximate methods for two neighbouring 
water molecules at the relative positions found in ice. For this purpose 
various models for the charge distribution in a water molecule have been 
assumed. 

The results of the most pertinent calculations on the subject have 
been summarized (Hobb’s, 1974; Eisenberg, 1969). The values computed 
for Eh range from 4 to 8 k cal mole'*. In all of these computations, the 
entire intermolecular energy has been attributed to hydrogen bonding. 

The positions of hydrogen atoms in ice and subject to the Bermal- 
Fowler (BF) rules these require that- 

1 . Each water molecule is oriented such that its two hydrogen atoms 
are directed approximately towards two of four oxygen atoms 
which surround it tetrahedrally, 

2. There is only one hydrogen atom on each 0-0 linkage and, 

3. Each oxygen atoms has two nearest neighbouring hydrogen atoms. 
Such that the water molecule as structural unit is preserved. 

An ice structure, which obeys BF rules, is termed ideal. Natural 
ice, however, does not behave ideally. The following major aromatic 
defects are found in natural ice (stacking faults), chemical defects; 
molecular vacancies (Schottky defects), interstitial molecules (Frenkel 
defects), ionized states and orientational defects (Bjerrum defects). 
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1 .7 STRUCTURE OF IHTATER AND AQUEOUS SOLUTIONS: 

In ascribing ‘structure’ to a fluid such as water the time periods ty 
and Xd just mentioned becomes especially relevant. In fact, we must 
consider three different time scales times t« Xy times intermediate to 
TyandTo, att » Tq. Assuming we were equipped with a camera which 
had shutter speeds less than Xy, we could obtain a relatively sharp picture 
of the actual position of a water molecules at any given instant. 

This would reveal the instantaneous water structure called the I- 
structure. If the shutter speed were between Xy and xd, each molecule 
would complete many oscillations while the shutter was open and the 
resulting somewhat blurred picture would provide information on the 
vibrationally averaged position of the water molecules in water, i.e. V- 
structure (Eisenberg and Kauzmann) of water, if the shutter sp eed were 
larger than Xd, the diffusionally averaged arrangement of the water 
molecules of D-structure could be found. Experimental studies which 
employ infrared or Raman spectroscopy, or neutron scattering techniques 
lead to information on the V-structure, while X-ray studies determine the 
D-Structure. 

Spectroscopic studies show that the frequency of oscillation for 
water molecules is slightly smaller in water than in ice, the period of 
vibration being Xy « 10“^^ second. Studies on self diffusion, viscosity, 
dielectric relaxation, and nuclear magnetic resonance relaxation show that 
a water molecule in water has a characteristic displacement period near 
0°C of Xd » sec. 

Information on the state of hydrogen bonds in water can also be 
obtained from infrared and Raman spectra (Walrafen, 1966, 1967, 1968a, 
1968b, 1972).Such spectra confirm that water molecules exists as entities 
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in water. They also give evidence that some 0-H groups in H 2 O are 
hydrogen bonded and, therefore, point towards a free, lone electron pair 
of a nighbouring H 2 O molecule, while other 0-H groups are non- 
directionally bonded to the surrounding water molecules and hence are 
disoriented with respect to neighbouring lone electron pairs (Kell, 1972a; 
Eisenberg, 1974).The latter are referred to as non-hydrogen bonded or 
‘broken’ 0-H groups. Estimates of the percentage of broken H-bonds in 
water as function of temperature are summarized in fig. (1.6). 

These experimental findings are supported by recent studies, which 
attempt to simulate the molecular structure of water by purely theoretical 
methods (Rahman, 1971; Stillinger, 1972; Popkie, 1973; Kistenmacher, 
1974a, 1974b). Applied conventional molecular dynamics to a system of 
216 water molecules, which interacted via a potential function developed 
(Ben Naim and Stillinger, 1972; Ben Naim, 1972).This pair potential 
function is considerably more complicated than the Stockmayer potential 
eq. (1.19). It is based on Bjerrum’s four point charge model for a water 
molecule and incorporates the linear bonding tendency between 
neighbours in a tetrahedral pattern such as that found through out the ice 
lattice. It has been argued (Stillinger, 1970, 1972) that this potential 
function also incorporates the principal features of non-additively; i.e. it 
takes into account the many body aspect of the problem. 

Recently, another pair potential function was developed (Clementi, 
1973; Fromm, 1975; Popkie, 1973; Kistenmacher, 1974a, 1974b; 
Abraham, 1974).This function is based on an analytically fitted Hartree- 
Fock potential (Kern, 1972), the point charge model and correlation 
energy corrections due to induced dipole interaction and short-range 
effects. It was^used in conjuction with the Monte- Carlo simulation. 
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fig. 1.6 : Variation with temperature of number of hydrogen bonds broken in 
water. 
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method (Barker, 1969), to study a group of 125 water molecules. In 
addition, the molecular dynamic study (Stillinger, 1972), predicted that 
the hydrogen bond rupture mechanism in water is characterized by an 
excitation energy of about 2.5 K cal mof'. This is in good agreement with 
the Raman and infrared spectra experiments (Walrafen, 1966, 1967, 
1968a, 1968b, 1972), who found that the energy necessary to rupture 
hydrogen bonds in water is about 2.55 Kcal Mole"'. Infrared spectra taken 
(Worley, 1966), imply the comparable value of 2.4 Kcal mor\ Finally 
(Bucaro, 1971), inferred 2.5± 0.1 Kcal mof’ from depolarized light 
scattering measurements. 

Considerable uncertainties still exist as to how X-ray, infrared, and 
Raman Studies should be interpreted in terms of the arrangement of the 
water molecules in water. We shall now describe briefly just the main 
features of some of the more prominent models put forward for the water 
structure. 

In the ‘quasi-crystalline model’ the water structure is assumed to 
resemble one of several possible forms; a broken-down ice. Ih structure, a 
quartz structure, a structure of octahedrally arranged, molecules (Van 
Eck, 1958), or a structural mixture of molecules arranged in a tridymite 
structure dispersed in a denser ice-III structure (John, 1966).In the 
‘interstitial model’ the water is visualized as consisting of a highly 
hydrogen bonded structure inside of which non bonded or partially 
bonded molecules occupy interstitial structure position (Forslind, 1952; 
Danford, 1962; Krestov, 1964).In the ‘Clathrate model’ water is assumed 
to have a structure similar to the clathrate structure of gas hydrates except 
that instead of a gas molecule a water molecule is held inside each cavity 
of a cage like, hydrogen bonded frame work of Pentagonal dodecahedron 
cages (Pauling, 1959, 1960; Frank, 1961). 
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The ‘flickering cluster model’ makes use of the partially covalent 
character of the hydrogen bond and assumes that H-bond formation in 
water is a cooperative phenomenon, in that the formation of a hydrogen 
bond between two water molecules reinforces that tetrahedral 
hybridization in oxygen atoms. 

The ‘mixture model’ picture water as a mixture of 0-, 1-, 2-,, 3- dnd 
4- bonded water molecules engaging in the formation of various sized 
clusters (Frank, 1957; Walrafen, 1972).Finally, the ‘bent bond model’ 
assumes that few, if any, bonds between water molecules are broken upon 
melting of ice, but instead become bent to various degrees (Pople, 1951). 

Experiments show that the structure of water is altered when water- 
soluble salts, in part dissociated into ions, are dissolved in water. The 
aqueous solution resulting from dissolving a salt in water would be an 
ideal solution if the dissolved salt molecules or ions in no way affected 
the water molecules (Haggis, 1952). 

In any real aqueous solution this is not the case. For example, some 
of the salt molecules or ions do not fit into the water ‘structure’ and 
therefore distort it, causing a size effect. This effect is called hydration. 
Since large ions have weaker local electric field than small ions, the 
hydration effect is greater for small ions. In addition, hydration is more 
pronounced for positive ions than for negative ones, since a positive ion 
tends to interact with both ions electron pairs, which blocks the formation 
of two H-0 group of a water molecule, which blocks the formation of 
only hydrogen bond. 
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1 .6 Water Vapour distribution in india: 

Sarkar et al. (1982) brought about an atlas on water vapour 
distribution over India as a part of their programme under the project 
“Anomalous, Microwave Propagation and Prediction of Radar 
Performance over India.” The list of sixteen regular radiosonde stations 
providing data is given in table (1.1). 

Water vapour density was calculated from randiosonde data of four 
year (1968-71) using fixed pressure level data from lOOOmb to 50mb 
level. Radiosonde flights were taken twice daily, one at 0000 GMT and 
the other at 1200 GMT corresponding to early morning and late evening 
in local time. Accuracy of relative humidity temperature and pressure 
obtained from radiosonde measurements are 5%, 0.25% and 2mb, 
respectively. Over the northern plains, the water vapour density is 
generally high during monsoon and post monsoon periods while the value 
are of moderate order during winter and summer months. 

The water vapour density values are obtained from two 
observations per day; hence these values do not strictly indicate the 
percentage of time in a year when a particular water vapour density can 
be expected. It seems unlikely that the extremes of diurnal density 
variation could always occur at OOOOGMT and 1200 GMT, and 
therefore cumulative distribution statistics have a tendency to be the case 
if the observations were made frequently, our data are especially useful 

V,..,.,..., 

for station to station in two regions, such as estimating the water vapour 
values for one region fig. (1 .7) to fig. ( 1.14). 
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Table-1.1 

List of Radiosonde Stations 


Minico 


Visakhapatnam 

3. Cochin 

4. Calcutta 

5. Bombay 

6. Madras 

7. Gauhati 

8. Ahmedabad 

9. Trivandrum 

10. Port Blair 


1 1 . Lucknow 


14. Nagpur 


15. Bangalore 


16. Srinagar 


Elevation (m 


2 




1666 





fig. 1.7 ; Water Vapour density profile over Northern Plains in Winter 
(OOOOGMT). 
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fig. 1.8 : Water Vapour density profile over Northern Plains in Pre monsoon 
(OOOOGMT). 
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fig. 1.9 : Water Vapour density profile over Northern Plains in Monsoon 
(OOOOGMT) 




altitude (Km) 
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fig. 1.10 : Water Vapoiir density profile over Northern Plains in Post 
monsoon (OOOOGMT). 
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fig. 1.11 ; Water Vapour density profile over Northern Plains in Winter 
(1200 GMT). 
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1.12: Water Vapour density profile over Northern Plains in Pre monsoon 
(1200GMT). 
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fig. 1.13 : Water Vapour density profile over Northern Plains in monsoon 
(1200GMT). 
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fig. 1.14 : Water Vapour density profile over Northern Plains in Post 
monsoon (1200GMT). 




1 .9 SUBFACE PBOPERTIES OF WATEB SUBSTANCE: 


Surface Tension:- Phases in contact are separated by a thin 
transitional region, generally only a few molecule thick; consequently, a 
useful abstracting regard. Such an interface as geometrical surface. To 
include surface effect, in complete analogy to the contribution-pdV of 
pressure volume work the internal energy contribution is given by ads, 

f 

where S denotes the area of the surface of separation and a the surface 
tension. The quantity a is an intensive thermodynamic variable and has 
dimension of energy per unit area or force per unit length. 

1.9 (i) Effect of Temperature On Surface Tension: The surface 
Tension of water decreases with increasing temperature, which is 
excepted on considering the effects of thermal agitation. This behaviour 
has been investigated experimentally by Dorsch and Hacker (1951) their 
results are shown in fig. (1.15). For the interval o to 40®C the data can be 
reproduced to within ± 0.02 dyne cm"’ by the following expression; 
o’w/a = 76.10 - 0.155 T (1.20) 

Where a^fa is in dyne cm"’, i.e. Erg cm"^, and T is in O^C. From 
fig (1.5), it appears that eq. (1.18), also represents a reasonable 
extrapolation for the temperature interval 0 to - 40’’C. 

1.9 (ii) Radius Dependence of Surface Tension: Surface Tension 
arises ^fom attractive forces between molecules near the surface, hence 
we may expect that only and attraction of average geometrical 
configuration of these molecules on a size scale comparable to the 
effective range of the attractive forces would significantly affect the 
Surface tension. Thus we expect a dependence of a^f.^ on size only for 

QD 



(thesis of devesh dwivedi^ 


extremely small drops, consistency of nearly a few tens or hundreds of 
water molecules. The surface tension of water decreases with decreasing 
radius of curvature of water surface. Tolman (1949a, 1949b),estimated 
that cTiv/v for a drop, which consist of 13 water molecules (equivalent to a 
drop radius of 4.6 x 10’^ cm, based on Vw = 30xl0'^‘*cm^ and pw = 1.0 gm 
cm'^) is 40% smaller than that for a plane water surface. Benson and 
Shuttle worth (1951) computed this surface tehsion of a small drop water 
molecules by counting the number of bonds, which had to be broken in 
order to cut off the group of molecules from the bulk water structure. In 
order to estimate the interaction energy between water molecules in water 
it was assumed that only the first and second nearest neighbour has to be 
considered. In this manner they predict that the surface tension of a drop 
of 13 water molecules is only 15% smaller than that for a plane water 
surface. 

Following Defay et al. (1966) a simple, approximately quasi 
thermodynamic derivation of radius dependence of Uw/i; is given by fig. 


(1.15) 




(<^w/v)oo 

i+{3iC7(Pw-pD] 


( 1 . 21 ) 


(-) 

Where (cTw/v) is surface tension of a plane water surface, 

f-) 

absorption and the term /(Pw ~ Pv) is independent of radius ‘a’ 
table (1.2) lists result of cr^^/v (3^) computed from eq.(1.21).Using 
» 0.87 X 10“^ mol cm'^ and approximatio « p^ = lg.cm~^. 
It is seen that the radius dependence becomes important for a « 
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fig 1.15 : Variation of Surface tension of water (against air) with temperature 
according to Gittens & Dorsch & Hacker (1951). 




Variation Of Surface Tension Of Water For Drops Of Various Radii 


a(cm) 

^Wjv! (^W/p) QQ 

[(<^w/ v)^ ®’iv/i;] / (s^wfv) ^ 

’ (%) 

00 

1.0 

0 

10-^ 

0.9997 

0 


0.9969 

0.3 

10'" 

0.9697 

3.0 

10'^ (140 molecules) 

0.7622 

23.8 

4.6x 1 0’**(1 Smolecules) 

0.5959 

40.4 

- 
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1.9 (iii) Surface Tension at Vapour Interface: The difference 
between the pressures insider water drop of radius ‘a’ and the pressure 
Csat-w of vapour with which it is in equilibrium is given by 

Pw-esat.w=^~ (1.22) 

Where now we have introduced the subscript w/v for the surface 
tension to emphasize that is the water vapour interface which is involved. 
Given that ^76 dyne cm'^ at O^^c, we see that the pressure 
difference is about 1.5 atm. For a = 1/rm , smaller drops have 
correspondingly larger internal pressure. 

For practical purposes one may replace by Uw/a > surface 
tension for a water-humid air interface. Experimentally have shown that 
increases by less than 0.05% if air at 1 atm is replaced by pure water 
vapour at saturation pressure (at the same temperature). 
l^(iv) Surface Tension at Ice-Vapour Interface and Surface 
Energy of Ice: Surface energy of ice is the energy required to form a 
unit area of low surface. For ideal crystalline ice energy may be identified 
with one half of the energy per unit area, which is needed to split an 
infinite crystal parallel to a particular crystallographic alone and separate 

Wr 

the two part by one infinite distance. It is natural to take — , as the 
surface tension or interfacial energy between the particular ice 
crystal face and water vapour or air, assuming the presence of such gases 
does not affect the surface energy. Disregarding the forces of interaction 
due to third or higher order nearest neighbour, energy per molecule 
required for cleavage of an ice crystal can be expressed as 
£■^ = 1 / 1 + 61/2 


(1.23) 


Where Ui & U 2 are the average interaction potentials between the 
molecules in the first and second interaction zones. The interaction 
potentials Ui & U 2 for intermolecular spacing’s of 2.76x10'* cm. and 
4.51x10'* cm, respectively, were computed (Reuck, 1957) on the basis 
(Rowlinson, 1951) offeree constants, which take account of the multiple 
electrostatic interaction forces, induction forces and repulsion forces. As 
a result, Reuck found NAEn,= 6.08 Kcal mole'* or = 4.22x1 0'^* erg per 
molecules, for an ice cystal in vacuum at 0*^K. 

The estimated values of Wc for different shapes of ice crystals are 

= 238 erg cm'^ for the basal faces of ice and erg cm'^ for 

the prism faces; the corresponding values for the surface energies are 

ri 5 ') 9 2 

= 1 19 erg cm'" and =126 erg cm* . Similar estimates were made 

by Me. Donald (1953b). 

Me. Donald (1953b) has pointed out that the surface energies thus 
computed pertain to “fresh” surface since molecules in a freshly cleaved 
surface will not remain in their original position but will relax into new 
equilibrium position, surface energy of an “aged” surface is somewhat 
less than that of fresh surface. Significance of such relaxation can be 
appreciated if we calculated for water at 0°C in the same way as was 
just done for ice and compare results with the known values of surface 
Tension. 

1.9 (v) Ice water Interface (State between ice and water): 

Consider a system of consisting of an ice crystal surrounding by super 
cooled water. For sufficient small particles of ice such a system can be in 
stable equilibrium so that it is again possible to speak of surface tension 



as interface energy of boundary separating the phases unfortunately it is 
very difficult to determine by experimental technique, and the 
results of numerous attempts show considerable spread, fig. (1.16) 
(Wood, 1970; Dufour, 1963; Jacobi, 1955a; Ketcham, 1969; Kubelka, 
1944; Turnbull, 1950, 1965). 

A reasonable relation for the variation of with the temperature 

f 

as given by 

28.5 + 0.25 T (1.24) 

Where a is in erg cm'^, and T is in ®C. Direct photographic 
evidence of the presence of steps at a growing ice-water interface has 
been provided (Ketcham, 1968;Hobb’s, 1969). 

Some of the observed steps had spiral forms appearing at a number 
concentration of about 10^ cm‘^, with step heights between 0.1 and 4/im, 
and spacing between the steps of 5 to 20jum. No explanation for this 
observation is currently available except to say that, as in the case of step 
formation at the ice-air interface, some sort of bunching mechanism may 
be operating. 

I.IONUCLEATON: 

The formation of a new phase at the expense of a metastable 
original phase (“Mother phase”) does not begin in a continuous manner, 
but rather takes place spontaneously as result of temperature and density 
fluctuations (heterophase fluctuation) in the original phase, provided that 
a critical supersaturation is exceeded. This spontaneous process is called 



0 ~ ^/w Cerg cm 



fig. 1.16 : variation with temperature of the surface tension of ice against 
supercooled water (1) Wood and Walton (1970) for spherical ice 
crystal, (2) Wood and Walton (1970) for Wulff crystal, (3) 
Dufour and Defay (1963) based on data of jacobi (1955 a), (4) 
Ketchman and Hobbs (1969), (5) Pruppacher (Unpub.) based on 
data of Kubelka and Prokscha (1944), (6) Pruppacher (Unpub.) 
based on (Turnbull, 1950, 1965). 
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nucleation. Thus the formation of critical nuclei of water or ice is called 
the nucleation phenomenon. Before the critical size the embryos are 
formed of different sizes. Beyond the critical size the nuclei are changed 
into drops by further condensation of water vapour. But below the critical 
size the embryos may be vapourized or dissociated. 

The nucleation phenomenon has been divided jn tvvo categories- 

1 . Homogeneous Nucleation 

2. Heterogeneous Nucleation. 

^ Homogeneous Nucleation of Water: Homogeneous 

Nucleation occurs when supersaturation ratio is several hundred 
percentages and it may be assumed that foreign substance are absent but 
in actual atmosphere only 0.01 and to 0.10 supersaturation (i.e. 1% to 
10%) is present in the clouds, hence homogeneous nucleation is not 
possible generally. Nucleation generally occurs with much more 
difficulty in the interior of a uniform substance, by a process called 
homogeneous nucleation. The liquids which are cooled above the 
homogeneous nucleation temperature are said to be supercooled. 
Quantitative comparison of experiment and theory for the growth of 
homogeneously nucleated water droplets had been explained by 
Alexander Fladerer et al. (2003). 

1.10 (ii) Heterogeneous Nucleation of Water: In the case of 
heterogeneous nucleation some energy is released by the partial 
destruction of the previous interface. For example, if a carbon dioxide 


C53 


bubble forms between water and the inside surface of a bottle, the energy 
inherent in the water bottle interface is released wherever a layer of gas 
intervenes and this energy goes towards the formation of bubble water 
and bubble-bottle interface (Nickolay, 2005). Heterogeneous nucleation 
is generally responsible for cloud formation. Such nucleation is assisted if 
suitable solid surface or some foreign materials like dust particles, cloud 
seeding particles and ions are present. The heterogeneous nucleation also 
takes place when external electric field is applied to the cloud. 

Thus heterogeneous nucleation may be of the following types, 

a. Ion induced nucleation 

b. The electric field induced nucleation 

c. Nucleation on certain solid surfaces of different shapes and 
materials. 

Daniel Beysensa, b, (2006) explained that Dew is the condensation 
of water vapour into liquid droplets on a substrate. It is characterized by 
an initial heterogeneous nucleation on a substrate and a further growth of 
droplets. 

1.11 CLOUD CONDENSATION NUCLEI (CCN): The drop 
formation in the atmosphere occurs via heterogeneous nucleation 
involving aerosol particles (A.P.). A.P. which are capable of initiating 
drop formation at the observed low supersaturation are called cloud 
condensation nuclei (CCN). 


Cloud condensation nuclei or CCNs (also known as cloud seeds) 
are small particles (typically 0.0002mm or 1/100* the size of a cloud 
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droplet. Water requires a non-gaseous surface to make the transition from 
a vapour to a liquid. In the atmosphere, this surface presents itself as tiny 
solid or liquid particles called CCNs. When no CCNs are present water 
vapour can be super cooled below O^C (32®F). The number of CCNs in 
the air can be measured and ranges between around 100 to 1000 per cubic 
cm. The total mass of CCNs injected into the atmosphere has been 
estimated at 2x1 0*^ kg over a year’s time. 

Aerosols play an important role in a variety of processes in the 
marine boundary layer. They transfer water vapour, heat and matter 
through the air sea interface, interact with the fields of temp and humidity 
by evaporation and condensation, and may act as condensation nuclei in 
the formation of clouds and fog (Van Eijk, 1996). 

J.P.D. Abbat and Broekhuizen et al. (2005), the ability of mixed 
ammonium sulfate/organic acid particle to act as ‘cloud condensation’ 
nuclei (CCN) has been studied in the laboratory using a continuous flow, 
thermal gradient diffusion chamber operated at supersaturation at 
between 0.3% and 0.6%. The organic acid studied were malonic acid, 
azelaic-hexanoic acid, cis pionic acid, oleic acid and stearic, and the 
particles were largely prepared by ‘Condensation’ of the organic vapour 
onto a dry ammonium sulfate core. For malonic acid and hexonoic acid, 
the mixed particles activated as predicted by a simple Kother theory 
model where both species are assumed to be fully soluble and the droplet 
has the surface tension of water. Three low solubility species, cis-pinonic 
acid, azelaic acid, oleic acid are well modeled where the acid was 
assumed to be either partially or fully insoluble. We observed no CCN 
behaviour that could be clearly attributed to lowering of the surface 
tension of the growing droplet by the presence of the organic constituents, 
some of which are highly surface active. 
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1.12 ICE FORMING NUCLEI (IN): 


The cloud glaciations generally begins at temperature which too 
warm for homogeneous freezing of water. For example at one occasion 
(Mossop, 1968), observed ice crystals in a long lived cumulative cloud 
whose top was probably never colder than -4®c and which was not seeded 
with ice particles from clouds at higher altitude. Such behaviour indicates 
that some fraction of the local Aerosol particles (AP) also, can serve as 
ice forming Nuclei (IN). 

1.13 NUMBER CONCENTRATION OF CCN: 

Twomey and Wojciechowski(1969), comprehensively studied 
CCN concentration over various part of the world. The results of these 
studies are summarized in fig. (1.17 a,b,c). The result shows no 
systematic latitudinal variation in concentration. The observation also, 
confirm that continental air masses are generally richer in CCN than are 
maritime air masses with in a particular air mass at flight level the 
variation of median CCN concentration was surprisingly small. At 
supersaturation between 0.1 and 10% the median concentration of CCN 
was found to range from a few ten to a few hundred cm in air over 
oceans and from a few hundred to a few thousand cm‘^ in air over the 
continents fig. (1.17). It is indicated that the concentration of CCN (i.e. 
Nccn) increases with increasing supersaturation Sy.w This behaviour can 
be expressed in the form 

^CCN — C^v.w 

Where k and C are approximately constant. 

r«3 


(1.25) 
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fig. 1.17 : Median world-wide concentration of CCN as a function of 
supersaturation required for activtion; (a) in air over oceans, (b) 
in air over continents, (c) all observation (Twomey, 1969). 
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At a given location the CCN concentration is magnitude, 
depending upon approximately of sources and an metrological factors 
such as wind direction, air mass type, precipitation and decreasing or 
increasing cloudiness (Twomey, 1959a; Rodke, 1969). 

Fig. (1.18) illustrates a typical time variation of CCN concentration 
(Twomey et al. 1970, 1971) showed that at given location, repeated 
patterns can be detected in the diurnal variation of the CCN 
concentration. This is a noon maximum and late evening maximum were 
observed during one year observation period at Robertson (N.S.W. 
Australia). Generally CCN concentration in maritime and modified 
maritime air masses which have been overland less than two days rarely 
exceed lOOcm'"’ are found in air which have been overland for several 
days. 

Variation of CCN at different values of supersaturation are 
discussed by Jiusto (1968) shown in fig. (1.19). Fig. (1.20) describes the 
comparison between total concentration of aerosol particles and 
concentration of cloud condensation nuclei activated at 0.7% 
supersaturation of various locations (Hoppel, 1973). In table (1.3), the 
values of Ailken particles per cm^ are given which were evaluated by 
Allee et al. (1970) and Twomey et al. (1964) at different places. 



■ 

This dip corresponds to the last storm cell recorded on the 
barograph 


Wind NW lOmph. light snow and fog-clearing 


Cu Breaking west 

4 . 

Wind NW5-1 Omph. light snow and fog. 

5 . 

Wind SW 5mph. Snow and fog. 

6 . 

Wind SW 2-5 mph light snow and fog. 

7 . 

Wind WS W 5- 1 0 mph Breaks in Cu. 


Wind WSW 5-1 Omph Fog 

9 . 

Wind WNW 5-10 mph Clearing to wast 

10 . 

Wind W 2-5 mph light fog 


fig. 1.18 : Variation with time of the CCN Concentration activated at 1% 
supersaturation during March 29, 1968, in air at observation 
station (2025m) in Olympic Mts., Washington State (Rodke, 
1969). 
















fig. 1.19 ; Variation of the CCN concentration required for activation at 
various locations, as a function of supersation (Jiusto, 1968). 
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Table-1.3 


Type of nuclei 

Location 

Number of Ailken 
Paticles 
(cm-’) 

Number ofCCN 
(cm’^) 

Washington D.C. (Alice, 

1970) 

7800 

2000 

6800 

2000 

5700 

5000 

50000 

7000 

Yellow Stone National 
Park (Wyoming) 

1000 

15 

Long Island (N.Y.) 

(Twomey, 1964) 

5700 

30 

6500 

150 

18000 

no 

51000 

220 




1 .14 A LABOEATORV STUDV OF ICE NUCLEATION DUE TO 


ELECTRICAL DISCHARGE: 

Experiments were performed (G.Mandal, 2002) in a walk- in cold 
room where the lowest attainable temperature is -30°C to investigate the 
effect of electrical discharge and high electric fields on ice nucleation. 
Experiments were carried' out at cloud temperatures varying from -9 to 
-12*’C. It was observed that ice nucleation does occur inside the cloud if 
there is an electrical discharge inside the cloud. However, not all 
electrical discharges resulted in ice nucleation. In no case was ice 
nucleation observed in the presence of high electric fields with no 
discharge. No ice nucleation was noted in the cloud when the Liquid 
water content (LWC) was greater than 2.1 gm'^ 

1.15 EXPERIMENTAL VRRMCATION OF HEIEROGENEOUS 

Water drop nucleation: 

Once water soluble and mixed A.P., grow with increasing relative 
humidity, it is appropriate to test the theory of nucleation of such particles 
in two ways: 

1. By experimentally determining the equilibrium growth of these 
particles and, 

2. By determining the number of particles of given size and 
composition which become activated at a given supersaturation. 
The agreement with theory is satisfactory as regards the shape of 
GCN Spectral actual concentration of activated drops. 
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Attempt to verify the prediction of water drop nucleation on water 
insoluble partially wetteable substrates have been made (Isaka, 1972; 
Mahata, 1975). 

These experiments were carried out in diffusion chamber where 
nucleation was forced to take place on plane substrate for which angle of 
contact of water had been measured separately. Considerable 
disagreement, while, the' critical supersaturation necessary for onset of 
drop nucleation are sufficiently lower for all wetting angles than those 
predicted by theory. 

The results of experiments carried out (Mahata, 1975) are 
reproduced in fig. (1.21). It is seem that the measured critical super 
saturations as a function of contact angle agree fairly well with those of 
Koutsky et al. (1965), up to a supersaturation of 10% corresponding to a 
contact angle of about 25°C. 

An overview of the two FEBUKO aerosol cloud interaction field 
experiments in the Thuringer Wald (Germany in October 2001 and 2002 
and the corresponding modeling project MODMEP is given by Hermann 
et al. (2005). Experimentally, a variety of measurement methods were 
deployed to probe the gas phase, particles and cloud droplets at three sites 
upwind, downwind and within an arographic cloud with special emphasis 
on the budgets and inter conversions of organic gas and particle phase 
within 30 cloud events three events (E I, E II, E III) are selected for 
detailed analysis. At various occasions an impact of the cloud process on 
particle chemical composition such as on the organic compounds content, 
sulphate and nitrate and also on particle size distributions and particle 
mass is observed. Moreover, direct phase transfer of polar organic 
compound from the gas phase is found to be very important for the 
understanding of cloud water composition. 



CRITICAL SUPER5ATURATI0N {%) 



Theory 

a]:j = 10*cm'^ sec’ 

3 . 2 '.} = Icm'^sec’ 

b:j = Icm'^sec'’ 

(With size correction to surface tension) 

fig. 1.21 ; Critical supersaturation for onset of water nucleation on a planjie, 
waterinsoluble, partially wettable substrate, as a function of 
contact angle of water on substrate. Comparison of experiment 
with theory (Mahata, 1975; Koutsky, 1965). 
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1.16 £L£CTE1CAL INflOMOeENETIES OTHER 'raAN IONS: 

The surface of the nucleation substrate consist of sharply defined 
boundaries between surface regions of different field sign, or of locations 
where the electric field vector in the substrate is oriented parallel to the 
surfaces. Growth of water cluster around electric inhomogeneties other 
than ions attract water molecules to the substrates and the dipole of a 
water molecule. Growth of water cluster at such sites can be aided if the 
difflisivity of water molecules of surface region with either an inward 
directed or an outward directed electric field is high as compared to the 
diffusivity over in a area where the electric field vector in the substrate is 
oriented parallel to the substrate, since than both the positive and negative 
ions of water molecules are partially tied down. 

Thomas^aluated both individual and combined effects of the 
horizontal and vertical variability of cumulus clouds on solar radioactive 
transfer are investigated using a two-dimensional (x-and z- directions) 
cloud radar data set. This high resolution dataset of typical fair weather 
marine cumulus is derived from ground based 94 GHz cloud radar 
observations. The domain averaged (along x-direction) radiative 
properties are computed by a Monte Carlo method. It is shown that (i) 
different cloud-scale resolutions can be used for accurate calculation of 
the mean absorption, upward and downward fluxes, (ii) the resolution 
effects can depend strongly on the solar zenith angle, and (iii) a few cloud 
statistics can be successfully applied for calculating the averaged 
radiative properties. 
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1.16 (i) Ice Forming Nuclei: Ice fonning nuclei exhibit three basic 
modes of action; 

1 . In the first mode, water is absorbed directly from vapour phase on 
to the surface of IN where, at sufficiently low temperature, 
observed vapour is transformed into ice. 

2. In the second mode, the IN initiates the ice phase from inside in 
supercooled water drop. 

3. In third mode in the action the IN initiates the ice phase at the 
moment of contact with the supercooled droplets. 

The first mode action is called deposition mode and AP which exhibit 
this behavior are called “deposition IN”; the second is called freezing mode and 
the corresponding ‘ AP’ are “freezing nuclei”, the third is called contact mode, 
participated in by contact angle. 

1.16 (ii) IN Concentration: The IN concentration measurements are 
reproduced as a function of temperature and location in fig. (1.22, 1.23, 
1.24), illustrates a day to day variation in IN concentration observed with 
one counting technique one particular site. Even at remote locations such 
as Antarctic, pronounced variation in the daily IN concentration do occur. 
Similar observation were made (Kikuchi, 1971), which were measured at 
an Antarctic station. IN concentration maxima as high as 30 at 20°G. 
Such high concentration counts were found to last up to two days fig. 
(1.21), also shows that the IN counts increase nearly exponentially with 
decreasing temperatures. 

Fletcher (1970a) gave a convenient statement of this behaviour as, 
Nif) - Aexp {fi.KT) (1.26) 

Where /? = 0.6 (°C)'', A = 10"^ where is the number 
concentration of IN active at a temperature warmer than T and where 


AT = To-T. 
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fig. 1.22 : Variation of the mean or median number conceutration of IN with 
temperature and geographic location. (1) Bracknell (England), 51* N, 0*W, 
(2) Clermont-Ferrand (France), 46* N, 3“ E, (3) Corvallis (Oregon, U.S.), 44“ 
N, 123“W, (4) Tokyo (Japan), 36*N, 140* E, (5) Tucson (Arizona, U.S.), 32* N, 
111"W, (Q Jerusalem (Israel), 32* E, (7) Palmbeach (Florida, U.S.), 2rN, 80* 
W, (8) Hawaii ^.S.), 20* N, 158” W, (9) Swakopmund (S- Africa), 34* S, 14* 
E, (10) Sidney' (Australia), 34“ S, 151* E, (11) Tasmania (Australia), 43“ S, 
147" E, (12) Antarctica, 78* S, 166* E. The dashed line represent Nm = 10^ 
exp (0.6 AT). 
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fig. 1.23 : Variation of the IN concentration with (a) relative humidity and 
with (b) supersaturation over ice. (Hulfman, 1973 a). 




Gagin (1972) and Huffman (1973a, 1973b) found that, at any given 
temperature the number concentration of IN increases with increasing 
relative humidity fig. (1.23a) and correlates logarithmically with the 
Supersaturation over ice, independently of temperature fig.(1.23b) 
according to the relation. 

Nin = CSfi (1.27) 

‘ Where K and C are constant. 

Huffman (1973b) observed that the vertical IN concentration 
profiles often exhibit a pronounced layer structure even in the lowest few 
Kilometers above ground. For example a strong concentration maximum 
was observed during a few days of summer time sampling at about 500 
meter above ground over St. Louis and over N.E. Colarado, no 
temperature inversion was present at the time over N.E. Colarado and the 
inversion over St. Louis was considerably above the IN concentration 
maximum, thus the concentration maxima appear not to be co-related 
with the temperature inversion. Bigg et al. (1970, 1963) studied range of 
median number concentration of IN as function of temperature at various 
geographic locations and the variation of the daily number concentration 
of ice forming nuclei at -20*^ C measured during December, 1961 to Jan., 
1962at Antarctica station fig. (1 .24) and fig (1 .25). 

1.17 ENHANCEMENT IIATIO OR ENHANCEMEST FACTOR: 

The enhancement ratio is defined as a ratio of ice crystal 
concentration determined at the cloud top temperature. This ratio can be 
as large as 10'^ to 10^ at temperatures between -5*^0 to -15®C. Also, 
tends to decrease with decreasing cloud top temperature and reaches unity 
between -25*^0 and -30°C. Fig. (1.26) (Hobb’s, P.V., 1974) represents 
some observations of over the cascade meters during winter of 1971 to 
1973, which were summarized by the best fit relation. 



number of ice-forming nuclei (m 
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Fig. 1.24 : Range of median number concentration of IN as function of 
Temperature for various geographic location; 44 stations. The 
dashed line represent Njn =10'^ exp (0.6AT) (Bigg, 1970). 
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fig. 1 .25 : Variation of the daily number concentration of ice forming nuclei at 
-20'^C measured during December, 1961, to January, 1962, at 
Antarctica station. (Bigg, 1963). 


maximum ice enhancement ratio, r 




CLOUD TOP TEMPERATURE i ^ C ) 


fig. 1 .26 ; Maximum ice enhancement ratio (Rm) in clouds over the Cascade 
Mts. (Washington State) during winter, 1971-72. (If value is 
below the detection threshold, it is show on the threshold line 
with an arrow attached.) Asaltostratus, S-stratus, Sc- 
stratocumulus, Cu-cumulus, C (S)-Cumulus with stratified tops. 
(Hobb’s, P. V., 1974). 




Log Rm = 5.02 + 0.204 Tc (1 .28) 

Where Tc is called cloud top temperature in ®C. 

1.17 (i) Ion Induced Nucleation: The observed anomaly of maximum 
enhancement ratio (~10'^ — 10^) can be explained when we consider 
the effects of ions on ice formation. The ratio can be written as. 



Where N;on is the concentration of ions (~10® to near 

the ground surface; IN, the concentration of ice forming nuclei at cloud 
top temperatures (~10^m“^); oc, the efficiency of ions acting as IN. As a 
rough estimate the value of oc can be taken as 0.1. For the above 
parameters the values of R[„ comes out to be (~10^). This value is in 
excellent agreement with the experimental observations. 

The ions provide a central force field, which is a simple form of 
attractive potential resulting in the formation of heterogeneous embryos. 
A basic understanding of ion induced nucleation, therefore, provides an 
important basis for further development in the theory of heteromolecular 
nucleation which is important in atmospheric process (Allen, 1969; 
Doyle, 1961). 

Condensation of water vapour on ions was also mentioned (Mason, 
1971). The possibility of ice germ formation in supercooled liquid, when 
ionizing radiations are passed through (Varshneya, 1969) is a milestone 
in the understanding of ions as centre of condensation. The pioneer works 
(Varshneya, 1971; Castelman, 1978; Mohiien, 1970, 1974; Chan, 1980a, 



1980b) has placed the belief on a firm footing. Ion induced nucleation 
process via a sequence of clustering reactions. These may be considered 
to commence with the third order reaction. 

A^ + B+M<r^A^-B + M 

M, being a third body. Here A- represents the ion (of either sign) 
about which molecules B is clustering. 

In the classical model of drop formation the Gibb’s function has 
been taken into consideration. The energy of formation of a nucleus is 
given by 

= -4nAG^ r^/3 + Ana^ (1-30) 

Where rw is the radius of water nucleus (assumed spherical), 
the surface free energy per unit area of water nucleus in water vapour; 
AGp, the volume free energy of condensate per unit volume per mole. In 
presence of ions the water vapour molecules form clusters with ion at the 
centre. In this ease the Gibb’s free energy is modified as eq. (1 .31). 

AGw = -(AnAG^Tw/S) + Anaw/v + (2'e)^/5rv(, (1.31). 

The last term in eq. (1.31) is due to the electrostatic energy of the 
cluster in the presence of an ion of charge (Ze), Z being the charge 
number, e the electronic number. 

Yuhang Wang and Russel (2005) positive matrix factorization 
(PMF) was used to infer the sources of PM 2.5 observed at four sites in 
Georgia and Alabama. 



1.17 (ii) Electric Field Induced Nucleation; Thunder clouds are known 
to have strong electric field. During a lightening discharge the maximum 
electric field (~10 e. s. u.) is produced near the channel. The electric field 
so generated affects the rated condensation of water vapour and ice 
deposition. It has been reported (Jalaluddin, 1962; Parmar, 1973, 

1 975)that the electric field induces nucleation in a liquid. 

The supercooled drops of a few millimeter diameter resting on a 
chromium surface froze immediately when electric fields of 20 to 60 
Kv/cm were applied. It has been observed that the supercooled water 
drops on a plastic surface froze when expressed to sparks from a 
Teslacoil. Sparks from the coil had no effect on the drop when they were 
not on a surface but suspended in the air (Evans, 1973; Singh, N., 1986). 

Now we can also write the phase change of water Vapour in 
following manners 

Condensation 

1. Water vapour ^ Liquid water (condensation) 

Evaporation 

Condensation Crystalization 

2 . Water vapour > Liquid water > Solid water (ice). 

Glaciation 

3. Water vapour > solid water or (ice) 

Nucleation is the key element in the understanding of many natural 
and industrial processes as well. It plays a prominent role in a wide 
variety of engineering process. 



1 .1 6 MODELS OF NUCLEATION: 


1.18 (i) Classical Liquid drop Model: Within meta-steady state bulk 
phase of water vapour of small molecular clusters of liquid water, which 
results from the ions agglomeration of water molecules; these are 
generally referred to as embryos if the vapour pressure is below the 
critical value required for nucleation. Such embryos have small binding 
charges and usually disrupted by thermal agitation. However, at the 
critical vapour pressure some embryos will reach critical size, at which 
point will be unstable equilibrium with the mother phase. A germ will 
proceed to grow spontaneously and thereby produce a microscopic 
phase, if a result of fluctuation in the mother phase; its size increase by 
even an infinitesimal amount. 

The workers description completely assumed that 

1. Pre-nucleation embryos may be regarded water spheres 
characterized by the usual macroscopic densities and surface 
tension. 

2. They are distributed according to the Boltzmann law. 

3. The i-mers are in thermal equilibrium, and the probability that they 
have the certain energy (A0)i , is just the probability for their 
existence, if we interpret (A0)i , as the energy of formation of i- 
mer. 

In this model embryos are considered as spherical of certain radius. 
The phase charge occurs at constant temperature. The total volume of the 
system considered (the mother phase + the condensed phase) remains 

on 


constant. 


Volmer (1939) and co-workers evaluated the free energy of 
formation of water clusters by ascribing macroscopic thermodynamics 
properties to them. Thus the free energy of droplets is described as the 
sum of the surface term. An attractive feature of this approach is that it 
permits ready visualization of the origin of the free energy barrier to 
nucleation in terms of the maximum in the above sum as a function of 
size. However, a very unattractive aspect is that the calculated size of 
critical nucleus i.e. The cluster size at a top of the free energy barrier is 
only about 100 molecules which leads to doubt the applicability of 
macroscopic concept in the present examples following standard methods 
(Frenkel, 1946). This spherical drop models leads to directly to a rather 
simple expression for the rate of homogeneous nucleation of droplets 
from vapour. 

The remarkable arrangement (Hirth, 1963; Volmer, 1939) of the 
macroscopic theory with observations of the critical supersaturation for 
appreciable nucleation rate of various liquids in cloud chambers stood for 
many years as the basis of our knowledge of nucleation. For example, 
referring to the data of Wilson (1899) and Powell (1928) who reported 
“fog limit” for homogeneous nucleation of water droplets at a 
Supersaturation ratio is about 5.0 at 273® K, agreement with microscopic 
theory is excellent (Feder, 1966). However, it has been pointed out 
(Feder, 1966; Dunning, 1969; Lothe, 1969; Frenkel, 1946; Hirth, 1963) 
that the external partition fractions for free translation and rotation had 
been neglected in the macroscopic theory. According to Lothe and Pound 
(1969, 1962, 1966, 1968) the Helmholtz free energy of formation should 
be corrected to include the following contributions; 

AF' = AFi + AFiT + AFiR - AFi_rep 
" '( " 69 


(1.32) 
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Where Fjx and FjR results from the translational and rotational 
degrees of freedom of the embryo. The term called replacement 

term by Lothe and Pound (1962) has been regarded as resulting from the 
deactivation of the vibrational degree of freedom, which the cluster 
would have possesses as bulk water (Lothe and Pound, 1962). Later, 
Lothe and Pound (1966) approximated AFi_fep as due to translation and 
vibrations about the centre of the embryo water mass considered again as 
a spherical partition of bulk water not yet ‘injected’ into the vapour with 
the introduction of these contributions and by expressing the result in 
terms of portion functions, the size distribution according to Lothe and 
Pound become 


^ r^i ^ 


(1.33) 


Where Zix and Zir are the translational and rotational partition 
functions of i-mer, Z\.k^ is the i-mer replacement partition function 
[« exp(S//c)], where S is the entropy per water molecules in the bulk 
state (Lothe and Pound, 1962). Nj is the uncorrected classical size 
distribution. Unfortunately, the correction factor 0ip turns out to be of the 
order of 10 , annihilating the classical theory and experiment. 

With the inclusion of this contribution the theory predicts a critical 
supersaturation ratio of about 3.0 per water vapour at 250®K in poor 
quantitative agreement with Powell’s (1928) observations. This situation 
stimulated a great deal of experimental and theoretical work in the field 
on the experimental side (Franck, 1956), two new techniques were 
developed to measure critical supersaturation ratio for homogeneous 
nucleation; 

1. The Diffusion cloud chamber (Katz, 1970; Frenkel, 1946) 

2. The Supersonic Nozzle method (Wegner, 1977) 


Work is still in progress, but it appears that virtually all substances 
with the possible exception of Argon (Hoare, 1980) and Ammonia 
(Dawson, 1969) follow the original macroscopic theory, which ignores 
the contribution from free translation and rotation. Since the innovative 
work of Wilson (1897) the expansion cloud chamber has been developed 
and used by numerous investigators to study homogeneous nucleation. 
Miller et al. (1983) presented experimental homogenous nucleation rate 
data extending over an exceptionally large Jw-Sv.w-T (Nucleation rate- 
supersaturation ratio- temp.) surface. Using the empirical nucleation rate 
formula (Miller et al., 1983). 

; = Sl^exp [328.124 - 5.58243T H- 0.03036572 - 5.0319 x 
10-572 - (999.814 - 4.100877 + 3.01084 x 10-272)^-25^^ (1.34) 

And adjusting only the sticking coefficient a good fit was obtained 
between classical nucleation theory and the experimental data (Hagen, 
1984) .Not surprisingly, the refinement (Lothe and Pound, 1969, 1962, 
1966, 1968) to classical approach has been questioned (Kikuchi, 1971, 
1969; Reiss, 1967, 1970) on theoretical side, much effort has been 
directed towards ascertaining the magnitude of another correction to the 
macroscopic theory, which tends to counter balance the effects of free 
translation and rotation. This is termed as replacement partition function, 
and it describes the free energy due to the six internal degrees of freedom. 
Efforts to calculate this quantity by classical phase internal methods 
(Reiss, 1977; Kikuchi, 1977; Nishioka, 1977) have not been entirely 
successful because of great difficulties in defining a cluster embedded in 
liquid. In fact, some treatments yields larger (Reiss, 1977; Kikuchi, 1977) 



replacement partition functions while others give smaller ones (Nishioka, 
1977). In summary one might conclude at the present time that the good 
agreement between the original macroscopic theory and experiments in 
liquids is fortuitous, but many appear genuine because of 

1 . A large replacement partition function, or 

2. The observation of potential energy due to macroscopic surface 
tension (Abraham, 1974; Binder, 1980). 

Kikuchi (1971) extended and refined the work of Reiss (1968) to 
obtain the following correction factor in place of 0^ p, 

0RK = Pw<^v/kT (1.35) 

Where 6^ is the volume per molecule in the vapour phase. 

At 303k, 0p;^=3xlO'‘. Although, this is more than 12 order of 
magnitude smaller than 0£, p, it still constitutes a significant correction to 
classical results. 

Unfortunately, the studies of Lothe and Pound (1969, 1962, 1966, 
1968) and of Reiss (1977), Kikuchi (1977) and Reiss et al. (1968) suffer 
from fairly obvious short coming ; either the energy of formation of 
germ, or the correction factor itself, or both have to be evaluated, in terms 
of macroscopic values of parameters like surface tension and 
density p^_ This tends to seriously compromise the microscopic point of 
view, which is taken when formulating modification to the classical 
model. Also, there are still partially unresolved conceptual problems, 
which arise in the attempt to merge the classical stationary drop model 
and statistical mechanics points of view. 



1.18 (ii) Computer Simulation and Monte Carlo Method: One might 
think that the issue could be settled by computer calculation using fairly 
realistic potential functions to calculate the free energy of cluster. Much 
work has been done in this area (Hale, 1982; Sjolander, 1978; Ward, 
1982; Abraham, 1974). In one of the more elaborate efforts (Miyazaki, 
1977), Monte Carlo methods were used to calculate the surface tension 
and chemical potential of bulk argon using Lehnard- Jones potential, and 
these results were applied to estimate the free energy of isolated argon 
droplets. The proposed approach is logical development of the statistical 
ions originally suggested for a single layer cluster (Thomas Ackermn, 
2005). 

The resulting free energies of formation were then compared with 
the actual Monte Carlo free energies of isolated argon clusters. It was 
found that the contributions from free translation and rotation were 
required to describe the Monte Carlo Cluster data. This means that, unlike 
the case of physical experiment on real liquids, the original macroscopic 
theory must be modified by the contribution from free translation and 
rotation in order to describe argon by L-J Potential. However, a major 
short coming of this and other computer calculations is that the pair wise 
potentials were used. This is a serious difficulty in view of the importance 
of three body potentials in describing surface properties. Further, a Monte 
Carlo procedure has been used to study the stability and structure of small 
water clusters absorbed on model basal and Prism faces of Hexagonal 
Agl (Ward, 1982) and to estimate the critical cluster size to calculate 
cluster free energies (Hale, 1982). 



1.18 (iii) Molecular Model: For many years the pre nucleated embryos 
have been treated classically (Becker, 1935) i.e. as liquid drop or as a 
solid spherical caps rather than as a collection of water molecules inspite 
of the fact that these embryos consists of typically from two to less than a 
hundred molecules clearly. Surface tension and density of Liquid are not 
applicable in the description of such small embryos or clusters. Strongly, 
however at least in the case of homogeneous nucleation of liquid water 
from vapour, this extrapolated macroscopic model has provided 
qualitative agreement with the experiment (Allen, 1969; Russell, 1969; 
Lothe and Pound, 1969). The liquid drop model also has other 
difficulties, which should be mentioned. There are terms missing from the 
free energy of formation which questionable are supposed to be present 
(Frenkel, 1955; Lothe and Pound, 1962), but have been cancelled by so 
called ‘replacement term’. There have been great deals of illuminating 
work directed towards evaluating the replacement terms (Dunning, 1965; 
Lothe, 1971; Lothe and Pound, 1962, 1966, 1968; Kikuchi, 1969; Reiss, 
1967, 1970, 1968) and these contributions have done a great deal to aid 
our understanding of nucleation process. Molecular model proposed by 
Plummer and Hale (1972) had several advantage such as to predict. 

1 . Molecular properties such as binding energy and external energy of 
the water cluster stored in the form of inter and intermolecular 
vibrations. 

2. The concentration of different sized water cluster or ice clusters. 

3. The detailed information directly related to physically meaning foil 
cluster properties such as cluster size, structure and stability. 

4. It can serve as a frame work for the incorporation of molecules of 

___ 


other substance. 


The basic assumptions of the Molecular model are as follows; 

1 . The cluster have a well defined structure, 

2. The internal properties of the molecules afe only slightly altered by 
cluster formation. 

3. Tne life time of clusters is sufficiently large that its vibrational 
spectrum can be characterized, and. 

4. The most Stable clusters are approximately spherical in shape. 

The assumed structure for water clusters in supersaturated vapour 
is that of a closed or partially closed clathrates composed of five 
membered rings. An example of such a structure for a 20-mer is show in 
fig.(1.27). These cluster forms fulfill the imposed criteria that the 
molecules associated by hydrogen bonding with bond angles, which are 
roughly tetrahedral, that the number of bonds be maximized, and the 
forms posses near spherical symmetry. This choice of geometry is 
supported by the studies of Searchy and Fenn (1974). 

However, these perfectly ordered clatharate structures could not be 
used to represent arbitrarily large i-mers, since it becomes difficult to 
maintain the closed ‘cages’ without grossly distorting the bond angles and 
lengths. This results in the occurrence of considerable bond strain for 
i > 80, the effect of which has been studied by Hagen (1973).However, 
nucleation of ice is not a homogeneous nucleation process. Hale and 
Plummer (1974b) have extended the molecular model to include 
configurational entropy contributions to the partitions function for the 
study of prenucleation embryos of ice in vapour. They assumed an ice-Ih 
structure composed of rings containing six water molecules each. A 
typical structure with 20 molecules is shown in fig. (1 .28). 
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fig. 1.28 : The 20 molecule ice Ih structure. From Hale and Plummer (1974 b) 
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fig. 1.29 ; The free energy of formation versus cluster size at Sv.w. = 5.0 and 
T=273 k for (i) The liquid drop model, AG/kT=9.3 1.6 

n^; (ii) The molecular model, A(j)/kT;(iii) Molecular model with 
strain and (iv) The least square fit to A(j)/kT=12.8nw^^ -2.4nw, 
From Hale and Plummer ( 1 974 a) . 



A comparison was made with the clatherate model and with a 
simple spherical cluster model using ice surface tension to evaluate 
energy of formation and to use it to calculate the nucleation rate and the 
exact critical supersaturation ratio. In this model, the Hamiltonian of the 
n- molecules cluster can be written as 

^(n) = + Hb (1.36) 

Where Ht and Hr , are the rigid body translational and rotational 
energy operators for the cluster. Hy, the vibrational energy internal to the 
cluster and Hb , inter-molecular binding energy. 

The partition function is given by 

= Y)-< ^/^n/kT >] (1,37) 

i 

Where the summation ‘i’ is overall possible states of the system k 
being Boltzmann constant using eq.(l .36), we can write- 

^(n) — (1.38) 

Where Zc , represents the contribution from configuration entropy. 

A comparison of {A(pjkT) (molecular model) and (A(l)/kT) (liquid 
drop model) was made (Hale, 1974a).It was found that at supersaturation 
ratio = 1, the energies of formation differed by about 5 to 10 kT (or 
0.1 to 0.2 ev at 273*^ K) over a wide range of temperatures. The molecular 
model predicts consistently higher energies of formation than the liquid 
drop model. This is in contrast to the corrected liquid drop model with 
replacement terms, which give consistently liquid drop model, fig. (1.29) 
shows the least square fit to {A4>/kT) for ^ — 5-0 and temp. 273° K. 



1.18 (iv) Semi-Molecular Model: The Cluster model (Chan, 1980a), 
with ion at the centre features two separate regions and is, therefore, 
discontinuous. It consists of a hydration shell containing discrete water 
molecules. The molecules are strongly bonded to the ion and their 
orientation is determined by the ion. Statistical mechanics is applied to 
calculate the excess free energy of the micro clusters. In the outer region, 
the water is taken to be a dielectric continuum and the modified Fisher 
droplet model for ion-water nucleation is used. 

This model has the advantage that it can be applied to both ends of 
the ions cluster spectrum. It predicts saturation ratio for condensation on 
ions close to values from cloud chamber experiments. A sufficient part of 
molecular approach has been included so that the model makes a 
distinction between positive and negative ions. 

1.18 (v) Extension of Molecular Model to Heterogeneous Nucleation: 
It is difficult to extend the molecular model to a study of heterogeneous 
nucleation. This is because very less is known about the mode and energy 
of interactions between clusters- of water molecules and substrate 
molecules. The major problem in determining the interaction energy is 
that it is specific to each solid and to each crystallographic arrangement. 
Further more it is a function of position on each crystallographic face, 
due to the heterogeneous nature of the surface. Inspite of these 
difficulties, preliminary studies (Plummer, 1973) have produced at least 
qualitative agreement with observation by predicting that the weaker the 
interaction between the surface cluster and solid substrate, the more likely 
it is that the water molecules will arrange in isolated multi-layered 
structures instead of a single layer. 



The conceptual difficulties associated with assigning surface 
tension and contact angles to microscopic embryos are avoided with a 
molecular model. The affinity of the surface for water and the strength of 
the water substrate interaction can be specifically included in this model. 
The model thus has the potential to predict which artificial nuclei should 
be the most effective in producing ice crystals or cloud droplets. The 
effect of temperature and supersaturation on nuclei activation can also be 
studied. The role of surface imperfection and contaminants can well be 
included. 

1.19 MODELS FOR HETEROGENEOUS WATER NUCLEATION: 

The classical model of nucleation theory is used for heterogeneous 
nucleation. The “replacement term” difficulties associated with the 
attempt to define the classical approach do not arise, since the nucleating 
substrate may be taken to be at rest. The modification has been made 
according to the geometry of surface. 

1.19 (i) Nucleation on water in soluble, partially wettables CCN: It is 

assumed that water embryos or germ or nucleus, nucleated in 
supersaturated vapour on a water insoluble partially wettable surface 
assumes the shape of a cap. 

1. Nucleation on a planar Substrate; Gibb’s free energy for the 
formation of spherical embryos is given by 


AGj j ~ ( 3 ^^vol "b ^w/v^ ^w/u) (1-39) 



Where ai is the radius of embryo, the surface tension of the 
water and = cos6, 6 being the angle of contact of the 
surface with water. 

2. Nucleation on a curved substrate: Fletcher (1958, 1959a) has 
extended the theory to include an amount of the effects of the 
finite size of the nucleating particle. In this extension the substrate 
is assumed to be a sphere of radius r^.. ' 

Accordingly, the energy of germ formation becomes 




AGas= -r 

3[RTpw-e 




(1.40) 


where 


2f{m,X) = 1 + ( 
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( 1 . 41 ) 

W[th0 = (l-2mX + X^) 2 , X = rN/a (1.42) 

Twomey (1964) modified the heterogeneous nucleation theory of 
Fletcher to consider nucleation on an insoluble, partially wettable 
spherically concave substrate of nucleus tn. The only change in theory is 
one of geometry. The volume energy in this case is given by changed 


value of 


2f(m,X) = 1- 


1 + mX 


X ^\2 


X + m\ fX + m 


+3mX^ - l) (1.43) 


and 0 = (1 + 2mX + X^y, X = 


(1.44) 


Mahata’s computations indicate that concave surface features can 
sufficiently enhance the capacity AP to serve as CCN. 


1 .20 £L£CTE1C FIELD AND ELECTRIFICATION OF CLOUDS: 


Here we shall describe various types of clouds and cloud particle 
acquire their characteristic state of electrification. 

1.20 (i) Electric State of cloudless Atmosphere; Under clear sky 
conditions, flat portions of the conducting earth carry a negative surface 
energy density which is approximately (To =-3.4x10'^ e.s.u. cm'l By, 
Gauss’s law there is a corresponding downward directed surface electric 
field of magnitude Eq = 4nao(.e.s.u) = 3 x lO'^Vm"^. The field 
strength in the atmosphere varies with height in accordance to Coulomb’s 
law (i.e. as d'^) where ‘d’, denotes distance from centre of the earth. 
However, in the rate of attenuation is much greater than this due to the 
existence of positive space charge, which rapidly screens out the surface 
field with increasing height. 

An approximate empirical description of the variation of the 
electric field strength with height is 

E(Vm-i) = 81 . 8e-^-S2z + + 10 . 27 e-®-i 2 iz (i 45) 

Where z is in Km. from Gauss’s law, 

V.E = - (1.46) 

£0 

The corresponding positive space charge density p (elementary 
charges cm'^) is 

p = 20.4e-^'^2z ^ 0.8e-°-372 + 0.069e-°-^2iz 47 ) 

Thus at ground level the charge density is about 21e cm'^, while the 
average over the first Kilometer of the atmosphere is about 5e cm’^. 


dl) 
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The small ions produced by cosmic rays and radioactive substances 
are generally singly charged molecules. They may become large ions by 
attaching themselves too much larger aerosol particles. Since the 
mobilities of large ions is usually 10 ’^ to 10 '"* times less than those of 
small ions, they contribute relatively little to the total ion current. 
Measurements show that the represented mobilities of the small 
negative ions are greater than the corresponding mobilities B+ of the 
positive ions according to Bricard (1969). 

B+ = 1.4cm^v”^sec“^ and B_ = 1.9cm^v”^sec“^ at STP. 

An extrapolation of these values to elevations Z< 10km in the 
standard atmosphere has been computed (Shreve, 1970) 

B+ = 1.4e0•l4^ B_ = (1.48) 

Where z is in km. Shreve also obtained expressions for the 
variation with altitude of the corresponding diffusivities (in cm'^sec'*) for 
example 

D+ = 3.6 X 10“2e°°^2z^ ^ 4 3 ^ 49 ) 

The vales for Z= 0 were obtained by substituting the corresponding 
mobilities and T=273K.In recent years the propagation effects on the 
electric field time derivative generated by return strokes in lightning 
flashes has been studied by Mahendra Fernando, Vernon Cooray (2007). 

1.20 (ii) Electrical state of the Atmospheric Aerosol: Aerosol particles 
acquire charge through Brownian deposition of ions. In term, such 

CEJ 


charged particles experience Brownian coagulation, which is enhanced or 
suppressed by electrostatic forces. Under fair weather conditions aerosol 
charging is approximately symmetrical, since small ions are created in air 
and have roughly equal mobilities. Therefore, in equilibrium particles of 
any size have a charge distribution approximately symmetrical about 
Zero charge, so that for every Vj particle bearing charge qj. Besides 
acquiring charge through ionic diffusion, aerosol particles polarized in the 
fair weather electric field and thereby received an additional ionic drift 
current. 

According to Brian A. Tinsleya (2005) et al. the electrical charges 
on aerosol particles and droplets modify the droplet-particle collision 
efficiencies involved in scaveging, and the droplet- droplet and particle- 
particle collision efficiencies involved in coalescence of droplets and 
particles, even in only weakly electrified clouds and aerosol layer. 

1.20 (iii) Cloud Electrification (Models for cloud electrification): 
There are a number of models for electrification of clouds. A few of them 
are given below; 

1. Convection charging 

2. Particle charging by selective ion capture 

3. Particle charging by thermoelectric effect. 

4. Particle charging by induction. 

5. Non-inductive charge transfer process. 

6. Electrification of thundercloud by an entrainment mechanism of 
charge separation. 

7. Electrochemical charge separation in clouds. 



1. Convection Charging: Vonnegut (1955) independently 
proposed that a convective cloud may operate as an electrostatic energy 
generator according to the following scenario: Initially, an up draft carries 
positive space charge from the lowest levels of the troposphere into the 
growing cloud, the electric field cloud soon acquires a negative charge 
screening layer at its edges due to cloud, the electric field cloud soon 
acquires a negative charge screening layer at its edges due to cloud 
particle capture of ions drifting from clear air to cloud under the influence 
of the main positive charge, finally, downdrafts carry the negative charge 
closer to the ground, thereby increasing the (reversed) electric field at the 
earth sufficiently to initiate positive point discharge which enhances the 
positive charge entering the cloud via the updraft. 

The time dependent model for cloud electrification has been 
developed (Mathpal, 1982) in which both convection and gravitational 
separation of charged precipitation particles work simultaneously towards 
its electrification. Not only the generated electric field but also the current 
density have been calculated for various values of precipitation intensity 
Po, liquid water content L, rebound mean angle probability <P>, 
relaxation time x, fractional constants of space charge and electric field 
f]f 2 and the vertical air motion U. It has been found that an electric field 
of about 4.0 X 10^ Vm"' and a current density of about 0.1 pAm“^ may 
normally be achieved for a vertical air motion U « 4ms“^, at a 
precipitation intensity of about lOmmh"’. However, occurrence of the first 
lightning flashes may also be explained even at precipitation intensity 
5mmh'^ ( U « 8ms“^ ). Thus it appears that the thundercloud 
electrification may be explained by the combined precipitative and 



convective (CPC) charging mechanisms and the vertical air motion in the 
cloud plays an important role in its electrification. The role of parallel 
electric fields in the mass-dependent effects in ion conic production had 
been studied by E.J. Lund et al.(1999). 

2. Particle charging by selecting ion Capture: Wilson (1929) 
described how an electrically polarized cloud particle may selectively 
capture ions of one sign as it falls. This happens, because while the lower 
surface of the particle may easily attract and capture ions which carry a 
sign opposite to the local surface charge, the upper surface is not as 
effective in this respect since ions attracted to it must first catch up with it 
in order to be captured. The net effect of this selective process is a large- 
scale separation of charge from due to the sedimentation of the charged 
cloud particles. This reinforces the existing field, so that its occurrence in 
the clouds would cause a field enhancement in qualitative agreement with 
what is expected for thunderstorm. 

Mathematical model for the “Wilson Process”, based on spherical 
particles in stokes flow has been worked out by Whipple and Chalmers 
(1944). The equilibrium charge for this process is proportional to the 
ambient field strength and particle surface area. Also, the process 
becomes ineffective for large fields, since then the ionic drift velocity 
greatly exceed drop terminal velocities. A simple estimate (Mason, 1971) 
suggests the mechanism may produce fields only as large as 100 Vcm"’, 
which is still three orders of magnitude below characteristics 
thunderstorm values. 


3. Particle charging by thermoelectric effects; Reynolds et al. 
(1957) obtained laboratory evidence that a Hail Pellet may become 
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charged as a result of collisions with ice crystals having a temperature 
different from that of the pellet. The physical basis of the charge transfer 
was suggested as being due to the diffusion of hydrogen ions down the 
temperature gradient existing in the region of momentary contact 
(Brook, 1961). Thus, Since ions have a greater mobility in the ice 
lattice than 0H“ ions, a temperature gradient maintained across a piece 
office will result in an excess of positive charge on the colder portion. 
Mason formulated a one-dimensional model for this process on the bases 
of an ideal ice structure and found that magnitude of equivalent surface 
charge density 0 on the ends of an ice rod having prescribed steady state 

temperature gradient dT/dx (®C. cm"’) is given approximately by 
dT 

cr « 5 X 10~^— e. s.u. cm“^ . Similar results follow from a more 

dx 

complete theory of the process in real ice having orientational defects 
worked out by Jaccard (1963). Apparently the results depend sensitivity 
on such factors as the time of contact, micro-topography of the areas of 
contact, presence of impurities in the ice, and the relative velocity 
impact. All of these factors are difficult to control in an experimental 
situation and hence their effects are not easily evaluated. 

4. Particle charging by induction: An uncharged drop polarized in 
the fair-weather electric field carries a surface charge density. 

ae. s. u = (4Tr)“^(Er)j.=a “ 3ECos0/4'n: (1.50) 

0, being the polar angle measured from the lowest point on the 
drop; thus the lower hemisphere is positively charged, and the upper 
hemisphere negatively charged. Therefore, if such a drop were to 
experience on its lower hemisphere a momentarily electrical contact with 
and subsequent separation from a similarly polarized smaller drop, there 


would result a net negative charge on the larger drop and a positive 
charge of equal magnitude on the smaller drop. The first to point out that 
such a process of inductive charge transfer, occurring throughout a cloud 
and followed by the large scale separation of charge through relative 
sedimentation under gravity, would serve to increase the in-cloud electric 
field in the sense normally observed in thunderstorms. 

5. Non-inductive charge transfer process: These processes occur 
during particle collisions in which charge is transferred independently of 
the local electric field strength. It is perhaps unusual for a process to be 
named in terms of what it is not however, a more succinct and yet all 
encompassing little has yet to emerge. Beared and Ochs (1986) use the 
term “interface charging” for charge transfers driven by differences in 
contract potential. 

Detailed laboratory Studies of the above charging processes have 
been carried out over the last ten year. Jayaratne et al. (1983) confirmed 
that the sign of graupel charging reverses as a function of temperature and 
that this “reversal temperature”, moves to higher (warmer) values with 
decreased cloud liquid water content. They also found that the amount of 
charge transferred depends on impact velocity and on the size of the small 
ice crystals used. 

Keith and Saunders (1990) extended the previous work by using 
larger ice crystals, up to 800pm diameter, and found that the charge 
transfer increases rapidly with crystal size for small crystals but increases 
at a lower rate for larger crystals. They suggested that high values of 
charge transfer are limited by the reverse charge transfer of some of the 
charge residing on the surfaces when the particle in the dark. A study of 



THESIS OF DEVESH DWIVEDI 


crystal-graupel interactions in the dark revealed light emission associated 
with this reverse charge transfer in the form of Corona (Keith, 
1988).From their charge transfer experiments; they formulated 
relationships between charge transfer, crystal size, and impact velocity 
for positive and negative charging situations. 

6. Electrification of thundercloud by an entrainment mechanism 
of charge Separation: Wagner and Telford (1981) proposed an 
entrainment mechanism for cloud electrification in which the micro-scale 
charge separation can take place at the interface. This conductivity in 
homogeneity has been postulated to result either from the one set of ice 
process or from the evapourative depletion of small ice or water particles 
in cloud particles, diluted by dry air entrained at the cloud top. They 
consider that the turbulence inside the thermals and the prevailing 
updrafts and downdrafts within the cloud are responsible for macro- 
physical charge separation. The ion pairs produced due to cosmic rays, as 
in the convective theory, provide the convective space charge in this 
mechanism. Sapkota and Varshneya (1988) examined the extent to which 
the leakage current produced by the electrical forces acting on the 
descending thermals influences the net charge separation caused by the 
entrainment mechanism. P.T. Tonev, a. and P.I.Y. Velinova (2007) 
examined Atmosphere-ionosphere vertical electric coupling about 
thunderstorms of different intensity. 

7. Electrochemical charge Separation in clouds: Pathak et 
al.(1980) presented a self-sustaining mechanism for the continuous and 
rapid lightning flashes. It has been shown that the initial discharge 
provides sufficient number of ions for the electrochemical charging 
mechanism proposed (Wahlin, 1973) to be operative for the rapid 
separation of charges within a thundercloud. 
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The mechanism is based on experimental observations by Wahlin 
and uses electrophillic property of NO 2 and O 2 molecule present in the 
air, when any surface is ventilated to ionized air, these molecules with 
electrons attached to them react with surface material transferring 
negative charge to it. The positive ions are left in surrounding air. They 
form an electrical double layer with charged surface at a distance equal to 
the average distance between ions in air (about 1mm). 

The Coulomb force does not bring them closer to oppositely 
charged surface. The surface is negatively charged because negative NO 2 
and O 2 ions can react with surface material to transfer negative charged 
and thus a net positive charge is swept with air. For a spherical surface of 
radius R the maximum charge that can be separated at equilibrium is 
given by 

qmax = Vo47I6oR[1 + (R/d)] (1.51) 

Where Vo, is the oxidation reduction potential of the surface 
material; d, the thickness of electrical double layer. 

1.21 EFFECT OF ELECTRIC FIELD OF MICROPfiV'SlCAL 
PROCESS: 

1.21 (i) Drop and Ice Crystal Nucleation: In the pr esent work w e have 

V”-* "" 

examine a number of physical parameters such as Gibb’s free energy, 
critical radius, equilibrium concentrations and nucleation rate, related 
with nucleation of water drop and ice crystal affected by the external 
electric field. 
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1 .22 CONDUCTION IN THE ATMOSPHERE: 

The surface of the Earth and highly conduction layer of ionosphere 
may be treated as the plates of spherical condenser. Ions caused b)^ 
various ionizing agencies e.g., Cosmic rays radio activity in air and 
ground etc. exist between the above mentioned plates. Possible influence 
of cosmic rays on climate through thunderstorm clouds were studied by 
LevT. Dormana, B, and Irena V. Dormanc (2005). With decreasing of 
cosmic rays intensity caused by increasing of solar activity or in some 
short periods of Forbush decreases, the intensity of secondary CR 
relativistic electrons decreases and the probability of formation of thunder 
storm is also expected to decrease. Cosmic ray-induced ionization in the 
atmosphere had been studied by LG. Usoskin (2004). 

Electric field changes of lightning observed in thunderstorms had 
been studied (William H. Beasley, 2000). Wilson’s (1920) Thunderstorm 
theory could explain reasonably some of the necessary conditions, such 
as explanation to unitary diurnal variation in potential gradient, for the 
fundamental question in atmosphere electricity. “How the atmosphere 
has a permanent electric field in spite of its conductivity.” The 
conductivity was measured by Gerdein’s Chamber described by Khera 
(1994). 

Gondot et al. (1988) have developed a new type of sensor to 
measure ionic electrical conductivity around and inside convection 
clouds. The resulting technical note is rather complicated and difficult to 
understand but it certain much material with regard the investigation it 
self, an attempt was made to utilize existing stations and network by 
combining their measurements in a sophisticated system which was 
called the “atmospheric electricity”. 
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1 .23 VfeRTlCAL STRUCTUEES OF TROPICAL PRECIPITATIN6 
CLOUD SYSTEMS: 

Krishna Reddy (2003) studied temporal variation of vertical 
structure of precipitating cloud in tropical India based on the observations 
of 1.357 GHz lower atmospheric wind profilers (LAWP) at GADANKI. 
The occurrence of deep convection in the tro| 3 ies plays an important role 
in global circulation, since it transports heat, water vapour and so on to 
the upper troposphere. 

Kandalgaonkar et al. (2003) analysed the data of point discharge 
current (PDC) measurements during a total of 65 thunderstorm at Pune 
(180® 32’N, 73®, 5rE) to study the (PDC) local diurnal variation and 
some issues related with the current. Majority of the studies were carried 
out sporadically during time bound field experiments to accomplish 
specific objectives of the study but the significant issue of the diurnal 
variation of PDC received a little attention (Williams, 1985) except a few 
studies (Whipple, 1936; Ette, 1980). 

In 80’s positive charge deposition by lightening (Holden et al, 
1983), phenomenon for the polarity of the thunderclouds and charge 
reversal microphysics (Jayratne, 1985) have been supported by the 
evidence now available from both field and laboratory. 

The laboratory result on charge reversal microphysics for the 
infrared charge transfer per collision is strongly dependent on crystal size 
over a range from 10-100 microns, and then levels off for larger crystal 
sizes (Marshall, 1988). 


Bourdeau and Chauze (1988) studied numerically two process that 
contributed to the limitations of hydrometer electricity charge in an 
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ambient field Frozen hydrometers the net charge is only limited by 
corona emission. Vonnegut (1988) proposed the mechanism to explain 
thunderstorm electrification differ greatly in how they might affect or be 
affected by the cloud and the clouds environment depending upon which 
mechanisms predominant, atmospheric electrification phenomenon may 
for some cases be in little influenced by environmental factors or exert 
only a minor influence on meteorological process. 

Thorpe (1981) in a discussion of thunderstorm dynamics, “It has 
long been established that although thunder and lightning are the most 
spectacular, they are the least important in the storm structure and 
evolution. This is not to say that understanding of the electrification does 
not represent an enormous challenge, lightning causes a lot of damage, 
but the storm dynamics would almost identical if thunder and lightning 
were not present.” 

During water vapour condensation the equivalence between 
supersaturation ratio and the external electric field has been established 
recently (Dwivedi, D.et al., 2008). 
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MODIFICATION IN THE NUCLEATION OF WATER 
VAPOUR CONDENSATION AND ICE GLACITION, DUE TO 
EXTERNAL ELECTRIC FIELD 


2.1 INTRODUCTION: 

There exists normally an electric field in the clouds. When the 
electric field becomes enough an electrical discharge (the bolt of 
lightning) occurs with clouds or between clouds and the grounds. During 
the strike, successive portion of air becomes a conductive discharge 
channel as the electrons and the ions of air molecules are pulled away 
from each other and forced to flow in opposite direction. During a 
lightning discharge the maximum electric field (~10 e.s.u.) is produced 
near the channel. 

Upward-propagating discharges from thunderstorms called “blue 
jets” and “blue starters” and high altitude near-cloud top discharges have 
been reported over strong convective storm systems (Glenn E.Shaw, 
1998) . The electric field so generated may affect the rate of condensation 
of water vapour. 

It has been reported (Jalaluddin et al., 1972; Parmar et al., 1973; 
Parmar et al., 1975) that the electric field induces nucleation in liquid. 
Schaefer (1953) , Salt (1961) and Pruppacher (1963) observed the 
freezing of water due to passage of electric discharge through it. It has 
been observed that electrical charges on aerosol particle and droplets 
modify the droplet particle collision efficiencies involved in scaveging 
and the droplet-droplet and particle-particle collision efficiencies 
involved in coalescence of droplets and particles, even in only weakly 
electrified clouds and aerosol layers (Brian A.Tinsleya et al., 2005). 

— di 


Poc (1967) , Schaefer (1968) , Garraud (1969), Roulleau et al. 
(1971), demonstrated that the electric field is mainly responsible for 
enhanced ice crystal production with in the supercooled fogs. G. Mandal 
and P . Pradeep Kumar (2002) observed experimentally that ice- 
nucleation does occur inside the cloud if there is an electrical discharge 
inside the cloud. However not all electrical discharges resulted in ice 
nucleation. In no case was ice nucleation observed in the presence of high 
electric fields with no discharge. Varshneya (1969, 1971) and Pruppacher 
(1973) found that when charged ions and aerosols come in contact with 
the neutral drops, the nucleation process is enhanced. Barlett et al. (1963), 
as well as Maybank et al. (1967) reported the tendency of generation of 
few ice splinters during ice crystal growth in strong electric field. 

Loeb (1963) discussed the disruptive effect of an electric field and 
suggested that the nucleation of the fragments of a drop is responsible for 
the observed enhancement in nucleation. Later Abbas and Latham (1969) 
gave a possible mechanism. Doolittle and Vali (1975) studied the effect 
of an electric field on sample of water at different time temperature 
cycles. No intrinsic effect was observed up to 6000 V cm-1. Evans (1973) 
experimentally demonstrated the effect of an electric field on the 
production of ice crystals in cloud chambers and argued that the 
accelerated charged water molecules moves to crystal tips, thereby 
increasing the nucleation rates. 

The effect of electric field on condensation and further 
polaraizability of water-vapour molecules have been discussed by Singh, 
N., et al. (2000, 2001). A calculation of electric and magnetic fields at the 
Earth’s Surface due to an ionospheric electrojet system is done by Risto J. 
Pirijola(1998). 
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Murino (1979) studied theoretically the effect of an electric field on 
the condensation of water vapour and concluded that under similar 
temperature conditions a bigger size of drop can be produced in a given 
time than the one obtained in the absence of an electric field. He observed 
the polarization of water vapour molecules in the electric field of the 
central dipole (embryo of water) alone. 

Singh, N., et al. (1986) discussed the effect of external electric field 
on the relaxation time in the nucleation process of water vapour 
condensation and ice glaciations. The similar study has been extended 
(Sharma, A.R., et al.,1992) and concluded that a small value of electric 
field is equivalent to very high supersaturation ratio to get a nucleus of 
given size under similar conditions of temperature. 

The effect of electric field is dominant over the Indian 
subcontinent (Manohar, G.K., et al., 1999). Shaw (1998) has been shown 
the effect of electrical discharge in the nucleation of water droplets. 
Kumar, A., et al. (1998) discussed the effect of orographic features on 
atmospheric electrical parameters of different cities of India. Recently 
Dwivedi D., et al. (2008) has established equivalence between the 
external electric field and supersaturation ratio during water vapour 
condensation. 

According to Pinksy et al. (1999), the inertia of cloud droplets 
plays a crucial role in the evaluation of the spectral during drop growth 
by diffusion. But, the electrical analogy does not apply to modeling dry 
deposition of particle (Venkatran, 1999). The fair weather electric field 
occurring soon after sunrise has been observed (Marshall et al., 1999). 
The changes in electric field due to lightening in thimderstorms have been 
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observed by Beasley et al. (2000). The most important features of quasi- 
electrostatic fields and currents, generated in the region between a 
thunderstorm and the ionosphere between lightning discharges, are 
theoretically investigated by P.T. Tonev et al. (2007). The electric field 
alignment of ice crystals in thunderstorms has been studied by Saunders 
et al. (1999). Although lightning flashes in thunderstorms are among the 
most ancient natural phenomenon known by human kind, several 
important physical processes related to them are not understood in detail 
yet (Lyons et al., 1999; Mac Gorman et al., 1998). In order to contribute 
to this understanding, recent effort has been done, for instance; the 
studies of the effects of thunderstorm and lightning on upper atmosphere 
(Sentman et al.,1998), and the development of a 3 dimensional detection 
of lightning path, as the recent methodology using GPS-based mapping 
system (Krehbiel et al.,2000). Electrical conductivity effects on lightning 
path simulation has been studied by Odim Mendes Jr. et al. (2005).The 
formulation of relaxation time required for the attainment of the steady 
state concentration of embryos of the critical size has been discussed 
previously(Kantrowitz, 1951 ;Probstein, 1 95 1 ; W akeshima, 1 954;Collins, 
1955) . Following Kinetic theory (Collins, 1955) inferred that the 
relaxation time is independent of the free energy of formation of the 
nucleus, but it varie& as the square of the radius of the critical nucleus. 

2.2 THEORETICAL CONSIDERATION: 

In the nucleation process of water vapour condensation and ice 
glaciations the fundamental quantity of importance is the Helmholtz or 
Gibb’s free energy of germ formation. But the relaxation time is found to 
be independent of the Gibb’s free energy. However, at normal electric 
fields the electrostatic energy of germ formation is very small compared 
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to the volume and surface energy terms. It becomes appreciable only at 
high electric fields at which break down in air might occur. Also the 
critical size of nucleus is produced quickly. Therefore, in this chapter we 
consider only the relaxation time of germ formation and thereby study the 
effect of an electric tield on the process of self-nucleation (homogenous 
case) only. While calculating the relaxation time of water phase we have 
taken the variation of surface energy and density of water according to 
the different temperatures given by Edward, W.W et al, (2003) as 
following. 


Tabel-2.1(A) 


T (Kelvin) 

Pwigrti/cTV?) 

<Tw/v(.(^ynes/crn) 

273 

0.999841 

75.7 

278 

0.999965 

74.9 

283 

0.999700 

74.2 

288 

0.999099 

73.6 

293 

0.998203 

72.8 

298 

0.997044 

72.0 

303 

0.995646 

1\2 


Similarly for ice phase Smitha R. P. et al. (1986) observed from the 
Antonow’s formula that for ice phase the surface energy of ice vapour 
interphase also changes according to the temperatures. 

^i/v ~ ^ijw ^w/v 








In which the value of {oii^ = 33 dynes/cm) which is independent 
of temperature. Thus in calculating for ice phase we will add 
<^i/w the surface energy of water vapour. Now the values of surface 
energy of ice phase will be as following: 

Table-2.1 (B) 



2.2.1 Homogeneous Nucleation Process: 

2.2.1 (i) In The Absence of An Electric Field: 


For the formation of a new phase a substantial supersaturation of 
the mother phase is required to provide the nuclei necessary for initiating 
the growth of this phase. In some very rapid process involving the 
formation of aerosols. There are indications that at supersaturation 
sufficient to provide nucleation according to steady state theories, there 
may be a measurement delay in nucleation. Probstein (1951), Wakeshima 
(1954) and Kantrowitz (1951) have shown that such delays in nucleation 
may be anticipated on the basis of the relaxation time required for the 
attainment of the steady state concentration of embryos of various sizes. 
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The frequency of growth of embryos beyond a critical size in the 
nucleation process corresponds to a free energy maximum. This energy is 

^w/v^n (2T) 

Where [ 1 ^;, is the chemical potential of vapour and [ii, chemical 
potential of liquid. (Tw/v> the surface free energy of water vapour 
interface; and surface area of embryo; ny^^, is the number of water 
molecules in the embryo. 

Corresponding to the free energy maximum for a critical size 
nucleus, we have, 


and thus the eq. (2.1) becomes 


A/-* _ ^q-w/p^w , 

~ „ ,1/3 ^ 

Sriw 


^w/vl^w^w 


( 2 . 2 ) 


Where = (2,3) 

rriw, being the molecular mass and the density of water. 

The growth of the embryo involves the step wise gain or loss of 
single molecules so that the embryo undergoes an analog of a random 
walk among its various sizes. Similar to the classical diffusion equation, 
the transformation equation of the set of differential equations 
representing this process is. 


dt dTivt 




+ J_(dc 

driw \ / 
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Where C = C (n^^, t) is the concentration of embryos of the size 
at the time t, and the (Diffusion Constant) D= D(n) is 

^ (2.5) 

Wher e is the frequency of collision of single molecules per unit 
area. In the case of self nucleation. 

Pw = 'n-i (kT ( 2 . 6 ) 

Where rij is the concentration of water vapour molecules; k, the 
Boltzmann constant; T, the temperature of the system. 

Case I - Steady State Case: For the steady state case eq, (2.4) 
becomes: 



With the assumption that the concentration of single molecules is 
constant and that for some embryo containing molecules, where,n(v>ni^, 
the probability of decomposition is zero. Hence the boundary condition 
are: 

C(l) = n, :C(n;) = 0 (2.8) 

The flow of embryos past the critical size is 

/ = (2.9) 

Where 




2 . /3mvp '\ *- 4/3 


f 
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Case II — Non Steady State Case: To calculate the mi nimum limit 
of relaxation time in the non-steady- state case, a constant value D* has 
been assumed disregarding its dependent on n^. 

The equation (2.4) then becomes 


^ = D-££ + D--i. 

dt dnlj dUu 


(c 

\ J 


( 2 . 10 ) 


The solution of this equation gives the time-dependent flux for 
small values of t. 
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The relaxation time may be written as 


T = InkTIyD* (2.12) 

The relaxation time is the characteristic time to achieve the quasi- 
steady-state (germ concentration equilibrium). In the absence of an 
electric field Collins (1955) evaluated the relaxation time. 

To = 9nkT(nwy^^/[iwl3w(Jw/v (2-13) 


The number of water molecules in a critical nucleus is given by; 


n; 


W 


= ( — ie — y 

Vl.958xl0-®/ 


(2.14) 


Where is the radius of the critical nucleus (Pruppacher and 
Klett, 1978) with 




'w 


Ppj; RT £tl 


(2.15) 






In this expression ^ is the supersaturation ratio of vapour over 
the plane water surface. 


2.2.1 (ii) In The Presence of An External Electric Field: 

Water is a strongly polarizable medium with an electric dipole 
moment of 1.83 xiO''* e.s.u. cm. (Eisenberg, 1969; Hobb’s, 1974). A 
perturbation by an external electric field induces an electric dipole 
moment on the embryo of water as well as on the surrounding water 
vapour molecules. 

The moment induced on the embryo is 

M = ErJ (2.16) 

Where is the radius of the water embryo and E is the inducing 
electric field. The moment induced on a water vapour molecules (Kittel, 
1966) is given by 

(2.17) 

Where (k=: 5 x 10“^^cm^) is the polarizability the electric field 
generated at any point (A) due to the embryo dipole, fig. (2.1), is 

Ej_ = [\3{M .X)X -X^M\X^] (2.18) 

Where X is the radius vector joining the centre of an embryo and 
the water vapour molecules, shown in fig. (2.1). 

However, the vapour molecules experience the external electric 
field as well as the field produced due to the embryo dipole. 
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Taking into account the total electric field (E + E^), the induced 
moment on a vapour molecule can be -written as 

M{=oc(E + E^) (2.19) 

Using eq. (2.18) and eq. (2.19), the potential energy of a vapour 
molecule becomes. 

U = 

= - a [dCos^e - l][(E^r^X^) + (E^r^X^)] (2.20) 

Where 9 is the direction of point (A) with respect to the direction 
of the electric field. The molecules can have any orientation ranging from 
0 to 2n. For simplicity we are replacing {3Cos^9 — 1) by its average 
value 

< scos^e - 1 >= — - 1) = i 

2jr “^0 '' 2 

The average potential energy is given by 

U = -f [(FV^/A3) + E^4/X^] (2.21) 

The negative sign in the expression for the potential energy 
shows that the potential is attractive and hence there is a flux of water 
molecules towards the embryo dipole. From the kinetic theory of gases 
(Engel, 1955) it is assumed that in a non-uniform field, the drift velocity 
is given by the transformation of potential energy lost in one mean free 
path into kinetic energy. The loss in potential energy in one mean free 
path is 



( 2 . 22 ) 
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jC*2-j-3 


a+A)3 (X+A)6 
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Where A is the mean free path (~10‘^ cm), from eq. (2.21) and eq. 
(2.22) substituting X = r^, the change in potential energy is given by 


dU = -(9 oc 


(2.23) 


(i.e. the change in potential energy of a molecule in a transition 
through a distance A onto the surface of the embryo of radius r^) 


The corresponding gain in kinetic energy is given by 


K.E.= 


(2.24) 


Where is the mass of one molecule and v its velocity. From eq. 
(2.23) and eq. (2.24) one can write from the law of conser\'^ation of 
energy that 


(9 a AE^/2r^) = 


Therefore 


V = (9 oc 


(2.25) 


The increase in mass (dm^) of the embryo in time (dt) is 


dm^ 


dt 


pySfiV 


(2.26) 


Where is the surface area of the water embryo of radius and 
Pi, is the density of the water vapour molecules, from eq. (2.25) and eq. 
(2.26) we have 
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dmw _ I' Pv\ 
dt XPw^ ym^Tw) 


(2.27) 


Rewriting this in terms of the increase in radius with respect to 
time, we have 


drw _ /M 
dt ypw-' 


(2.28) 


Integrating within the limits = 0 to (cr itical radius of 
nucleus) and t = 0 to t = (relaxation time), we get 


= [3pvi9 cc AE^ (2.29) 

Putting, 2 w 10“^cm,oc= 5 x = 3.0 x 

Pw = 1 ctnd Pt; = 1 0~^ (a t ~10°c) we get 

r* = (3.18 X (2.30) 

From eq.(2.30) the relaxation time in the presence of an electric 
field is, 

Tt, =rf^/3.18xlO-'^E ( 2 . 31 ) 


Thus the relaxation time for the growth of a nucleus in the presence 
of an electric field varies inversely with the applied electric field. Under 
the combined effect of a n electric field and diffijsio n the reduced 
relaxation time To A/ is given by 



‘^ON 


_ — .ZsIiL 

To+Tw 


(2.32) 
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Where Tq is the relaxation time in absence of electric field with 
only diffusion of water vapour molecules. 

2.2.2 Ice Nucleation in the presence of an Electric field: 

In the atmosphere the water vapour changes its phase via either of 
following thermodynamic processes 

Condensation 

1. Water vapour ^ Liquid water. 

deposition 

2. Water vapour > Solid Water (ice) 

Condensation freezing 

3. Water vapour > Liquid water >ice 

The proc ess concerned in the present discussion is of type (if or fiil 
only. For the purpose of type (iii) a more rigorous mathematical treatment 
is required. 

The above derivations are applicable equally well to the ice 
formation from the vapour phase (type iii) by introducing the following 
replacements 

Pw Pv ^w/v ^v.i 

2.3 DISCUSSION OF THE RESULTS: 

In homogeneous nucleation in absence of external electric field 
Collins (1955) evaluated the relaxation time for the growth of a nucleus is 
given by eq. (2.13) 

Tq = ! Pw /v 

— ' 
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In homogeneous nucleation in presence of external electric field 
taken into account the induced dipole by eq.(2.31) 

Tn =r^^^^/3.18xlO-\E 

Murino (1979) also, obtained an expression for the drop growth in 
presence of an electric field without considering dipole contribution. With 
some manipulations the relaxation time on the bases of Murino’s theory 
can be written as 


X 10 


(2.33) 


The relaxation time, in the absence of an electric field, varies 


directly with the rad ius of the nucleus, but in the presence of an electric 
field, it varies . In the absence of an electric field one can write 

from eq. (2.14) 

. n; = (C/1.958 X 10-8)3 
And from eq.(2. 1 5) we get 




,W 


From eq. (2.6) we have 

Where, rii is the concentration of water vapour molecule; k, the 
Boltzmann constant; T, the temperature of the system. Here 


n.of molecule N 

unit volume 22.4X10^ 


(2.34) 
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Putting the values of in eg. (2.31) and putting 

m^= 3xl0-^^gm.; Mw =18; N - 6.023x10^^; k = 1.38xlF^erg K'Vi 
R-8.317xl0^erg K'' mole'*; = 2.1792x10*^®. We have 


To = 3782.026 x 10"^ 




(2.35) 


From eq. (2.35) we can calculate the different values of Tq at 
different a, p, T. 


Similarly the values of Tjv and have been calculated from eq. 
(2.31) and eq. (2.33) respectively. 


The variation of and Tq with supersaturation ratios are given in 

table 2.2 t o table 1.1 . _ The relaxation time for the formation of critical 

— *' 

nucl eus jn absence of electric fie ld is determined by temperature and 
supersaturation. Typical values of relaxation time (t;^,T;i^,Tojv) at E = 2, 
4, 6 e.s.u. are given in table (2.2) to table (2.7) for both water and ice 
phase in Homogeneous nucleation. Graphs between variation of critical 
radius and relaxation times with supersaturation ratios and E = 2, 4, 6 

'■ t 

(e.s.u.) have been plotted from fig. 2.2(a) to fig. 2.2(e iii) at T=278K for 
water phas e. Similarly for T=298K graphs are plotted from fig. 2.3(a) to 
fig. 2.3(e iii). Similarly graphs have been plotted fo r ice phas e also from 
fig. 2.4(a) to fig. 2.4(e iii) for T=278K and E = 2,4,6 e.s.u. and fi-om fig. 
2.5 (a) to fig. 2.5 (e iii) for T=298K and E=2,4,6 e.s.u. 


2.4 CONCLUSION; 

In the present study, we have calculated the relaxation time (tq) in 
absence of electric field during homogenous nucleation process. Also, the 
composition has been made with the relaxation times calculated in 
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presence of external electric field during these nucleation process using 
Murino’s expression (t^) and the expression obtained by eq. (2.31) 
forr^v- Th® reduced relaxation time (tqjv) has also been calculated both 
for water vapour condensation and ice glaciations as the function of 
supersaturation. Thus from these results it may be concluded that under 
the combined effect o f ex ternal electric field, the droplets and ice crystals 

are nucleate'd in very less time. From the tables we conclude that the 

' - ^ 

radius of critical nucleus of embryo is inversely proportional to the 
applied supersaturation ratios. Thus if supersaturation ratio is increased, 
s mall drops of water ^ idi^ crystals are formed. Similarly the relaxation 
times decreases with increase in applied electric field and is inversely 
proportional to the applied supersaturation ratios. 




THESIS OF DEVESH DWIVEDI 


Table-2.2 


Value of Critical radius (rj^ i4) and relaxation time in water vapour 
condensation at T=278 K and E= 2, 4, 6 e.s.u. as a function of 
supersaturation ratio(5|7^^), for homogenous nucleation. Suffix w 

stands for water phase 




To 

ill Sec.) 


TJV 

CuSec.) 

178009.6918 


Tm 

(fiSec.) 

314483.7889 


Tqjv 

(jx Sec . ) 
105396.6620 




48254.8399 
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Table-2.3 

Value of Critical radius (r^J^ A) and relaxation time in water vapour 
condensation at T=288 K and E= 2, 4, 6 e.s.u. as a function of 
supersaturation ratio(5,7^), for homogenous nucleation. Suffix w 

stands for water phase 









1113.6296 


227.1151 


27.3291 


15.9865 


10.0863 


9.9368 


To 

ill Sec.) 

E 

(e.s.u.) 

(At 5ec. ) 

TjVf 

(At Sec . ) 

To/V 

(At Sec . ) 


2 

164657.1213 

290894.2476 

100567.9930 

2293.6307 

4 

82328.5606 

145447.1238 

62434.6465 


6 

54885.7071 

96964.7492 

45269.4048 


2 

58432.4399 

103230.6438 

40265.4886 

576.2640 

4 

29216.2200 

51615.3219 

23838.5018 


6 

19477.4800 

34410.2146 

16931.1520 


2 

5381.60582 

9507.5036 

4470.6935 

23.9680 

4 

2690.8029 

4753.7518 

2442.0198 


6 

1793.8686 

3169.1679 

1679.7823 



396.8576 



0.3470 


0.1188 


0.0473 


0.0459 


112.3182 


74.8788 


100.5014 

50.2507 

33.5005 

50.3667 

25.1834 

16.7889 

49.2510 

24.6255 

16.4170 


198.4288 


132.2859 


177.5525 

88.7762 

59.1842 

88.9812 

44.4906 

29.6604 

87.0101 

43.5051 

29.0034 


73.1553 


95.3472 

48.9282 

32.9075 

48.2931 

24.6541 

16.5520 

47.2378 

24.1117 

16.1871 













THESIS OF DEVESH DWIVEDI 


Table-2.4 

Value of Critical radius (r^ A) and relaxation time in water vapour 
condensation at T=298 K and E= 2, 4, 6 e.s.u. as a function of 
supersaturation ratio(S^ jy), for homogenous nucleation. Suffix w 

stands for water phase 




To 

Sec . ) 


Tw 

(juSec.') 


tm 

(ji Sec . ) 


151833.7055 268239.5376 


1.005 2104.8292 2140.5773 


1055.0326 537.8101 


75916.8503 

50611.2335 

53881.7484 

26940.8742 

17960.5828 


134119.7688 

89413.1792 

95191.0888 

47595.5444 

31730.3630 


TfliV 

in Sec . ) 

95634.7698 

58676.4891 

42321.3170 

38050.9649 

22301.6644 

15773.1536 






92.6744 

46.3372 

30.8915 

46.4442 




163.7248 

81.8624 

54.5749 

82.0514 




88.2742 

45.2104 

30.3866 

44.6753 


22.7077 




40.1169 




22.2701 





2 


4962.4888 

8767.0636 

4177.5877 


4 


2481.2444 

4383.5318 

2268.1684 


6 


1654.1629 

2922.3545 

1556.6718 




-4 


, ^4"'/ 



,! S ’- 1 s ,i 


1 

23.2221 

41.0257 

22.7713 


) 

15.4814 

27.3505 

15.2797 


) 

45.4154 

80.2338 

43.6981 




11 
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Table-2.5 

Value of Critical radius and relaxation time in ice deposition 
at T=278 K and E= 2, 4, 6 e.s.u. as a function of supersaturation 
ratio(S^, i), for homogenous nucleation. Suffix i stands for ice phase 


s. 

No. 

S • 

rl(A) 

To 

ifiSe.) 

E 

(e.s.u.) 

Tw 

(jttSec.) 

Tm 

(^fiSec.) 

Tow 

(fiSec.) 





2 

350490.5271 

619199.9000 

153685.9000 

1. 

1.005 

3676.3958 

4208.8848 

4 

175245.2635 

309600.0000 

106838.5000 





6 

116830.1757 

206400.0000 

81879.5100 





2 

124379.7808 

219737.6000 

65235.7000 

2. 

1.01 

1842.7706 

1057.4627 

4 

62189.8904 

109868.8000 

42791.8700 





6 

41459.9260 

73245.8700 

31838.2000 





2 

11455.3312 

20237.7500 

8127.6400 

3. 

1.05 

375.8171 

43.9821 

4 

5727.6656 

10118.8800 

4754.3760 





6 

3818.4437 

6745.9170 

3359.8970 





2 

478.1629 

844.7544 

418.6970 

4. 

1.5 

45.2226 

0.6368 

4 

239.0814 

422.3772 

223.2292 





6 

159.3876 

281.5848 

152.1830 





2 

213.9281 

377.9397 

192.9670 

5. 

2 

26.4535 

0.2179 

4 

106.9641 

188.9698 

101.4538 





6 

71.3094 

125.9799 

68.8176 





2 

107.2110 

189.4061 

98.6953 

6. 

3 

16.6903 

0.0867 

4 

53.6055 

94.7031 

51.3885 





6 

35.7370 

63.1354 

34.7379 





2 

104.8361 

185.2105 

96.5662 

7. 

3.05 

16.4429 

0.0842 

4 

52.4181 

92.6053 

50.2657 





6 

34.9454 

61.7368 

33.9755 
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TabIe-2.6 


Value of Critical radius (xlA) and relaxation time in ice deposition at 
T=288 K and E= 2, 4, 6 e.s.u. as a function of supersaturation 
ratio(5|7 i), for homogenous nucleation. Suffix i stands for ice phase 


s • 

^V.l 

rlHA) 

To 

in Sec. ) 

1.005 

3505.9873 

3943.4947 


990.7846 



([I Sec . ) 
326406.1904 
163203.0952 
108802.0635 
115832.8893 
57916.4447 
38610.9631 


10668.1657 

5334.0828 


(ji Sec. ~ 


■Tojv 
Sec . ) 


576650.9000 148869.3000 

288325.5000 102239.3000 


192217.0000 

204638.1000 

102319.1000 
68212.7000 


18847.0900 

9423.5460 


77853.5000 

62805.1300 

40724.2500 

30130.9000 


7723.3090 

4479.9860 


3155.0540 



2 


786.7062 

393.2867 


mi 

lllllllllll^^ 


393.3531 

208.8418 



6 

148.4351 

262.2354 

142.1671 


■ 


25.2273 


15.9167 


15.6808 


0.2042 


0.0813 
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TabIe-2.7 


Value of Critical radius (vlA) and relaxation time in ice deposition at 
T=298 K and E= 2, 4, 6 e.s.u. as a function of supersaturation 
ratio(Sy f), for homogenous nucleation. Suffix i stands for ice phase 


rKA) 


To 

(jiSec.') 


Tfj 

(^5ec.) 


Tm 

(fi Sec.) 


Tow 

(fiSec.] 


1.005 3337.4799 3690.4361 


1672.8910 927.2049 


341.1716 


41.0536 


38.5644 


0.5584 


303159.2778 535581.4000 143838.7000 

151579.6389 267790.7000 97553.1000 

101053.0926 178527.1000 73803.8800 

107583.1774 190063.6000 60298.0900 

53791.5887 95031.8100 38640.7600 

35861.0591 63354.5400 28429.6500 

9908.3703 17504.7900 7317.1020 

4954.1852 8752.3940 4208.9160 

3302.7901 5834.9290 2954.0730 

413.5904 730.6763 368.3410 

206.7952 365.3382 194.8282 

137.8635 243.5588 132.4402 

185.0387 326.9016 169.1463 

92.5193 163.4508 88.3680 


38640.7600 

28429.6500 

7317.1020 

4208.9160 

2954.0730 

”3683410“ 

194.8282 

132.4402 

169.1463 

88.3680 
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Supersaturation ratio Sy.w 

fig. 2.2 (b) Variation of Calculated relaxation time Toif^Sec.) with 
supersaturation ratio Sy.w at T=278K. 
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Supersaturation ratio Sy.w 

fig. 2.2 (c-I) Variation of relaxation time Sec.) with supersaturation 
ratio Sv.w ^^t T=278K and E=2 e.s.u. 
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fig. 2.2 (d-I) Variation of relaxation time x^QiSec.) with supersaturation 
ratio Svw at T=278K and E=2 e.s.u. 
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Supersaturation ratio Sy.w 

fig. 2.2 (e-I) Variation of relaxation time Sec.) with supersaturation 

ratio Sv.w T=278K and E=2 e.s.u. 
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Supersaturation ratio Sy.w 

fig. 2.2 (c-II) Variation of relaxation time Sec . ) with supersaturation 

ratio Sv.w at T=278K and E=4 e.s.u. 
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fig. 2.2 (d-II) Variation of relaxation time Sec.) with supersaturation 

ratio Sv w at T=278K and E=4 e.s.u. 
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fig. 2.2 (e-II) Variation of relaxation ti 
ratio Sv.w at T=278K and 



Relaxation Time Sec,) 


THESIS OF DEVESH DWIVEDI 


30000 H 


T = 278 K 
E = 6 e.s.u. 



0 0.5 1 1-5 


2 2.5 3 3.5 
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fig. 2.2 (c-III)yariation of relaxation time Sec . ) with supersaturation 
■ ratio Sv w at T=278K and E=6 e.s.u. 
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fig. 2.2 (d-III) Variation of relaxation time tmCm Sec.) with supersaturation 
ratio Sv.w 3^1 T=278K and E=6 e.s.u. 
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Supersaturation ratio Sy.w 

fig. 2.2 (e-III) Variation of relaxation timeTojvCM Sec.-) with supersaturation 
ratio Sv.w at T=278K and E=6 e.s.u. 



Supersaturation ratio Sy.w 

fig. 2.3 (a) Variation of critical radius of water mokcule (^4) with 
supersaturation ratio Sy.w at T=298K. 





Relaxation Time ) 


THESIS OF DEVESH DWIVEDI 


140000 i 


T = 298 K 
E = 2 e.s.u. 


120000 I 


100000 ! 


80000 1 


60000 I 


40000 


20000 


...j 


0.5 1 1-5 


Supersaturation ratio Sy.w 

fig. 2.3 (c-I) Variation of relaxation time XfiifiSec.) with supersaturation 
ratio Sv.w at T=298K and E=2 e.s.u. 
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fig. 2.3 (d-I) Variation of relaxation time t^iy-Sec.) with supersaturation 
ratio Svw at T=298K and E=2 e.s.u. 
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fig. 2.3 (c-II) Variation of relaxation time x^iiiSec.) with supersaturation 
ratio Svw at T=298K and E=4 e.s.u. 
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fig. 2.3 (e-II) Variation of relaxation time roNQ^Sec.) with supersaturation 
ratio Sv.w at T=298K and E=4 e.s.u. 
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fig. 2.3 (d-III) Variation of relaxation time TmCiW 5ec.) with supersaturation 
ratio Sv.w at T=298K and E=6 e.s.u. , 
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fig. 2.3 (e-III) Variation of relaxation time towCm^^c.) with supersaturation 
ratio Sv.w 'at T=298K and E=6 e.s.u. 
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i . 

fig. 2.4 (a) Variation of critical radius of ice molecule ri*C-4)with 

supersaturation ratio Sv.i at T=278K. 
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fig. 2.4 (b) Variation of Calculated relaxation time To OuS'ec.) with 
supersaturation ratio Sy.i at T=278K. 
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fig. 2.4 (c-I) Variation of relaxation time Sec.) with supersaturation 

ratio Sv.i at T=278K and E=2 e.s.u. 
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fig. 2.4 (e-I) Variation of relaxation time with supersaturation 

ratio Sv.i at T=278K and E=2 e.s.u. 
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fig. 2.4 (c-II) Variation of relaxation time TNiliSec.)with supersaturation 
ratio Svi at T=278K and E=4 e.s.u. 
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fig. 2.4 (d-II) Variation of relaxation time 5ec.) with supersaturation 
ratio Sv.i at T=278K and E=4 e.s.u. 
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fig. 2.4 (e-II) Variation of relaxation time ToNO-^ec.)with supersaturation 
ratio Sv ; at T=278K and E=4 e.s.u. 
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fig. 2.4 (d-III) Variation of relaxation 5ec.)with supersaturation 

ratio Sv.i at T=278K and E=6 e.s.u. 
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fig. 2.5 (a) Variation of critical radius of ice molecule with 

supersaturation ratio Sv.i at T=298K. 
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fig. 2.5 (b) Variation of Calculated relaxation time Sec.)with 

supersaturation ratio Sv.i at T=298K. 
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fig. 2.5 (c-I) Variation of relaxation time %(/iSec.) with supersaturation 
ratio Sv i at T=298K and E=2 e.s.u. 
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fig. 2.5 (d-I) Variation of relaxation time TM(/^5'ec.)with supersaturation 
ratio Sv.i at T=298K and E=2 e.s.u. 
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fig. 2.5 (c-II) Variation of relaxation time %(M5ec.)with supersaturation 
ratio Sv.i at T=298K and E=4 e.s.u. 
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fig. 2.5 (e-II) Variation of relaxation time with supersaturation 

ratio S„.- at T=298K and E=4 e.s.u. 
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fig. 2.5 (c-III) Variation of relaxation time Tw(pt 5ec.)with supersaturation 
ratio Sv.i at T=298K and E=6 e.s.u. 
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fig. 2.5 (d-III) Variation of relaxation time T^CiU^ec.) with supersaturation 
ratio Sv : at T=298K and E=6 e.s.u. 
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EQUIVALENCE OF THE EXTERNAL ELECTRIC FIELD 
WITH SUPERSATURATION RATIO DURING WATER 
VAPOUR CONDENSATION AND GLACIATION 


8.1 INTRODUCTION: 

Water in atmosphere exists in all three forms-vapour, liquid and 
solid (ice). The phase changes take place due to change in tempeiature in 
its cycle through the atmosphere. First we will briefly discuss the basic 
stages of the cycle and examine the details later. 

1. Water is transferred from the earth surface, primarily the oceans, 
into the atmosphere in the first stage of evaporation liquid water 
undergoes what is know as a change of phase; it vapourizes and 
becomes water vapour, a gas. 

2. This water vapour condenses and forms cloud. In this second stage 
the water vapour undergoes the revers e^ charge o f phase and again 
becomes water or ice. 

3. In the third stage droj^of water fall back to earth from the clouds. 
Supersaturation ratio is defined as 

~ ^sat.w 

Where ‘e’ is vapour pressure of water vapour over water and 
e^at.w. the maximum possible saturated values. 

Air for which relative humidity = 1(100%) is saturated, if 
Ov > 1 = (> 100%), it is said to be supersaturated, corresponding to 

Supersaturation 5^ > 0. 
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According to Goff (1965). = 0.01«C) = 6.11 12inb, and 

^sat.w (^0 ~ 0 C, p— latm) = 6.1067 mb. The former result is in excellent 
agreement with the most recent experimental values obtained as U.S. 
National Bureau of standards (Guider, 1977). The temperature variation 
of saturation pressure is shown in fig. (3.1). For T< O^C, e sat. w>esat.i- 
The maximum difference of e^at.uf - esat.i, can be calculated by the 
figure. The variation of Bsat.w - ^sati is illustrated in fig. (3.2). The 
interesting point to emphasis here is that air saturated with respect to 
water, with the consequence that supercooled water drops and ice crystals 
cannot coexist in equilibrium. 

It is also worth emphasizing that at sufficiently low temperatures, 
air may be ice supersaturated but water sub saturated. This is illustrated in 
fig. (3.3) similarly fig. (3.4) shows the ice supersaturation, which are 
required for water saturation to occur also. 

We shall study the equilibrium of (i) a pure water or aqueous 
solution drop surrounded by water vapour or humid air and (ii) an ice 
crystal in humid air. We shall see later that the relationships provided by 
these case studies are needed in order to formulate the conditions for 
which cloud drops and ice crystals are nucleated in the atmosphere. 

8.2 EQUILIBRIUM BETWREN A PURE WATER DROP AND PURE 
WAim Vapour OR HUMID AIR: 

The Kelvin equation has been derived by W. Thomson according 

. ,-i Saw f2Vw.o(^wla\ 1 -y 


PRESSURE (mm Hg) 
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fig. 3.1 : P-T phase diagram for bulk water substance. 
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fig. 3.2 : Variation of Aesat= esat.w” esat.i with temperature. 



Sw = e/e 
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fig. 3.3 : Water Supersaturation as a function of saturation ratio with 
respect to ice at temperature below O^C. 


Sw = e/e 
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where Mw is the molecular weight of water; , the surface free energy 
of water air interface; a the radius of water droplet; p^, the density of 
water; T, the temperature; R, the universal gas constant and = — is 
the molar volume of pure water in bulk. 

According to this equation at any given temperature saturated 
vapour pressure over the surface of water drop is larger than that over flat 
surface, and increases as the decreasing radius. 

Since Kelvin equation assumes equilibrium between the drop and 
its environment. We have e^.w = partial pressure of vapour in the 
environment. 


Hence we can say that equilibrium requires an environment 
supersaturation of 



1 > 0 


Also the Kelvin equation may be expressed in terms of the 
saturation ratio 5^.^ = ^a.wl^sat.w i^^ 


{tiS = 

l^iLUpyyf 


2Vw,o^w/a ^^w^wfa 

RTa “ RTpw.a 


(3.2) 


A numerical evaluation of eq. (3.2) is plotted in fig. (3.5) for 20®C 
and -20® C. Dufour (1963) in deriving the Kelvin equation assumed that 

compressibility of water is negligible. 
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fig. 3.5 ; Variation of the equilibrium vapour pressure over a water sphere 
and over an ice sphere with sphere size, for different 
temperatures. For (1) and (2), ordinates represents eo/esat.w for 
(3) and (4), ordinates represent eo/esau. 

— 
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8.8 EQUILIBRIUM OF ICE crystal IN HUMID AIR: 


There are three equilibrium situations involving the ice phase. 

1 . An ice particle in Humid air 

2. An ice particle and separated supercooled solution drop in 
Humid air. 

3. An ice particle in supercooled solution drop in humid air. 
These three cases are illustrated in fig. (3.6) 

In present case we analyze only case (1) and assuming that ice 
particle is spherical in strict analogy to the derivation of eq. (3.1) and 
define the saturation vapour pressure over ice, varies with radius 
according to 


‘'V.W 


^sat 


^ = exp = exp 

ati ^ ^ ^ RTpja / 


(3.3) 


This result is calculated for T— 263*^K. Shown in fig. (3.7) (a,b). 
The behaviour is seen to be very similar to that for a pure water drop, 
except that at small radii , Sjr.j » 

8.4 ELECTRICAL STATE OF ATMOSPHERE: 

We shall describe some observed and predicated modification in 
the behaviour of isolated and interacting cloud particle in consequence of 
the ambient electrical field. 
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fig. 3.6 ; The three basic equilibrium types involving ice crystals. 
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fig. 3.7 : Variation of Critical radius of water and ice molecule and 
Gibb’s energy function for water and, ice molecule at T = 263 
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Work carried out prior to 1963 suggested qualitatively that electric 
fields and the charges enhance ice nucleation (Pruppacher, 1963). 
Gabarashivili and Gliki (1967) and Gabarshivili and Kartsivadze 
(1968,1969) found that supercooled drops containing particles of quartz 
or naphthalene were nucleated to ice at significantly wanner temperatures 
when the particles carried a net negative charge than when they were 
neutral or carried a net positive charge. Abbas and Latham (1969), 
Morgan and Langer (1973) observed that charge nuclei produced during 
corona discharges or by sparks promoted ice nucleation of supercooled 
drops . The negatively charged ions are more effective than the positively 
charged ions (Singh, N., et al., 1985). Pruppacher (1973) found that the 
freezing temperature of supercooled water drops of 100 to 350 ixm radius, 
freely suspended in the air stream of a wind tunnel, was raised 
considerably when contacted by predominantly negatively charged 
amorphous sulfur particles which, when unchanged, are known to be poor 
ice forming nuclei. 

The effect of external electrical fields on ice nucleation has been 
studied under essentially two different experimental conditions. 
Pruppacher (1963) observed in laboratory experiments that millmeter 
sized drops, forced in an external electric field to deform and thus to 
rapidly spread over a solid surface, froze at temperatures up to 10°C 
warmer than when the drops were unaffected by the field. These results 
are consistent with the observations (Doolittle, 1975), who found that an 
electric field had no effect on the freezing of supercooled drops if they 
remained motionless with respect to the supporting surface. Other studies 
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have dealt with drops in free fall. For example (Dawson et al, 1973), 
observed millimeter sized drops suspended in the air stream of wind 
tunnel at temperatures -8 to -15^C, and detected no electro freezing effect 
for external field upto 4 KVcm Coalescence between drops did not alter 
this outcome. In contrast to these studies, however (Smith et al., 1971; 
Abbas and Latham, 1969) found that millimeter sized drops falling in air 
of -5 to -12°C through intense fields froze if disruptions caused small 
filaments to be drawn out from or between drops. 

The theory of droplet growth in presence of external electric fileld 
has been applied to water vapour condensation (Singh, N.,et al., 1986). 
The radius of critical nucleus is found to depend on electric field and 
relaxation time. Further in presence of electric field, Gibb’s free energy 
of formation of water molecules about a critical nucleus is found to be 
much less than that in absence of the field, but equilibrium concentration 
of nuclei is very much large. From the observations we come to know 
that a small value of electric field is equivalent to very high 
supersaturation ratio to get nucleus of given size under similar conditions 
of temperature. 

It is well know fact that there exists a vertical electric field in free 
and undisturbed atmosphere. Also, intense electric fields are generated in 
thunder clouds (Sapkata et al, 1989). The intense electric fields in 
thunder clouds are not widespread and the charge and intense electric 
fields seem to be concentrated in relatively very small volume. Under 
cloud condition the rate of growth of electric field varies considerably 

and has been calculated by Winn et al (1974). 
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Moore et al. (1962) and Levin and Ziv (1974) reported a raingush 
following lightning stroke. Murino (1974) has shown that action of a 
constant and uniform electric field accelerated the condensation of water 
vapour by a factor depending upon the intensity of electric field. Changes 
in scaveging of particles due to weak electrification has been studied 
(Brian A. Tinsleya et al., 2006). Mahendra Fernando (2007) studied 
propagation effects on the electric field time derivatives generated by 
return strokes in lighting flashes. 

In this chapter, we have considered the resultant effect due to an 
external electric field and the field due to the central dipole. Following 
Kinetic theory (Collins, 1955), inferred that the relaxation time is 
independent of the free energy of formation of the nucleus but it varies as 
the square of the radius of the critical nucleus. 

3.5 THEORETICAL CONSIDERATION: 

3.5.1 Water Vapour Condensation: 

In atmosphere under suitable conditions the phase change of water 
vapour may take place from, 

1 . Water vapour to liquid water (Condensation) 

2. Water vapour to ice (deposition) 

3. Water vapour to ice via liquid phase. 

It has been reported that the last mode of phase change is less 
effective than in (2) and rigorous mathematical consideration is required. 
Therefore, we consider the first two cases only. In the classical model of 
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drop formation the Gibb’s function has been taken into consideration. 
However (Abraham, F.F., 1968; Dufour et al., 1963), have shown that the 
Helmholtz free energy is the proper thermodynamic potential and Gibb’s 
function is only its approximation. For practical purposes, the resulting 
difference turns out to be negligible. This energy of the formation of 
nucleus is given by 

+ (3.4) 

With AGy, = the volume free energy of condensate per 

Mw 

unit volume per mole where r^ is the radius of water nucleus; 
surface free energy per unit area (surface) of water vapour interface; p^,, 
the density of water; T, the temperature; , the molecular weight of 
water; R, the universal gas constant and 5^,^, the supersaturation ratio of 
water vapour over water surface. 

The Gibb’s free energy at a given temperature and supersaturation 
increases first with increasing radius of the condensed drop attains a 
maximum and then decreases monotonically. The critical radius 
corresponding to the maximum, Gibb’s free energy is significant in the 
sense that above this value the embryo will grow to form a drop because 
it comes to the lower state in doing so. This critical radius obtained by 
setting. 

£(^ = 0 (3-5) 

dryj 

— C£L>- 



I 
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From equation (3.4) and (3.5) we have the radius of critical 
nucleus as 


“^^w/v ^^w^w/v 




S’p Vi 


(3.6) 


Substituting equation (3.6) back into (3.5) the critical Gibb’s 

energy of formation of a critical nucleus is 


AG, 




w 


^w/i? 


AG* - 

3iPwRTinS^,y,)2 


ag* - 


'w 


3 (AGv )2 


(3.7) 


The number of water molecules in a critical nucleus is given by 


K = 


AnryfpwN 

3Mw 


(3.8) 


N. 


Where N is Avogadro’s number. The critical supersaturation ratio 
can be obtained from eq. (3.5) for a neutral drop, it is given by 


£n 5 ^.vi 7 = 


PwRT 


(f) 


(3.9) 


Which is Kelvin classical formula, also sometimes called the 
Thomson, or Gibb’s Thomson formula. The expression for equilibrium 
concentration for the critical nuclei is 

C(n;) = C(l)„e;tp[-AG;/kT] 




(3.10) 
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Where C(7i^,) is the equilibrium concentration of critical nuclei 
having n^i/Water molecules in it; C (l)o, concentration of monomers 
(water vapour molecules) and k Boltzmann constant. 

3.5.2 Nucleation Rate (J): 

When homogenous nucleation of water from the vapour occurs, 
what is observed is the formation of a cloud of small drops. Thus the 
experimental quantity of interest is the rate at which drops appear in the 
system as a function of the prevailing saturation ratio let us denote 
this rate by J, measured as the number of drops appearing per unit volume 
and per unit time. We shall make the traditional assumption that J 
corresponds completely to the rate of germ formation; i.e., it is the 
nucleation rate. 

Further the rate of nucleation of critical nuclei is given by 
Pruppacher and Klett (1978), 


I 


W 


PwL TT J V RT 7 



(3.11) 


Where is condensation coefficient 

An inspection of the above equation shows that C{nl,) and are 
extremely sensitive to supersaturation. This behaviour enables us to 
define a critical supersaturation ratio at which drops suddenly appear in 
vapour. By convention this is taken to correspond to 7^,, = 1 germ cm 
sec'' . Pruppacher et al. (1978) calculated the pre-exponential factor for 
oc^rr 1 is 9 xlO^' cm'^ sec'' and 3.7 x 10^^ cm'^ see' at -30°C and 30''G 
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respectively. Although pre-factor changes with temperature and 
supersaturation ratio, but it has been assumed constant for present 
calculations, because the nucleation rate is mainly dependent on the 
exponential term. 


Hence eq. (3.1 1) become 


Jw = 9xl0"'exp[:^] 


(3.12) 


The computed results for, rv^,n;, , AG^ and are given in table (3.1). 


3.5.3 Water Vapour Ice Nucleation: 


The equation analogous to eq.(^, eq.(3_T^, eq.(3.11) for 
homogeneous water condensation hold good for spherical ic e germ 
nucleation from water vapour with 


(Zc ^di^v.w ^v.i'Pw Pi’^w/v ^tjv OTld 6sat.w ^sat.i' 


The suffix ,‘i’ stands for ice and being deposition coefficient. 


The computed results for rl,nl, AG^ and are given in table (3.2) 
with Pi = 0.917 gm cm~^,aify = 100 dynes cm~^ and = 1. 

From equation (3.9) shows that varies approximately as 

(r^*)“^ for neutral drop. 
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3.6 EFFECT OF EXTERNAL ELECTRICAL FIELD ON NUCLEATION 

DURING Water Vapour condensation and ice glaciations : 

As water is strongly polarisable medium, an electric dipole 
moment on both the water droplet and surrounding water vapour. The 
moment induced on the drop is given by (Kittle et al., 1966). 

r 

M = ErS (3 13 ) 

Where r^ is the radius of droplet and E the inducing electric field. 
The moment induced on water vapour molecules is given by 

M = aE (3.14) 


WTiere a is polarizability. 

Water is strongly polarizable medium. An external electric field 
induces an electric dipole moment on both the water droplet and 
surrounding water vapour. The expression for rate of increase in the 
radius of water droplet has been derived as 

^ = (^) (3.15) 

Where A is the mean free path of water vapour molecule; E, the 
external electric field; m^, the mass of water vapor; a, is 
polarizability; Pi;, the density of water vapour; the radius of water / 

n ucleus in presence of electric field. Integrating the above equation 
within limits r.^ = 0 to = r^* (critical radius of w ater nucleus) and 
t=0 to T(relaxation time), we get 
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r’* = 

'w 


3pvT[9aAE2/i^/2 


2Pw 


2/3 


(3.16) ^ 


The expression for the Gibb’s free energy for the formation of 
critical water nucleus. 


AC 


4irn«?<riv/i< 



(3.17) 


And the equilibriuni concentration of the critical nuclei is given by 


C(n;;) = C(l)oexp(-AG^7/cT) (3.18) 

Where k is Boltzmann constant and C (l)o the number of water 
molecules (monomers) at temperature T. 

3.7 EQUIVALENCE BETWEEN ELECTRIC FIELD AND 

SUFERSATURATION RATIO: 

From above we see that in absence of electric field the main 
parameters to the nucleation of water droplets and ice glaciation are 
temperature and supersaturation. By thermal fluctuation at given 
supersaturation the water droplets are nucleated. Hence, supersaturation is 
the main parameter in the nucleation phenomenon. 

Similarly in presence of electric field the main parameter to the 
nucleation phenomenon is the magnitude of electric field. 

Thus, it appears that supersaturation and external electric field are 
analogous to each other. In the absence of electric field the radius of 
critical nucleus has been calculated (Gupta, N. et al., 2007). Recently the 
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equivalence between electric field and supersaturation ratio has been 
studied briefly (Dwivedi, D. et. al., 2008). I ^present s tudy, we discussed 
the equivalence between the two parameters. If the same size of the water ^ 

nucleus is obtained m the two cases ; in absence of electric field and in ' 
presence of external field, we may equate eq.(3.6) with eq.(3.16), Thus; 

2A^w(^w/p fsppTE . 

p^RTinS^,^ ~ [ 2p^ J (3.19) ' 

r 2Mwcrw/v 1 ^/^ 2 Ph, ^ _ r 

IpwRT tn Sp.v^j 3 p„T V aA y ^3 ^ 

lRT tn Sj; 9 ppT \p,^,aAy 


or 

E = — iL_ X / AWlyy ^ / ^W^^W/V ^^^^1 

[9pvT Kp^aV \RTinSv.w) \ 

eq. (3.20) may be written as 


4 

X 


f \3/2' 

9ppT 

Kp^ax) 

Ur in Sv.w/ 


(3.20) 


(3.21) 


Where E = Egq, is equivalent electric field corresponding to 
supersaturation ratio 5i, VI,. 

From eq. (3.20) we can write 

g2/3 ^ X X (3.22) 

\9pvrJ Kp^aX) \RT in Sy,,^J ^ ' 

From eq. (3.22) we can write 

(3.23) 
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From eq.(3.23) we can calculate (equivalent 

supersaturation ratio at a given electric field) and the values of Eeq has 
been calculated from eq. (3.21). 

3.6 RESULT AND DISCUSSION: 

In table (3.1) and table (3.2) the cri+ical properties of uncharged 
water and ice molecules at T=273K as a function of supersaturation ratios 
have been calculated. The graphs of these tables have been plotted from 
fig. 3.8 (a-c) to fig. 3.9 (a-c). From the graphs it is clear that at high 
supersaturation ratio the critidalradius of the molecule decreases and the 
no. of molecules also decreases. The Gibb’s free energy decreases 
accordingly increase in supersaturation ratios. From eq. (3.21) the values 

of equivalent electric field Eg. in table (3.3 a) to table (3.3 i) and critical 

- — — 

radius have been calculate^for different values of supersaturation ratio 
at given temperature. Similarly the values of {S.^,w)eq (3-4) are 

calculated by eq. (3.22) and are calculated by eq. (3.16). 

In Present conditions the constant values used are: X = 10“^ cm; 
Ptj = at 10°C; a = 5 x = Ijmvi, = 3 x 

10-23^. ;i? = 8.31 X Kf ergK-^mole-'^-, (3 = 1.1267 x = 

1.38 X lO-i^er^ = 6.023 X 10-23. 

From the tables, it is clear that at a given relaxation times and 
temperature, the equivalent electric field decreases with increase in 
supersaturation ratio. Also, the equivalent electric field decreases with 
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increasing temperature. Corresponding to a given supersaturation ratio at 
a constant relaxation time. 

Table (3.3a) to table (3.3i) represents the equivalence of 
supersaturation ratio with the electric field at T=278K, 288K, 298K and 
relaxation times t = 5, 6, 7 (xlO ^sec.). The equivalence depends on 
temperature and relaxation time as well. The graphs are plotted from fig. 
(3.10a) to fig. (3.10i) in supersaturation ratio and equivalent electric field. 

Similarly, Table (3.4) compares the equivalence of electric field 
with supersaturation ratio as function of temperature and relaxation time. 
The graphs have been plotted from fig. (3 . 1 1 a) to fig. (3 . 1 1 c). 

From Table (3.4) it is clear that at a given relaxation time and 
temperature, the equivalent supersaturation ratio decreases with increase 
in electric field. 

3.9 CONCLUSION: 

The nucleation rate of water condensation increases very much 
due to applied external electric field. Hence the electric field is very much 
effective in nucleation process. A small external electric field is 
equivalent to be large value of supersaturation ratio, which otherwise 
never exist in the clouds. In electric field induced nucleation, critical 
nuclei of smaller size are formed having smaller number of water 
molecules and require less Gibb’s free energy for there formation and 
equilibrium concentration of critical nuclei is enhanced as compared with 
the homogeneous nucleation. 
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Also, from the equivalence of supersaturation ratio with external 
electric field, it appears that electric field is more efficient than 
supersaturation. Although, there are a number of mechanisms giving 

possible explanation for rainguish after lightning discharge, creating 

ionization of medium and a large number of ions are formed. Most of 
these ions get attached with water embryos of critical nuclei. Further 
condensation of water vapour on critical nuclei have less radii, requiring 
less energy of formation (Sapkata, et.al., 1989). The number of critical 
nuclei is enormously increased using eq. (3.8). Hence, for a given amount 
of liquid water a large number of droplets of smaller size are formed 
immediately following lightening. Hence, Lightening seems to be cause 
of raingush and not the effect. 
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Table-3.1 

Critical Properties of Uncharged water molecules at T=273K, as a 
function of supersaturation ratio ^ 


s . 

No 


rUA) 


ergs. 

AG*^/kT 


1. 

1.2 

64.40691969 

37410.19383 

128.5183564 

3411.32761 

-3360.8376 

2. 

1.4 

34.89966953 

5951.909441 

37.73485919 

1001.615416 

-951.12542 

3. 

1.6 

24.98442381 

2183.737266 

19.33921229 

513.3304744 

-462.84047 

4. 

1.8 

19.97794534 

1116.46428 

12.3652171 

328.2161994 

-277.7262 

5. 

2 

16.94123585 

680.8126814 

8.891812202 

236.0198599 

-185.52986 

6 

2.2 

14.89334802 

462.5613956 

6.m0257%4 

182.4076494 

-131.91765 

7. 

2.4 

13.41312301 

337.8952551 

5.57390958 

147.9510957 

-97.461096 

8. 

2.6 

12.28951252 

259.8940517 

4.67917639 

124.201741 

-73.711741 

9. 

2.8 

11.40496155 

207.7178637 

4.029839556 

106.9660656 

-56.476066 

10. 

3 

10.68872976 

170.9899803 

3.539584614 

93.95298121 

-43.462981 

11. 

3.2 

10.0956555 

144.0774249 

3.157687508 

83.81609353 

-33.326094 

12. 

3.4 

9.59552673 

m.mim 

2.852579413 

75.71745535 

-25.227455 

13. 

3.6 

9.167350763 

107.8757721 

2.603681168 

69.11082358 

-18.620824 

14. 

3.8 

8.796075268 

95.29260112 

2.397054886 

63.62623788 

-13.136238 

15. 

4 

8.470617926 

85.10158518 

2.22295305 

59.00496497 

-8.514965 

16. 

4.2 

8.18263291 

76.71346079 

2.074370086 

55.06105234 

-4.5710523 

17. 

4.4 

7.925711309 

69.71193195 

1.94615131 

51.65767665 

-1.1676767 
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Table-3.2 


Critical Properties of Uncharged ice molecules at T=273K, 
function of supersaturation ratio ^ 


1.2 94.91426167 109788.317 377.1649546 lOOi 1.27979 

1.4 51.4304423 17467.16747 110.7411179 2939.457394 

1.6 36.8186858 6408.650017 56.75510746 1506.479468 

1.8 29.44081073 3276.506262 36.28840795 963.2215307 

~2 24.96571642 1997.992282 26.09494892 692.6514021 

2.2 21.94781458 1357.486609 20.16744818 535.3147577 

2A 19.76645791 991.6268161 16.357845 434.1945375 

I 18.11063179 762.7153893 13.73205665 364.4969118 

16.80709946 609.5930639 11.82643706 313.915089 


AGl/kT 


15.75161331 501.8071345 10.32767714 275.7253741 


14.87762017 422.8264113 I 9.266917965 1245.9764815 


14.14059762 363.0487945 8.371512172 222.2092736 


13.50960942 316.5848197 7.641066359 202.8206816 


£nJi 

-9960.7898 

-2888.9674 

-1455.9895 

-912.73153 

-642.1614 

-484.82476 

-383.70454 

-314.00691 

-263.42509 


-225.23537 


-195.48648 


-171.71927 


-152.33068 


12.9624735 

279.6567789 

7.034676779 

186.7249769 

-136.23498 

12.48285821 

249.7490352 

6.523737231 

173.1628505 

-122.67285 

12.0584646 

225.1322673 

6.087688338 

161.5885846 

-111.09858 

11.67984808 

204.5847643 

5.71 1402567 

151.6006415 

-101.11064 
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Table-3.3 (a) 

Values of Ecq Corresponding to supersaturation ratio 5^ ^ for given 
function o f t = 5 X 10~^ sec. and T=278K 


S. No. 

^V.W 

r;;,(x 10"^ cm.) 

Eeq(x ^O^v/cm.) 

1. 

1.2 

64.02416305 

258.4834456 

2. 

1.4 

34.69226821 

103.1017204 

3. 

1.6 

24.83594669 

62.45079161 

4. 

1.8 

19.85922066 

44.65397978 

5. 

2 

16.84055769 

34.86996997 

6. 

2.2 

14.80484002 

28.7423615 

7. 

2.4 

13.33341167 

24.56568064 

8. 

2.6 

12.21647856 

21.54447546 

9 

2.8 

11.33718429 

19.26081583 

10. 

3 

10.62520891 

17.47523747 

11. 

3.2 

10.03565916 

16.04115827 

12. 

3.4 

9.538502553 

14.86405042 

13. 

3.6 

9.112871145 

13.88032751 

14. 

3.8 

8.743802061 

13.04569876 

15. 

4 

8.420278843 

12.32839611 

16. 

4.2 

8.134005261 

11.70505727 

17. 

4.4 

7.878610491 

11.15812736 

18. 

4.6 

7.649118222 

10.67416543 

19. 

4.8 

7.441582701 

10.24270909 

20. 

5 

7.25283342 

9.855494585 

21. 

5.2 

7.080292365 

9.505909083 

22. 

5.4 

6.921840798 

9.188598576 

23. 

5.6 

6.775720573 

8.899182246 

24. 

5.8 

6.640459938 

8.634041082 

25. 

6 

6.514817016 

8.390159218 
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Table-3.3 (b) 

Values of Ecq Corresponding to supersaturation ratio 5^^ for given 
r■^^, as the function of t = 6 x 10“^ sec. and T=278K 


I 


No. 

€ 

ri^(x 10 ®cm.) 

Eeq(x 10^ v/ cm.) 

1. 

1.2 

64.02416305 

215.4028713 

2. 

1.4 

34.69226821 

85.91810037 

3. 

1.6 

24.83594669 

52.04232634 

4. 

1.8 

19.85922066 

37.21164982 

5. 

2 

16.84055769 

29.05830831 

6. 

2.2 

14.80484002 

23.95196792 

7. 

2.4 

13.33341167 

20.47140053 

8. 

2.6 

12.21647856 

17.95372955 

9. 

2.8 

11.33718429 

16.05067986 

10. 

3 

10.62520891 

14.56269789 

11. 

3.2 

10.03565916 

13.36763189 

12. 

3.4 

9.538502553 

12.38670868 

13. 

3.6 

9.112871145 

11.56693959 

14. 

3.8 

8.743802061 

10.87141563 

15. 

4 

8.420278843 

10.27366343 

16. 

4.2 

8.134005261 

9.754214389 

17. 

4.4 

7.878610491 

9.29843947 

18. 

4.6 

7.649118222 

8.895137858 

19. 

4.8 

7.441582701 

8.535590908 

20. 

5 

7.25283342 

8.212912155 

21. 

5.2 

7.080292365 

7.921590902 

22. 

5.4 

6.921840798 

7.65716548 

23. 

5.6 

6.775720573 

7.415985205 

24. 

5.8 

6.640459938 

7.195034235 

25. 

6 

6.514817016 

6.991799348 
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Table-3.3 (c) 

Values of Eeq Corresponding to supersaturation ratio 5^, for given 
as the function of x = 7 x 10“^ sec. and T=278K 


s . 

No. 


1. 

1.2 

2. 

1.4 

3. 

1.6 

4. 

1.8 

5 

2 

6. 

2.2 

7. 

2.4 

8. 

2.6 

9. 

; 2.8 



^(x IQ-^ cm.) 

64.0241 6305~ 
34.69226821 
24.83594669 
19.85922066 
16.84055769 
14.80484002 
13.33341167 
12.21647856 
11.33718429 
10.62520891 
10.03565916 
9.538502553 
9.112871145 
8.743802061 
8.420278843 
8.134005261 
7.878610491 
7.649118222 
7.441582701 
7.25283342 
7.080292365 
6.921840798 
6.775720573 
6.640459938 
6.514817016 


Eeq(x 10^ v/cm.) 

184.6310325 

73.64408603 

44.60770829 

31.89569984 

24.90712141 

20.53025821 

17.54691474 

15.38891104 

13.7577256 

12.48231248 

11.45797019 

10.61717887 

9.914519649 

9.318356255 

8.805997223 

8.36075519 

7.970090974 

7.624403878 

7.316220779 

7.03963899 

6.789935059 

6.563284697 

6.356558747 

6.167172202 

5.99297087 
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Table-3.3 (d) 

Values of Ecq Corresponding to supersaturation ratio for given 
as the function of t = 5 X 10"^ sec. and T=288K 


S. No. 


r;;,(x 10~® cm.) 

Eeq(x ^O^v/cm.) 

1. 

1.2 

60.78316285 

251.8965789 

2. 

1.4 

32.936093 

100.474406 ' 

3. 

1.6 

23.57871341 

60.85937426 

4. 

1.8 

18.8539168 

43.51607398 

5. 

2 

15.98806312 

33.98138756 

6. 

2.2 

14.05539657 ' 

28.00992735 

7. 

2.4 

12.65845416 

23.93967976 

8. 

2.6 

11.59806189 

20.99546317 

9. 

2.8 

10.76327883 

18.76999745 

10. 

3 

10.08734473 

17.02992052 

11. 

3.2 

9.527638879 

15.63238559 

12. 

3.4 

9.055649093 

14.48527367 

13 

3.6 

8.651563792 

13.52661872 

14. 

3.8 

8.301177545 

12.71325859 

15. 

4 

7.994031562 

12.01423479 

16. 

4.2 

7.722249583 

11.40678032 

17. 

4.4 

7.479783283 

10.87378769 

18. 

4.6 

7.261908261 

10.40215844 

19. 

4.8 

7.064878503 

9.981696788 

20. 

5 

6.885683997 

9.60434957 

21. 

5.2 

6.721877231 

9.263672464 

22. 

5.4 

6.571446723 

8.954447899 

23. 

5.6 

6.432723325 

8.672406691 

24. 

5.8 _ 

6.304309789 

8.41402205 

25. 

6 

6.185027102 

8.176354964 
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Table-3.3 (e) 

Values of Eeq Corresponding to supersaturation ratio for given 
r;^ as the function of t = 6 x 10"^ sec. and T=288K 


S.No. 


rj^(x 10 ®cm.) 

£'eq(x 10^ r/cm.) 

1. 

1.2 

60.78316285 

209.9138157 

2. 

1.4 

32.936093 

83.72867168 

3. 

1.6 

23.57871341 

50.71614522 

4. 

1.8 

18.8539168 

36.26339498 

5. 

2 

15.98806312 

28.31782296 

6. 

2.2 

14.05539657 

23.34160612 

7. 

2.4 

12.65845416 

19.94973314 

8. 

2.6 

11.59806189 

17.49621931 

9. 

2.8 

10.76327883 

15.64166454 

10. 

3 

10.08734473 

14.19160043 

11. 

3.2 

9.527638879 

13.02698799 

12. 

3.4 

9.055649093 

12.07106139 

13. 

3.6 

8.651563792 

11.27218227 

14. 

3.8 

8.301177545 

10.59438216 

15. 

4 

7.994031562 

10.01186232 

16. 

4.2 

7.722249583 

9.50565027 

17. 

4.4 

7.479783283 

9.061489745 

18. 

4.6 

7.261908261 

8.668465364 

19. 

4.8 

7.064878503 

8.318080656 

20. 

5 

6.885683997 

8.003624642 

21. 

5.2 

6.721877231 

7.719727053 

22. 

5.4 

6.571446723 

7.462039916 

23. 

5.6 

6.432723325 

7.227005576 

24. 

5.8 

6.304309789 

7.011685042 

25. 

6 

6.185027102 

6.813629137 


dD' 
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Table-3.3 (f) 

Values of Eeq Corresponding to supersaturation ratio for given 
as the function of t = 7 x 10“^ sec. and T=288K 


S.Nu; 

^V.W 

r;,(x 10"^ cm.) 

Eeq(x 10^ 17/cm.) 

1. 

1.2 

60.78316285 

179.9261278 

2. 

1.4 

32.936093 

71.76743286 

3. 

1.6 

23.57871341 

43.47098162 

4. 

1.8 

18.8539168 

31.08290999 

5. 

2 

15.98806312 

24.27241968 

6. 

2.2 

14.05539657 

20.00709096 

7. 

2.4 

12.65845416 

17.09977126 

8. 

2.6 

11.59806189 

14.99675941 

9. 

2.8 

10.76327883 

13.40714104 

10. 

3 

10.08734473 

12.16422894 

11. 

3.2 

9.527638879 

11.16598971 

12. 

3.4 

9.055649093 

10.34662405 

13. 

3.6 

8.651563792 

9.661870513 

14. 

3.8 

8.301177545 

9.080898996 

15. 

4 

7.994031562 

8.581596277 

16. 

4.2 

7.722249583 

8.147700232 

17. 

4.4 

7.479783283 

7.76699121 

18. 

4.6 

7.261908261 

7.430113169 

19. 

4.8 

7.064878503 

7.12978342 

20. 

5 

6.885683997 

6.860249693 

21. 

5.2 

6.721877231 

6.616908903 

22. 

5.4 

6.571446723 

6.396034214 

23. 

5.6 

6.432723325 

6.194576208 

24. 

5.8 

6.304309789 

6.01001575 

25. 

6 

6.185027102 

5.840253546 
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Table-3.3 (g) 

Values of Ecq Corresponding to supersaturation ratio 5^ ^ for given 
as the function of t = 5 x 10“^ sec. and T=298K 


S. No. 

^V.W 

10”® cm.) 

Eeq(x 10^ v/ cm.) 

1. 

1.2 

57.58170636 

243.9717028 

2. 

1.4 

31.20134503 

97.31339753 

3. 

1.6 

22.33681975 

58.94468767 

4. 

1.8 

17.86087874 

42.14702206 

5. 

2 

15.14596992 

32.91230481 

6. 

2.2 

13.31509714 

27.12871171 

7. 

2.4 

11.99173186 

23.1865175 

8. 

2.6 

10.98719058 

20.33492842 

9. 

2.8 

10.19637564 

18.1794777 

10. 

3 

9.556043069 

16.49414504 

11. 

3.2 

9.025816995 

15.14057772 

12. 

3.4 

8.578686968 

14.0295549 

13. 

3.6 

8.195884888 

13.10106003 

14. 

3.8 . . 

7.863953527 

12.31328889 

15. 

4 

7.57298496 

11.63625696 

16. 

4.2 

7.315517771 

11.04791352 

17. 

4.4 

7.0858222 

10.53168928 

18. 

4.6 

6.879422681 

10.07489787 

19. 

4.8 

6.692770504 

9.667664295 

20. 

5 

6.523014194 

9.302188735 

21. 

5.2 

6.36783515 

8.972229615 

22. 

5.4' 

6.225327835 

8.672733512 

23. 

5.6 

6.093911015 

8.399565555 

24. 

5.8 

5.972261036 

8.14930991 

25. 

6 

5.859260982 

7.919120029 


( 195 > 
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Table-3.3 (h) 

Values of Eeq Corresponding to supersaturation ratio for given 
rl, as the function of t = 6 x 10“® sec. and T=298K 


S. No. 

^V.W 

rj^(x 10 ®cm.) 

Eeq{><10^v/cm.) 

1 

1.2 

57.58170636 

203.3097523 

2 

1.4 

31.20134503 

81.09449795 

3 

1.6 

22.33681975 

49.12057306 

4 

1.8 

17.86087874 

35.12251838 

5 

2 

15.14596992 

27.42692067 

6 

2.2 

13.31509714 

22.60725976 

7 

2.4 

11.99173186 

19.32209792 

8 

2.6 

10.98719058 

16.94577368 

9 

2.8 

10.19637564 

15.14956475 

10 

3 

9.556043069 

13.74512087 

11 

3.2 

9.025816995 

12.6171481 

12 

3.4 

8.578686968 

11.69129575 

13 

3.6 

8.195884888 

10.91755003 

14 

3.8 

7.863953527 

10.26107408 

15 

4 

7.57298496 

9.696880798 

16 

4.2 

7.315517771 

9.206594598 

17 

4.4 

7.0858222 

8.776407733 

18 

4.6 

6.879422681 

8.395748225 

19 

4.8 

6.692770504 

8.056386913 

20 

5 

6.523014194 

7.751823946 

21 

5.2 

6.36783515 

7.476858012 

22 

5.4 

6.225327835 

7.227277927 

23 

5.6 

6.093911015 

6.999637962 

24 

5.8 

5.972261036 

6.791091592 

25 

6 

5.859260982 

6.599266691 
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Table-3.3 (i) 

Values of Eeq Corresponding to supersaturation ratio for given 
as the function of t = 7 x 10“® sec. and T=298K 


S. No. 

^V.W 

r;;,(xlO~^cm.) , 

feq(x 10^ V /cm.) 

1. 

1.2 

57.58170636 

174.265502 

2. 

1.4 

31.20134503 

69.50956967 

3. 

1.6 

22.33681975 

42.10334834 

4. 

1.8 

17.86087874 

30.10501576 

5. 

2 

15.14596992 

23.50878915 

6. 

2.2 

13.31509714 

19.37765122 

7. 

2.4 

11.99173186 

16.56179822 

8. 

2.6 

10.98719058 

14.52494887 

9. 

2.8 

10.19637564 

12.98534122 

10. 

3 

9.556043069 

11.78153217 

11. 

3.2 

9.025816995 

10.81469837 

12. 

3.4 

8.578686968 

10.02111064 

13. 

3.6 

8.195884888 

9.357900022 

14. 

3.8 

7.863953527 

8.795206352 

15. 

4 

7.57298496 

8.311612112 

16. 

4.2 

7.315517771 

7.891366798 

17. 

4.4 

7.0858222 

7.522635199 

18. 

4.6 

6.879422681 

7.196355622 

19. 

4.8 

6.692770504 

6.905474497 

20. 

5 

6.523014194 

6.644420525 

21. 

5.2 

6.36783515 

6.408735439 

22. 

5.4 

6.225327835 

6.194809651 

23. 

5.6 

6.093911015 

5.999689682 

24. 

5.8 

5.972261036 

5.82093565 

25. 

6 

5.859260982 

5.656514306 
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Table-3.4 

Values of (■Sv.w)eq Corresponding to E as the function of t and T 


equivalent supersaturation 


(X 10-^5ec. ) (X IQ^VJcm . ) 



' w 


1 at 

111 Ctrl.) 

273K 

280K 

290K 

3i0OK 

4.953 

10.023 

9.462 

8.756 

8.145 

5.489 

8.002 

7.597 

7.084 

6.637 

6.000 

6.704 

6.392 

5.996 

5.649 

6.490 

5.806 

5.557 

5.237 

4.956 

6.962 

5.153 

4.946 

4.681 

4.446 

7.419 

4.659 

4.483 

4.257 

4.056 

7.862 

4.272 

4.119 

3.923 

3.748 

8.293 

3.961 

3.827 

3.654 

3.500 

5.593 

7.699 

7.316 

6.831 

6.407 

6.198 

6.308 

6.024 

5.662 

5.344 

6.775 

5.392 

5.169 

4.885 

4.633 

7.329 

4.748 

4.566 

4.333 

4.127 

7.862 

4.272 

4.119 

3.923 

3.748 

8.378 

3.906 

3.775 

3.606 

3.455 

8.878 

3.618 

3.503 

3.355 

3.222 

9.365 

3.384 

3.282 

3.150 

3.032 

6.198 

6.308 

6.024 

5.662 

5.344 

6.869 

5.269 

5.055 

4.780 

4.537 

7.509 

4.574 

4.403 

4.184 

3.989 

8.122 

4.078 

3.937 

3.755 

3.593 

8.713 

3.707 

3.587 

3.432 

3.295 

9.285 

3.420 

3.316 

3.182 

3.061 

9.839 

3.191 

3.099 

2.981 

2.874 

10.110 

3.004 

2.922 

2.816 

2.721 



Critical radius n; (A') 



Supersaturation ratio Sy.w ^ 

fig. 3.8 (a) Variation of Critical radius of water molecules {A') with 
supersaturation ratio Sv.w atT=273K. 
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fig. 3.9 (a) Variation of Critical radius of ice molecules r*(i) with 
supersaturation ratio Sv.i at T=273K. 
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fig. 3.9 (c) Variation of natural log of rate of nucleation ^nJiv 
supersaturation ratio Sy.i atT=273K. 
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Supersaturation ratio Sv.w ^ 

fig. 3.10 (d) Variation of supersaturation ratio Sy.w with equivalent electric 
field at T=288K and T = 5xlO'Sec. 
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Supersaturation ratio Sy.w 

fig. 3.10 (e) Variation of supersaturation ratio Sy.w with equivalent electric 
field Eeq(x lO^V/cm) at T=288K and T = dxlO'^sec. 
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fig. 3.10 (i) Variation of 
field 1 
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—♦—equivalent supersaturation 
at273K 

- equivalent supersaturation 
at280K 

equivalent supersaturation 



0 0.5 1 1.5 2 2.5 3 

E (X lO^V/cm) ► 

fig. 3.11(a) Variation of equivalent Supersaturation ratio (Sv.w)eq 
Corresponding to E (x 1 0^ l^/ cm) as the function of T = 
5x 10’^ sec. at various temperatures. 
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fig. 3.11(b) Variation of equivalent Supersaturation ratio (Sv.w)eq 
Corresponding toE(xlO®V/ cm) as the function of T = 
6x10'^ sec. at various temperatures. 
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fig. 3.11(c) Variation of equivalent Supersaturation ratio (Sv.w)eq 
Corresponding to E(x 10^ V/cm) as the function of T = 
7x 1 0'^ sec. at various temperatures. 
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CHAPTER-4 


EQUIVALENCE BETWEEN TEMPERATURE AND THE 
ELECTRIC FIELD DURING WATER VAPOUR 
CONDENSATION AND GLACIATION 


4.1 INTRODUCTION: 


Atmosphere is filled with a large number of gases and water 
vapour.The amount of water vapour in an atmosphere is constrained by 
the restriction of partial pressures and temperature. On earth the 
condensing substance is typically water vapour, which for ms small 
droplets or ice crystals, typically 0.1 mm in diameter. Above the earth 
surface the temperature varies with height due to thermal fluctuation at a 
certain supersaturation. The water vapour is nucleated to form water 
droplets or ice crystals under certain thermodynamical conditions. In this 
condition the necessary energy of formation of the droplet is provided by 
the phase change of vapour. 

The same energy of formation of the droplet may be provided by 
some other external agency for example-ions, external electric field or 
some foreign material acting as clouds seeds. The phenomenon is called 
cloud seeding. 

In the present study the equivalence between critical temperature 
corresponding to unit nucleation rate per unit volume has been studied 
with external electric field applied to the water vapour. 

4.1.1 Water in the Atmosphere: 

The amount of water vapour present in the atmosphere is a 
complicated function of (1) The amount, which enters the atmosphere 
through evaporation and sublimation, (2) Its transport by motions of 

: v '' ' ^ , ;■ 



~~ • ( thesis of DEVESH DWIVEDl^ 

various scales throughout the troposphere and lower stratosphere, and 
(3) The amount which leaves the atmosphere intermittently and almost 
exclusively as a flux of rain, hail, and snow. The fact that all three phases 
of water contribute to this cycle of the prevailing terrestrial temperature 
and pressures is most fundamentally consequence of the molecular 
structure of water, which permits the strong association of water 
molecules through hydrogen bonding. 

'Fhis is also the principal determinant for the relatively small 
amount of atmospheric water only about 1.3xl0’^ Kg. or about 10'^ of the 
total surface store inspite of the presence of extensive water surfaces for 
evaporation and sublimation. 

Since the earth’s surface is the primary source of water vapour, we 
expect a decrease in the water vapour mixing ratio Wv with height. From 
a slightly different point of view, one may also attribute this expectation 
to the observed decrease in temperature with height in the troposphere. 
Since the maximum possible (saturation) mixing ratio Wvsat decreases 
with decreasing temperature on the average water is squeezed out of air 
parcels during their ascent. Thus Wv is not a conservative quantity in the 
troposphere, but in fact decreases with height. An example of the 
observed vertical distribution of mixing ratio of water vapour in the lower 
troposphere is shown in fig. (4.1) and can be seen to be approximately 
exponential. 


^218 )■ 
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TEMPERATURE »C) 



fig. 4.1 ; Variation with height of the mixing ratio of water vapor in upper 
torosphere at lower stratosphere over Alaska on March 25, 1968. 




4.1.2 The Standard Atmosphere: 
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The balloon surroundings have shown a wide scatter of 
temperatures at different pressures: warmer in summer than in winter, 
warmer near the equator than near the pole, and also changing from day 
to day with the passing whether the systems. But it is possible and of 
much practical value to form same average of the temperature structure 
with height, or temperature large rate. Such averaging was first carried 
out in Europe and later the United States. One called the temperature 
structure so obtained the standard atmosphere (Table 4.1). 

4.1.3 Thermal Layers of the Standard Atmosphere: 

The entire temperature structure of the atmosphere has been 
approximately defined by the 1962 U.S. standard atmosphere. Our 
atmosphere is divided into five layers starting from the earth’s surface fig. 
(4.2). 


1. Troposphere 

2. Stratosphere ' 

3. Mesosphere 

4. Thermosphere 

5. Exosphere 

1. Troposphere: This layer is the most important layer of 
atmosphere. Its average height is 13Km., the air we breathe exists 
here. Almost all the weather phenomenon like rainfall; fog and 
hailstorm occur in this layer. 

2. Stratosphere: Above the troposphere lies the stratosphere. It 


extends up to a height of 50Km., this layer is almost free from 
clouds and associated weather phenomenon making conditions 
most ideal for flying aeroplanes. It contains a layer of ozone gas. 
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3. Mesosphere: It is the third layer. It lies above the stratosphere. It 
extends upto the height of 80Km., Meteorites bum up in this layer 
on entering from the space. 

4. Thermosphere: Here temperature rises very rapidly with 
increasing height. Ionosphere is a part of this layer. It extends 
between 80-400 Km., this layer helps in radio transmission. In fact 
radio waves transmitted from the earth are reflected back to the 
earth by this layer. 

5. Exosphere: This upper most layer of the atmosphere is known as 
exosphere. In this layer gases like He and H 2 float into the space 
from here. 

4.1.4 Water Vapour: 

In the atmosphere, water exists in all three phases, vapour, liquid 
and ice, at different temperature and supersaturation. Phase change takes 
place due to thermal fluctuations. In all phase changes at any temperature, 
a high supersaturation (several hundred percent) is required for 
homogeneous nucleation. In presence of foreign particles, heterogeneous 
nucleation takes place at much lower supersaturation ratio (about 1%) 
(Singh, N., et al., 1985). 

Murino (1979) gave the theory of drop growth in external electric 
field has been extended to nucleation process in water condensation and 
ice glaciations (Singh, N, et al., 1986). In presence of external electric 
field water vapour molecules acquire much earlier the necessary Gibb’s 
free energy for formation of critical nucleus. It has been shown that in 
electric filed induced nucleation; nucleation rates have increased by 
several orders. Inamdar (1998) discussed the effective role of external 




electric field on the relaxation time in nucleation process of water vapour 
condensation and ice glaciations. Pasko (1999) studied about electrical 
analogy applied to modeling dry deposition of particles, occurring soon 
after sunrise has been observed to be enhanced (Marshall, 1999). 

The effect of electric field on water vapour condensation and also 
the polarizability of water vapour molecules have been discussed (Singh, 
N., et al., 2000, 2001). Singh, N. and Kumar (2002) and Sharma et al. 
(1992) have obtained equivalence between applied external electric field 
and the supersaturation ratio. But,, to author’s knowledge the equivalence 
between temperature and the external electric field h as not b een dealt so 
far. 

It is a well-known fact that there exists a vertical electric field in 
free and undisturbed atmosphere. Also, intense electric fields are 
generated in thunderclouds (Sapkota, 1989). The intense electric fields in 
thunderclouds are not widespread and the charge and intense electric 
fields seem to be concentrated in relatively small volumes. Although, 
under actual cloud condition the electric field varies considerably, its rate 
of growth has been calculated by Winn (1974). Moore et al. (1962) and 
Levin and Ziv (1974) reported a raingush following lightning stroke. 
During a lightning discharge maximum electric field (about 10 e.s.u.) is 
produced near the channel. 

Thus, the electric field plays an important role in cloud physics. 
This role could be effectively verified in absence of supersaturation. 
Murino (1979) has shown that a droplet acquires a particular size under 
very low supersaturation. Singh, N. et al. (1986) have shown that in the 
resultant effect on a droplet due to an external electric field and the field 
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induced due to the central dipole. The rate of nucleation in water vapour 
condensation and ice glaciation is about 100 times more near breakdown 
for dry air, as compared to that in absence of electric field. 


In this chapter we have considered the resultant effect due to an 
external electric field and the electric field due to the central dipole. 
These theoretical considerations are applied to the nucleation process in 
water vapour condensation and ice glaciations in order to estimate the 
critical size of the nucleus. 


Following kinetic theory, Collins (1955) inferred that the 
relaxation time is independent of the free energy of formation of the 
nucleus, but it varies as the square of the radius of the critical nucleus. 


4.1.5 Structure of Water Vapour: 

Experiments indicate that water molecules in water vapour tend to 
interact and form crystals, in contrast to ideal gas behaviour dimmers as 
well as higher order polymers are considered to be present in water 
vapour, do in small concentration only. The experiments involving 
molecular beam techniques (Lin et al., 1973; Searcey, 1974) suggested 
that in oily supersaturated water vapour clusters of upto 180 water 
molecule may be present. Clusters of 21 water molecules seemed to 
exhibit particularly large stability. 21 water molecules can be arranged in 
the form of a pentagonal clode cahedron with a molecule at each comer 

and a single molecule in the center of the ‘cage’. 


Studies on the possible and more likely clusters types have been 
reviewed, by Rao (1972) and Kell (1972 a). Recent theoretical studies of 
the formation of water clusters have been carried out (Kistenmacher et 
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al., 1974a, 1974b; Abraham, 1974). Kistentenmacher et al. (1974a, 
1974b) found two possible stable configuration for the dimmers, a cyclic 
form and an open form which was more stable. For the trimmers and 
tetramers the cyclic form seemed to be somewhat more stable than the 
open structures. For the large cluster the authors suggested not a single 
structure but a statistical distribution of different configurations, since 
many configurations with significantly different geometry were found to 
passes nearly the same energy. 

The Potential energy of interaction, U, between a pair of water 
molecules has the general character of being strongly repulsive at very 
close separations and weekly attractive at longer range. One widely used 
and relatively simple expression for it is due to Stockmayer (1941). 


U = 


fi^f c ccr^^ 
r24 


(4.1) 


Where r is the separation of the molecules; ju, is the dipole moment 
of an isolated water molecule; a, is the collision diameter (the molecular 
separation at which U=0 at {i = 0); c is an adjustable constant and f is a 
known function of the mutual orientation of the two molecules. 


The first term on the right side of eq. (4.1) is just the dipole-dipole 
contribution of the interaction energy, and may be attractive or repulsive, 
depending on the dipole orientations. The second term represents 
contribution firom: (1) the interaction energy between a permanent dipole 
of one molecule and the dipole it induces in the other (dipole-polarization 
or induction interaction), (2) the net energy arising fi-om momentary, 
fluctuating dipoles interacting with the corresponding induced dipoles 
(polarization- polarization or dispersion interaction). Even though the 
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time average of these dipole fluctuations may be zero, the energy 
contribution is proportional to their mean square, which is finite and 
positive. Both (1) and (2) are usually as Vander Wall’s interaction, which 
by its nature can be seen to bring about an attractive force between the 
molecules. The third term in eq. (4.1) represents the short-range repulsive 
forces, which loosely may be ascribed to the overlap of electronic orbital, 
which are incompatible according to the Pauli Exclusion Principle. 

There is a little doubt that the Stockmayer potential or similar ones, 
such as (Rowlinson, 1949, 1951) potential, portray with fair accuracy the 
interaction between pairs of water molecules at large separations in dilute 
water vapour (Singh, N. et al., 2002) . This is evidenced by the fact that 
values for the second viral coefficient computed via eq. (4.1) can be made 
to fit experimental values. On the other hand, the same potential functions 
yield values for the third viral coefficient of water vapour, which disagree 
substantially with experiment. Partly, this is due to the approximate 
nature of eq. (4.1), and partly because three-body interactions should be 
included also, since other molecules in the system can significantly 
modify the interaction of a given pair. In particular, the Stockmayer 
potential is insufficiently “directional” in character to account for the 
geometry of cluster formation in water vapour. A recent more 
complicated potential function, which has proven to be predictive value 
in this respect, is described briefly. 

4.1.6 Vapour Pressure Over Pure Droplets: 

It is found from measurement that the vapour tension of a concave 
surface of a liquid is lower than that of a flat surface. Convex surfaces 
show higher vapour tensions than flat surfaces. Droplets therefore require 
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a certain degree of supersaturation in the ambient space in order to remain 
in equilibrium. Physically, one might indeed expect the escaping 
tendency of molecules from a spherical droplet to be greater than the 
smaller droplet because the component of the binding force in the 
tangential direction is diminished. Lord Kelvin (Wiliam Thomson) first 
deduced the relation, which today can be presented in more modem 
thermodynamic from the point of view of the surface free energy. 

In the supersaturated state the transition surface lies above that 
depicted in plate 2.1, each isotherm being at a higher pressure and 
therefore at a higher free energy. The elevation of the free energy above 
that for a plane surface of pure water depends on the saturation ratio and 
on the surface free energy. 


The problem is approached by considering a system initially 
consisting of vapour at supersaturation, with saturation ratios, which 
changes of final state, and forming from the vapour a single embryonic 
droplet of radius r containing ‘g’ molecules. Note that 


4 3 « 

g = -nr^Zi 


(4.2) 


Where Zl is the number of molecules in a cubic centimeter at the 
liquid, and that Zxr^pjm^ where N is Avogadro’s number 
(6.022x10^^). Its molecules are considered in place of moles, Nk is 
substituted for R, where k is Boltzmann’s constant (1.381 xlO^Srg K ^ 
molecule'^ or effectively the gas constant per molecule). The free energy 
differential dG of equation becomes 


dG = gkTdiinp) -^nr^ZikTdi-inp) 


(4.3) 


- ( 228 > 
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As a molecules go from the vapour at pressure p to the liquid at 
tension po the change in chemical potential or free energy of the system is 
given, after the manner of equation as, 

t^Po-l^p = Ip° 9 ^Tdi-enp) = GkT 

= -{nr^Z^kT-en^ = nr^ Z j^kT dinS (4.4) 

Where S is now the saturation ratio. As the transition is made from 
vapour only to an embryo droplet of radius r, the surface free energy goes 
from zero to 4nr^a. The elevation of the total free energy of the system 
of g molecules in the presence of these two effects is therefore their sum, 
or 


AG = 4nr^ - J nr^Zj^kTd-enS (4.5) 

The significance of this expression may best be appreciated by 
plotting AG against ‘r’ for several values of S at constant temperature, 
and therefore also, constant a (fig.4.3). 

4.1.7 Equation of State for Water Vapour: 

Unlike other atmospheric constituents, water appears in all three 
phases, solid, liquid and vapour. In the vapour phase water in the 
atmosphere behaves as an ideal gas to a good approximation. Its eqation 
of state is 

(4.6) 




B — pv^v^ 
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Where e, is vapour pressure; is vapour density; and R^, is individual 
gas constant for water vapour (461 JKg-'K'’). This equation sometimes 
appears in the form 

rr, 

e-Pv~T (4 7 ) 

Where£ = - = ^ = 0.622 

Rv rn 

Water vapour has a specific heat at constant pressure equal to 1850 
JKg''K'‘. The specific heat at constant volume equals to 1390 JKg''K'‘. 

4.1.8 Equilibrium Vapur Pressure e^ and latent Heat 
L Over a range of Temperature: 

Values of Cs and L over a range to temperature are given in Table 
4.2. At temperature of 0®C and colder the latent heat of sublimation and 
the equilibrium vapour pressure ei over ice are also tabulated. 

4.2 IHEOeETlCAL CONSIDERATIONS: 


4.2.1 Equivalence between Electric Field and Temperature: 

In presence of an external electric field, the rate of increase of 
radius with respect to time is given (Singh, N., 1 986) as 




(4.8) 


Where is the density of water vapour; p^ the density of water; 
cc the polarizability; A, the mean free path of water vapour molecules; E 
the applied external electric field; m^, the mass of water molecule; 
the radius of water droplet. 
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TABEL- 4.2 


Latent heat of condensation (sublimation) and saturation 
vapour pressure for various temperatures 
(Source, Smithsonian Meterological Tables) 


T(®C) 

e* (mb) 

e; (mb) 

L (calg-') 

Ls (calg‘) 

l| 

■ 

'f 
^ j: 

i 

' ''1 

: I 
. 1 

-40 

0.189 

0.128 

621.7 

678.0 

-35 

0.314 

0.223 



-30 

0.509 

0.380 

615.0 

678.0 

-25 

0.807 

0.632 



-20 

1.254 

1.032 

608.9 

677.9 

-15 

1.912 

1.652 



-10 

2.863 

2.600 

603.0 

677.5 


603.0 


677.5 
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Integrating eq. (4.8) within limits nv = 0 to r*, the critical radius 
of nucleus and t = 0 to t = Tyy the relaxation time, we have 

r* = (4.9) 

Following liquid drop model, in homogeneous process of water 
condensation, the nucleation rate (Pruppacher, 1978) (i.e., number of 
critical nuclei per unit volume per sec.) is 

Where is condensation coefficient ; density of water; N, 
Avogadro’s number; Mw, molecular weight of water; <Jw/v, surface 
energy of water- vapour interface; Csat.vv the saturated vapour pressure 
over water; R, the universal gas constant; Sp_w th® supersaturation ratio of 
water vapour over liquid water; AG^, the Gibb’s free energy for 
formation of critical nucleus; k, Boltzmann constant and T, temperature 
of water vapour. 

The pre-exponential factor has been calculated by Pruppacher 
(1978) for ac=l is 9x10^' cm'^ sec'* and 3.7x10^^ cm'^sec'* at -30°C and 
30°C, respectively. Although prefactor changes with temperature and 
supersaturation ratio, but it has been assumed constant for present 
calculations, because the nucleation rate is mainly dependent on the 
exponential term. 


Hence, eq. (4. 1 0) becomes 
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For critical nucleation ratej^=l, hence 

The value of AG^ in homogeneous nucleation process is 


AG = 

Subsisting eq. (4. 13) in eq. (4.4), we have 

(fn s,,„y = 1 ^ 4 /. 

R 1 pl.x3k£n9xl0^^ 


(4.13) 


(4.14) 


For critical nucleation rate (J - 1) we use critical temperature T = 


Tc hence 


167rMwcrj^,yry 

{■in Sp,wf p^x3kin9xl0^^ 


(4.15) 


Alternatively, radius of critical nucleus in homogeneous 
nucleation with supersaturation is 


p\^RT fti Sp^ 


(4.16) 


And in presence of external electric field, the same critical radius 
is obtained in electric field E and relaxation time as 


^ V^Pv'^N j 

l'2aXE^\ 

1 2Pw ' 

\ J 


172^2/3 


(4.17) 
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Equating eq. (4.16) and (4.17) at critical temperature T=Tc, we 

have 

K 

Where constant K 

Eq. (4.18) gives equivalence between critical temperature Tc and 
the applied external electric field. 

4.3 RESULTS AND DISCUSSION: 

From eq. (4.18), the values of critical temperature Tg have been 
computed at different values of supersaturation ratio for a given value of 
external electric field and relaxation time. At E=2,3,4 e.s.u. and T/v = 2 x 
10“^ sec. The value of Tg at different values of supersaturation ratios are 
given in Table (4.3). 

From the table, it is evident that at a given value of electric field, 
the values of critical temperature decrease with increase in super 
saturation ratio. For example, at an external electric field of 4 e.s.u. and 
relaxation 2x10'^ sec., the critical temperatures are 2.0990xl0'*K, 
0.814271 xloX 0.552134538 x 10Y corresponding to the supersaturation 
ratios 1.2, 1.6, 2. Simil^ly the values of critical temperatures have been 
calculated at E=2,3 e.s.u. also and are given in table (4.3). 


rwy 
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TABEL- 4.3 


Calculated values of critical temperatures for different 
supersaturation ratios at Relaxation time TA^=2xlO'^sec. 


^v.w 


1.2 


1.4 


1.6 


1.8 


2.0 


2.2 


2.4 


2.6 


2.8 


3.0 


3.2 


3.4 


3.6 


3.8 


4.0 


E=2 e.s.u. 


3.332107614 


1.805542870 


Te(x 10^) K 


E=3 e.s.u. 


2.542874029 


1.292575227 


1.033563849 


0.876458947 


0.770510973 


1.377887092 


0.986419515 


0.788756839 


0.668863359 


0.588009923 


0.693931173 


0.635800911 


0.590038453 


0.552984027 


0.522301186 


0.496426903 


0.474275116 


0.455067088 


0.438229473 


0.529568598 


0.485206906 


0.450283614 


0.422005795 


0.398590405 


0.378844631 


0.361939653 


0.347281185 


0.33443168 


E=4 e.s.u. 


2.099096262 


1.137420734 


0.814271368 


0.651104425 


0.552134538 


0.485391497 


0.437149246 


0.400529476 


0.371700933 


0.348358108 


0.329029129 


0.312729352 


0.298774601 


0.286674302 


0.276067269 


i, 
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4.4 CONCLUSION: 

The study of equivalence between external electric field and the 
critical temperature suggests that in electric field induced nucleation 
process, the critical temperature required to critical nucleation = 1) 
decreases with increase in electric field as compared with homogeneous 
nucleation. Thus, the nucleations of water condensation increase very 
much due to applied external electric field. Hence, the electric field is 
very much effective in nucleation process. A small external electric field 
is equivalent to very large value of supersaturation ratio, which, 
otherwise, never exists in the clouds. In electric field induced nucleation, 
the critical nuclei of smaller sizes are formed having smaller number of 
water molecules and require less Gibb's free energy for their formation 
and equilibrium concentration of critical nuclei is enhanced as compared 
with the homogeneous nucleation. Also, from the equivalence of 
supersaturation ratio with external electric field, it appears that electric 
field is more efficient than supersaturation. Although, there are a number 
of mechanisms giving possible explanation for raingush after lightning. In 
present study, we may add possibility to explain the raingush 
phenomenon. A high electric field produces lightning discharge, creating 
ionization of the gaseous medium and a large number of ions are formed. 
Hence, for a given amount of liquid water a large number of droplets of 
smaller size are formed immediately following lightning. Hence, 
lightning seems to be cause of raingush and not the effect. 
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I CHAPTERS 

EFFECT OF IONS ON ELECTRIC FIELD INDUCED 
NUCLEATION PHENOMENON AND STUDY OF 
THUNDERSTORMS 


5.1 INTRODUCTION: 

Gagin (1972) and Huffman (1973) found that at any given 
temperature, the number concentration of ice nucleating (IN) particles 
increase with increasing relative humidity and co-relates with 
supersaturation over ice according to the relation 

Nin = CS«i (5.1) 

Where is the supersaturation ratio of water vapour over ice; C 
and K are arbitrary constant. 

Considering that the aerosol particle (AP) concentration decreases 
with increasing height, one might expect the same from concentration of 
IN. However (Bigg, 1961, 1963), found that in air over Australia the 
concentration of IN was significantly higher at 13 to 27 km altitude then 
at the ground Analogously (Rosinski, 1967), found over Colorado that on 
certain occasion the IN concentration was highest near the jet stream. If 
one compares the concentration of IN and cloud ice^ particles, it is found 
in many clouds, particularly at relatively warm temperatures, the 
concentration of ice particles may exceed by many orders of magnitude, 
the concentration of IN determined at the cloud top temperature 
(Pruppacher, 1978). 

— rsn,. — 


5.2 IONS, THE POTENTIAL CONDENSATION NUCLEATION 


CENTERS: 

Studies ot ion-induced nucleation were first made by Wilson 
(1897). Since the time of Wilson, numerous research papers have been 
devoted to the subject (Castleman, 1979a, 1979b, 1972; Suck, 1982; 
Varshneya, 1969, 1971) ion induced nucleation. 

The ion provides a central force field, which is a simple form of 
attractive potential resulting in the formation of heterogeneous embryos, 
A basic understanding of ion induced nucleation therefore provides an 
important basis for further developments in the theory of heteromolecular 
nucleation, which is important in atmospheric processes (Allen, 1969; 
Doyle, 1961). 

Below saturation the chemical bond energy of hydrated ions are too 
weak to permit the growth of critical embryos. The prominent forces 
resulting in the stabilization of small clusters about simple ions such as 
those isoelectronic with the noble gas structures are electrostatic in order 
to assess the importance of electronic structure on bounding investigation 
(Castleman, 1972) on metal ions (e.g. Pb"^ which have non-closed shell 
electronic configurations, were compared with alkali metal series (Dzidic, 
1970). 

Condensation of water vapour of ice germ formation in super 
cooled liquid, when ionizing radiations are passed (Vashneya, 1969), is a 
milestone in the understanding of ions as centers of condensation. This 
basis has been confirmed theoretically (Varshneya, 1971). The pioneer 
worker (Castleman, 1979; Mohnen, 1970, 1974; Chan, 1980a, 1980b) has 
placed the belief on a firm footing. 
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5.3 ION SPECIES IN TEOPOSP!i£ft£; 

The sequence of chemical reactions that alters the ion species 
between their creating by cosmic rays or radioactive decays and their 
removal by either ion-ion recombination, attachment to a particle, or 
some other loss mechanism is quite complex. 

Ion lifetimes in clean air are expected to be the order of 100 sec 
(Ferguson, 1979). However, during the first one-second interval after 
ionization, most of the ions in the troposphere have experienced enough 
collisions with neutral species, whose concentrations are in parts per 
billion (ppb) range or higher, to have established very stable core ions 
due to attachment. 

Perkins and Eisele (1984) reported first mass spectrometric 
measurements of naturally occurring positive and negative ions at ground 
level under a variety of weather conditions. Positive ions were observed 
with masses 54, 60, 80 and 94 amu, mass 54 amu is believed to be 
NH^. 2 ( 7 / 20 ) the negative ions were observed to have the masses 60 and 
125 amu and were believed to be IVO 3 and NO^ (HNO 3 ) . Several ions 
that appear to be H 2 O cluster of the above ions were also observed. 

Mohnen (1970) reported the terminal ions in troposphere as being 
hydrated 0^ and CO 4 .The degree of hydration depends upon 
temperature and the mean number of water molecules attached to these 
ions. As soon as the water vapour concentration is less than the 
concentration of carbon dioxide (above 8km). The terminal ion is found 
to be hydrated . 
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Using an improved aircraft borne mass spectrometer with high 
mass resolution (Hauck, 1984), reported positive ion composition 
measurements in the lower stratosphere and upper troposphere (Table 
5.1). 

The ions 37±1, 55±1 and 71±1 may readily be identified as proton 
hydrates (//+(//2 6>)„;n = 2,3,4) which are expected (Arnold, 1982) to 
form quickly at the tropopause from primary positive ions produced by 
galactic cosmic ray ionization. 

Legrand and Delmass (1984) analyzed 100 snow samples representing a 
continuous sequence of 10 years of precipitation at South Pole. The results 
show that the ionic budget of deposited smoke at this central Antarctic location 
is satisfactorily balanced when the concentration of eight major ions 
,Na^ ,K'^ ,NH'^ ,SOl~ .NO^ ,Cl~) are taken into account. 
Na^ ,K'^ ,Mg^^ and Cl~ are typical constituents of sea salt particles 
which form one part of atmospheric aerosol over Antarctica, but 
NH^ excess SOl" and excess CT are linked to an important 
acid aerosol fraction which is most probable gas derived. 

Theory of ice formation charged particles has been supported 
indirectly by many workers. Bigg (1961) and Bigg et al., (1963) found 
experimentally that IN concentration is sufficiently higher (~10^) in upper 
and lower atmosphere than at ground. The number concentration of cloud 
ice particles have been observed in excess than the number concentration 
of IN at cloud top temperature. The above findings can be explained if we 
consider the effects of ions on condensation and ice nucleation. 
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Table- 5.1 

Positive ion masses obtained from high-resolution spectra at 


11.28 km altitude (Hauck, 1984) 


Mass 

Tentative Ion IdentiHcation 

28 ± 1 

HCO^ 

H^HCN 

32 ± 1 

H^CHsOH 

37± 1 

H*(H20)2 

50 ± 1 

H^CHsOH.HjO 

55 ± 1 

H"(H20)3 

59 i: 1 

H"(CH3)2C0 

64 ± 1 

H^'CHsCOH.HsO 

H^C 2 H 50 H.H 20 

irCH 30 CH 3 .H 20 

H11C00H.H20 

73 ± 1 

H"(H20)4 

77± 1 

ir(CH3)2C0.H20 

88 ± 1 

H^CH 30 H.(H 20)2 

144 ± i 

r(CH3)2C0.(H20)3 

126 ± 1 

H*(H20)7 
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However, Detweiler and Vonnegut (1980) did not find any effect 
of ions on condensation process. Their estimate of the concentration of 
ice particles cannot be taken to be very reliable. Recently Dwivedi D. et 
al., (2009) studied the variation of relaxation time in water vapour 
condensation and ice glaciation as a function of supersaturation ratio in 
presence of ions. 

5.4 IONS IN THE ATMOSPHERE: 

Near the ground a, |3 and y radiations from Thoron and Radon 
decay and cosmic radiation are the principle ionizing agents causing 
charge separation initially into electrons and singly charged positive ions. 
At altitudes up to 60 km cosmic ray background is still the major 
ionization source; if we neglect disturbances known to occur after high 
solar activity. 

The fi-ee electrons will, after thermalization attach to neutral 
molecules. The positive and negative ions now undergo a series of so- 
called ion-molecule reactions covering simple charge transfer of the type 
(Mohnen, 1974). 

A-^ + B-^A + B* 

A~ -h B A B 

or attachment reactions (^[=3*^^ body) 

A^ B M A"*" . B ~h M 


A~ + B + M A~.B + M 



or switching reactions 


. B C —> . C -{■ B 

A-.B + C A-.C + B 
or proton transfer reactions 
H'^.A + B -^A'^ +H.B 
or clustering reactions 

A-^.B + B + M A+.B 2 + M 
A-^.B2 + B + M A-^.Bs +M 
And 

A^.B + B + M t;A-^.B2+M 

Each ion-molecule reaction is characterized by a rate constant. 

5.5 DESmUCmON OF SMALL IONS IN THEATHOSPiffiHE: 

Particles in the atmosphere may be formed through the coagulation 
of two or more particles. But this process does not produce any new 
particles. Infact it decreases the total number of particles. New particles 
are produced in the atmosphere through a number of mechanisms. The 
wind stresses and other natural or artificial mec^ipii'cal stresses can 
detached the particle from the surface. For example, fragments can be 
detached from solids through abrasion cracking, weathering, etc. and 
from oceans through breaking of waves, bursting of bubbles. 
Supersaturated vapour can produce solid (i.e. ice) and liquid elements 
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directly from the water phase by nucleation. In addition, human activity is 
another effective producer of particles. (Table 5.2 gives a general picture 
of existing particles in the atmosphere.) 

In general the concentration of large nuclei over oceans is one 
order less than that of over land. About half of the small condensation 
nuclei carry a net electrical charge. 

The small ions (created by the atmospheric radioactivity and 
cosmic rays) will be destroyed by the attachment and uncharged large 
nuclei. When a small ion gets attached to uncharged large nuclei, it will 
become a large ion. 

5.6 AK INTRODUCTION ON ATMOSPHERIC ELECTRICAL 
CONDUCnVtiY: 

Deshpande and Kamra (2002) reported the atmospheric electrical 
conductivity measurements made over the Indian Ocean with a Gerdien’s 
apparatus mounted MB polarbird during XVI Indian Scientific 
Expedition in 1996-1997. They also reported simultaneous three hourly 
measurements of aerosol concentrations of 13-1000 n.m. size and some 
meteorogical parameters. Latitudinal variation of conductivity along the 
Cruise route show a minimum at ~28® S. Further, the variation in the 
conductivity in 10‘^-20‘^S and 60°- 70®S Latitudinal belts show opposite 
trends on the outward and return cruises, which fall near to the reset and 
withdrawal phases of the north-east monsoon season. 
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Table -5.2 


Some Characteristic Features of Atmospheric Particles 


I 

i 

Type 

Concentration 

no/m^ 

Range of Sizes 
Diameter in 
Meters 


Gas 

molecules 

Almost 25x10^*^ 

(2.8-6.5)xl0''° 


Small ions 

(l-7)xl0' 

(1.5-1.0)xl0''' 


Large ions 



Small Aitken 

(lOLlO") 

(l-40)xl0'** 


Approximate 
Terminal 
Velocity m/sec 


Condensation 
Nuclei 
Large nuclei 


Giant nuclei 


Dry haze 




Droplet 


Drizzle 
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From these measurements in the free troposphere of the equatorial 
pacific, Clarke (1993) showed that nucleation in the upper troposphere 
and subsequent growth by coagulation and vapour diffusion could be a 
major source of new particles in that region. These measurements support 
the conclusions of a large scale model of aerosol nucleation and dynamics 
in the free troposphere and marine boundary layer developed by Raes 
(1995), which suggests that injection by mixing with the free troposphere 
is a major source of new sub-micron particles in the marine boundary 
layer. For interpreting their mid-pacific formation occur in outflow 
regions of the frontal or convective clouds in the free troposphere where 
conditions of new particle formation are met (Raes, 1993). These newly 
formed particles grow by coagulation and condensation. Transport of 
aerosols from the free troposphere to the marine boundary layer occurs by 
convective mixing in the vertical and general circulation in the horizontal. 
These particles grow by vapour deposition and therefore the growth is 
faster in the tropical marine boundary layer where gaseous precursor 
concentrations are large. 

Their measurements of the atmospheric electric conductivity along 
the cruise track during the months of December and March do indicate a 
secular decrease in the conductivity which is most likely a consequence 
of an increase in the background air pollution. 

5.6.1 Diurnal Variation of the Atmospheric Electric 
Field and Conductivity at Maitri, Antarctica: 

Deshpande and Kamra (2002) made surface measurements of 
atmospheric electric conductivity at the Indian Station Maitri (70 
45’52”S., 11° 44’ 3”E, 117 m above mean sea level), Antarctica from 
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January 10 to February 24, 1997. The total conductivity is 2.1x10'*'* S.m'* 
and does not show any significant diurnal variation fig.(5.1) shows the 
diurnal variation of hourly averaged values of conductivity of both 
polarities averaged for 20 fair weather days at Maitri. Where the 
convention was that the fair weather field pointing downward is negative. 
The electrical conductivity of both polarities does not show any 


significant variation except for somewhat higher values during the 1500 
to 2400 UT period. The conductivity of positive polarity is generally 
always higher than that of negative polarity. The mean values of positive 
and negative conductivity are +1.2 and -0.9 Sm'*, respectively. The 
difference may possibly be due to the manifestation of the electrode 
effect. Kandalgaonkar et al., (2003) analyzed the data of point discharge 
current (PDC) measurements during a total of 65 thunderstorm of Pune 
(180° 32’N, 73°, 5rE) to study the (PDC) local diurnal variation. 

5.6.2 Variation in Atmospheric Electrical Conductivity 
Beneath A Thundercloud: 

Singh et al. (1985) has developed an expression for atmospheric 
electrical conductivity on the assumption that thunderclouds produce 
intense electric field at the ground surface creating large scale ionization. 
Atmospheric electrical conductivity has been found to remain unchanged 
unless the ground electric field exceeds a certain threshold value to 
initiate corona near the ground surface. Effect of height of thunderclouds 
on the atmospheric conductivity has also been discussed. 


Considering dimensions of thundercloud (Singh, N., 1984; 
Mathpal, 1983; Singh, N., 1985) D=l. 25km., L=1.50 km., H=0.50 km. 
ttip hf^tinvimir at atmosohcric electrical conductivity beneath a 
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UT 

fig. 5.1 : Diurnal variations of the potential gradient, the conductivity of 
both polarities and the air-earth current averaged for 20 fair 
weather days at Maitri. Continuous curve with solid circles 

shows the Carnegie curve for the diurnal variation of potential 
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thundercloud have been determined. The calculations have been 
performed assuming = Eom = 100 Vm'* , q = 10x10 PI m^ sec'’, 
N=1.5xl0'“’ m'^ Tj = 4.6 xlO'*^ m'^sec', q o=1.8xlO’‘^ sec"' and a = 
1.6 xlO m^ sec '. The results for various concentrations of large ions 
have been shown in fig. (5.2). This figure shows that atmospheric 
electrical conductivity remains unchanged unless ground electric field, 
Bout exceeds its threshold value to initiate corona near the ground surface. 
As ground electric field exceeds its threshold limit a sudden increase in 
the conductivity has been noted. Fig.(5.3) shows the nature of 
conductivity variation for different cloud base heights from the ground 
surface. A lower thundercloud produced electric field with higher growth 
rate at the ground surface, which involves into higher maximum value of 
ground electric field. The threshold value of ground electric field to 
produce corona ions also reaches earlier in case of lower clouds. 
Consequently lower base thunderclouds affect conductivity to the large 
extent fig. (5.3). 

A thundercloud constitutes a rapid and widespread perturbation in 
the electric field at the ground surface. Anderson (1977) also pointed out 
that investigations of the atmospheric phenomenon near the ground must 
consider the electrode phenomenon (increased positive ions density) 
because almost all measurements near the ground will be affected by it to 
some degree therefore, taking this aspect into account efforts have been 
made to gain a physical insight into processes controlling the electrical 
state of the atmosphere near the ground. 
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fig. 5.2 : Variation in conductivity for various concentration of 
large ions. 
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fig. 5.3 . Variation in conductivity by different cloud base heights 
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The electrical aspects of atmosphere near the ground surface are 
principally characterized by the different conditions of the ionization 
equilibrium existing near the ground. The general ionizing agencies, 
which are supposed to produce ions in the lower atmosphere, are: (1) the 
penetrating radiation from the radioactive substances in the earth crust, 
(2) radiation from the radioactive dust in the atmosphere and (3) cosmic 
rays. 

The statistically inhomogeneous Markovian model and the 
stochastic radiative transfer equation have been used to derive equations 
for the mean radiance of solar radiation. It was demonstrated that in 
extreme cases the obtained equation agree with corresponding equation 
previously derived for (i) the statistically homogeneous broken clouds 
and (ii) the vertically inhomogeneous overcast clouds (Evgueni,2003). 

5.7 CONTRIBUTION OF CLOUD AND PRECIPITATION PARTICUES 
TO THE ELECnaCAL CONDUCTlVtlY AND THE 
RELAXATION TIME THE AIR m THE THUNlKRSrORMS: 

Kamra (1979) calculated electrical conductivity caused by the 
charged cloud and precipitation particles in thunderstorms. It is found that 
the conductivity may be 1 or perhaps 2 orders of magnitude higher than 
the clear conductivity at the same altitude. Since the conductivity 
increases with the precipitation intensity, liquid water content, electric 
field and the electrical charge on the particles. 

The value of electrical conductivity within thunderstorms is a 
matter of considerable dispute. Since small ions produce by cosmic rays 
should soon be captured by the cloud droplets. Gunn (1954), Phillips 



(1967), and many others have suggested that the conductivity of cloudy 
air should be much less than that of clear air. 

Incorporating polymerase chain reaction based identification, 
population characterization and quantification of microorganism into 
aerosol science is studied by Jordan Paccia et al. (2006). 

There are many sources and types of charge carriers in 
thunderstorms. Small ions are produced by cosmic rays and are further 
supplemented by point discharge current flowing from the ground and the 
conduction currents flowing from outside the cloud. Hydrometeors 
themselves may produce small ions by corona in intense electric fields. 
Further, the cloud and precipitation particles may themselves become 
charged by capturing these small ions or various processes of charge 
generation in thunderstorms. 

A unique feature of conductivity caused by cloud and precipitation 
particles is that in spite of its high values, the corresponding relaxation 
time of the cloudy air for the existence of space charges may not 
necessarily be small, this being specially so far small pockets of space 
charge and for the individual charged particles. Growth of large volumes 
of charge producing intense electric fields may, however, be significantly 
influenced by dissipation currents caused by this conductivity. This effect 
is now recognized and included in the calculations of the ambient electric 
field caused by the growth of the main charge centers of the 
thundercloud. (Kamra, 1970, 1971, 1975; Mason, 1972; Paluch, 1973; 
Ziv, 1974; Scott, 1975; Griffiths, 1975) 
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6.6 VtolATlONS IN ATM0^4ffiRlC AEROSOLS AND 
ELECTRICAL CONDUCTTVllY AT ROORKEE DURING THE 
TOTAL SOLAR ECLIPSE OF OCTOBER 1995 

Singh et al. (1999) made measurements of some meteorological 
and electrical parameters with the total solar eclipse of October 24, 1995 
continuously from October 22-27, 1995 at Roorkee (29*’ 52’ N, 77° 53’ 
52” E, 275m above sea level). Roorkee observed (90-92%) maximum 
obscurity of eclipse, as it was close to the path of totality. The even lasted 
from 7.10 am to 9.30 am (1ST; 0140-0400 hours GMT). Considerable 
increase of aerosol and positive and negative conductivities were 
recorded during the eclipse with respect to those made on any other day. 
Increase in relative humidity and decrease in temperature were also 
observed during the eclipse. 

5.8. 1 Measurement of Electrical Conductivity: 

The technique for the measurement of atmospheric conductivities 
makes the use of the Gerdien condenser (Mukku, 1982). This instrument 
fig.(5.4), has been used widely for the measurement of electrical 
conductivity, ion mobility and density. It has been used for balloon-bom 
and ground-based measurements. The electrical conductivity of the air 
can be determined from current- voltage characteristics of the condenser. 

A positive potential at the outer electrode makes electrons and 
negative ions move to the central electrode, thereby constituting a current. 
This gives the negative conductivity. The negative potential on the outer 
electrode similarly provides us the positive condudcivity. The total 
conductivity cr, in terms of positive (oj,) and negative (an) conductivity, 
is given by 


<^ 4 > 


(7 = 0-p + <T„ 


(5.2) 
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From the measurement of output current of EMA fig. (5.4), the 
positive and negative conductivities can be obtained from (Mossop, 1955) 

ip ~ (5.3) 

Where, 

Wp - Mobility of positive ions 
Wc - Critical mobility of ions 


M - Aspiration rate 
e - Electronic charge 
rip- Number of positive ions 
(7p- Conductivity of positive ions 
Also, 




IrcLV 


(5.4) 


Where, 

Rq - Radius of outer electrode 
Ri- Radius of inner electrode 

V - Voltage applied across the electrodes of the capacitor 
L - Length of cylinder 
Similarly, one may write 


MeinjiWn+neWe) 


‘■n 


We 


M(Tn 

We 




(5.5) 



B = Battery 

E = Electrostatic shield 


EMA = Electrometer Amplifier 
MCR = Multichannel Chart Recorder 
P = PercePQX 
T = ration 

fig. 5.4 : Experimental set-up for conductivity measurement. 
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Where is the current due to negative ions and electrons, 71 ^ the number 
of negative ions and Ug that of electrons, and Wg are the mobilities of 
negative ions and electron, respectively. 

• ^*^ 3^ prese nt case, si nce the applied voltage was the same in each 
case (±30V), the critical mobility is given by 

Wg = 2.6624 X (5-6) 

5.8.2 The Observed Variation in Electrical Conductivity: 

Singh et al. (1999) observed the variation in electrical conductivity 
of atmosphere fig. (5.5), at the time of eclipse. The variation in negative 
conductivity is found to be larger than for positive conductivity. The 
duration of peak negative conductivity is larger than for positive 
conductivity. Table 5.3 shows the measurements of atmospheric 
conductivities by other workers during various solar eclipses. The 
observed variations in positive and negative conductivity are in good 
agreement with the present results. During the solar eclipse, both the 
aerosol concentration and electrical in which the increase of aerosol 
concentration normally occurs with the decrease in atmospheric electrical 
conductivity. 

Nitrogen budgets of agriculture fields of the Changjiang River 
basin from 1980 to 1990 (Xiang Bao, 2006).The increased ionization 
during the solar eclipse is responsible for the increase of both positive and 
negative conductivities. Singh et al. (1985) have shown that the 
atmospheric ions increase the cloud condensation nuclei (CCN) 
formation. Singh, N., Ph.D. Thesis, (1985) has shown that the 
atmospheric ionization plays an important role in aerosol formation. 



Conductivity xIO ohm 
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fig.5.5 ; Time variation of atmospheric conductivity (a) tor positive 
conductivity ond (b) for negative conductivity. 
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TabIe-5.3 


A Summary of Measurements of Electrical Conductivity by Various 

Investigators 


Solar Eclipse 
period 

Investigator (s) 

Place 

Remarks 

7 March 1970 

Anderson and 

Dolezalek 

(1972) 

Waldor Annex 
of the Naval 
Research Lab., 
Washington 

DC (USA) 

Clear decrease in 
negative conductivity 

10 July 1972 

1 

j 

1 

Lane-Smith and 
Markson (1977) 

Malignant 

Cove Nova 
Scotia (USA) 

The negative 
conductivity began 
to rise up to High 
values about 5 min, 
after the totality 

16 Feb 1980 

Nizamuddin et 
aL(1981) 

Nagarampalam 

Visakapatnam 

(India) 

Both types of polar 
conductivities 
increase but increase 
in positive 
conductivity is more 
pronounced than 
negative conductivity 

18 March 1988 

Dhanokar 
Deshpande and 
Kamra(1989) 

Pune (India) 

Remarkable changes 
in atmosphere 
conductivity during 
the period of eclipse 



24 Oct 1995 


Singh et al. 
(1999) 


Roorkee 

(India) 


Increase in both 
positive and negative 
conductivity 



( thesis OF DEVESH DWlVEDl) 

The increased ionization would increase the concentration of 
aerosols. Thus, the atmospheric conductivity and aerosol both has been 
found to be increased during the total solar eclipse because of the excess 
atmospheric ions created by galactic cosmic rays. Forbush decrease 
enhanced by obstruction of solar wind particles by the moon. 

5.9 USING ATMOSPHERIC ELECTRICAL CONDUCTIVITY AS AN 
URBAN AIR POLLUTION INDICATOR: 


Guo et al. (1996) investigated the atmospheric electrical 
condectivity as an air, pollution indicator through correlation analysis 
based on two year of daily averaged observations of the atmospheric 
electrical conductivity; weather classification index, and the 
concentration of aerosol partieles, carbon mono-oxide, sulfur dioxide and 
nitrogen dioxide at a single site in downtown Toronto from 1973-1974. 
The main findings are following: (1) The electrical conductivity was 
negatively and almost linearly correlated to the individual pollutant 
concentrations, (2) the effects of the individual gaseous pollutants (i.e. 
carbon monoxide, sulfur dioxide and nitrogen dioxide) on the electrical 
conductivity were superimposed with unequal weighting factors, (3) The 
multi-correlation of the electrical conductivity versus the superimposed 
concentration of the gaseous pollutants and the particulate pollutant 
concentration (i.e. aerosol particles was significant, suggesting that the 
atmospheric electrical conductivity can be used as a composite air 
pollution indicator for detecting total air pollution and (4) The weather 
conditions did not seem to have any significant impact on the above 
correlation, although they had influenced both electrical conductivity and 
pollutant concentrations. 
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In the Retezat Mountains concentrations of O3, NO2 and SO2 in 
summer season 2000-2002 were low and below toxicity levels for forest 
trees. While NH3 concentrations were tow in 2000, the 2001 and 2002 
concentrations were elevated indicating possibility for increased N 
deposition to forest stands. More than 90% of the rain events were acidic 
with pH values <5.5, contributing to increased acidity of soils. Crown 
condition of Norway spruce (Picea abies) and European beech (Fagus 
sylvatica) was good, however, defoliation described as > 25% of foliage 
injured increased from 9.1% in 2000 to 16.1% in 2002. Drought that 
occurred in the southern Carpathians between fall 2000 and 2002 and 
frequent acidic rainfalls could cause the observed decline of forest 
condition. Both Norway spruce and European beech with higher 
defoliation had lower annual radial increments compared to the trees with 
low defoliation. Ambient O3 levels found in the Retezat did not affect 
crown condition of Norway spruce or European beech (Akula Venkatram, 
2006). 

Fundamentally, A is related to concentrations of air pollutants 
through destruction processes for small air ions due to recombination (or 
neutralization) of small ions with opposite signs of charges and 
attachment (or adhesion) of the small ions to large particles, as described 
in the following ion budget equation (Mathpal, 1982): 

■ ■ 

^ = q — an^ — ^nN ' * (5.7) 


Where, q is the small ion formation rate by radioactive and cosmic 
sources, a the recombination coefficient, p the attachment coefficient, N 
the number density of large particles and t the time. In destroying small 
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air ions, both molecular (or gaseous) and aerosol (or particulate) air 
pollutants with various chemical characteristics and sizes are believed to 
be active contributors for reducing both and (Paoletti, 1989; Retails, 
1991; Sheftel, 1972). 


The present study is intended for consideration of both individual 
and collective effects of important air pollutants based on experimental 
measurements at a single in Toronto for two consecutive years. Its 
investigation will cover (1) variation patterns of A and air pollution, (2) 
statistical relationships between A and concentrations of important air 
pollutants, and (3) influences of the weather on the relationships. 


6.10 CONOUCTiyimr AND IONIC CONCENTRATION: 


5.10.1 Sources of Ions in the Atmosphere: 


Uranium and Thorium are distributed in the earth’s crust. The 


Radon (Rn) and 'Thoron (Tn) gases, thus formed will diffuse through 
interstices or soil capillaries into the atmosphere where they decay. The 
decay products soon attach themselves to the fine particles (aerosols) in 
the atmosphere. Hence, the fine particles found in the atmosphere possess 
a small but measurable amount of natural radioactivity. Thus the chief 
ionizing agents in the lower atmosphere are the penetrating radiation from 
the radioactive substances in the earth’s crust, radiation from the 
radioactive substances in the atmosphere and the cosmic rays. 


The impact of urban street lay out on local atmospheric 
environment is studied by Xie Xiaomin (2006). 


The y? radiation from the earth’s surface will give about 0.3x10^ 
pairs of ions (PI) per m^ per sec, with radiation about 3.2 x io^ PI m ^ sec \ 


IS 
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a activity of the air gives about 4.4 xlO^ PI m'^ sec'’ with /? and y 
radiations about 0.03x10^ PI m‘^ see’' and 0.15xio^ PI m‘^ sec’, 
respectively. Thus the total ionization rate from the radioactivity is about 
8x10 PI m sec (Bricard, 1965). Table 5.4 gives an average vertical 
distribution of soil radiation (Israel, 1971). 

Contribution ot radioactivity towards ionization over oceans is 
insignificant. The electric cosmic ray component at the ground level will 
give about 2x10^ PI m ^ sec The ionization due to solar cosmic rays is 
essentially a high altitude phenomenon and is important only at the time 
of solar disturbances. On the whole the ionization rate at the ground level 
is about 10x10® Plm'^ sec"'. 

The characterization of atmospheric dry deposition particulates by 
Cobb (1970). 

Table - 5.4 

RADIOACTIVE RADIATION FROM THE SOIL AS A 
FUNCTION OF ALTITUDE 


Altitude above the ground (m) 

___ ^ ^ _ 
_ 

__ . 

__ _ _ 

1000 ' 


Soil radiation in percentage 
97 

’ ^ 83 

_ _____ __ 

_ ^ _ 


0.1 
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5.10.2 Conductivity in Land and Sea Breeze: 

Land and sea breeze is a consequence of simultaneous 
involvement of processes like heating expansion of air, development of 
pressure gradient, forces, accelerations and air motions. There are a 
number of studies on concentration of particulate matter and radioactive 
content in the land and sea winds (Misaki, 1977; Akula, 2006). Misaki et 
al. (19771 could notice drastic reduction in radioactive content in the sea 
air. Morita and Ishikawa (1971) recorded continuous increase in 
conductivities as their relationship, moved away from shore. Blanchard, 
(1966) had obtained high values P.G. and positive space charge in the sea 
wind opposite to the land wind. 

5. 10.3 Atmospheric Electric Conductivity: 

In this section we describe how in our model we construct global 
values of the atmospheric electric conductivity. This includes taking 
cognizance of the distribution of aerosol particles at the earth’s surface, 
its orography, the ionization by the cosmic rays and radioactivity and the 
corona ions beneath the active thunderclouds. The ionization rate by 
cosmic rays and radioactivity has been calculated from the empirical 
model of Makino and Ogawa (1985). 

5.10.4 Electrical Conductivity in and Above the 
Thundercloud: 

There is a considerable uncertainty with regard to an appropriate 
value for conductivity within the Thundercloud. Inside the cloud the ions 
may have extremely low mobility. 
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Size fractionated speciation of nitrate and sulfate aerosol in a 
subtropical industrial environment was found (Sudhir Kumar Pandey, 
2006). 

The regions of shielding charge that are observed at cloud 
boundaries are evidence to suggest these low conductivities. Supporting 
this are the measurements of Gish and Wait (1950) around the 
thunderstorm cells Freier (1962), Evans (1969) and Kamra (1979) on the 
other hand suggested that the conductivity within the thundercloud might 
be as high as 10 to 18 times that of the clear air conductivity. In this 
present study the conductivity has been calculated. 

Atmospheric environmental protection was observed in China, 
current status, developmental trend and research emphasis (K. Madhavi 
Latha, 2006). 

The effective attachment coefficient to the cloud particles 
calculated from (Volmer, 1939). 

= 4jiR,DN, (5.8) 

Where Rc is the mean radius of the cloud particles; Nc , its number 
density and D, the diffusion coefficient of ions. 

Variation of D with height is assumed as followed (Makino, 1985) 

D = 4.2 X lQ~^exp[9.16 x 10“^(Z - ZG)]m^S-^ (5.9) 

The altitude profile of atmospheric conductivity (Ac) can be 
approximated by 
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~ Zq^/S(;] (5.10) 

Where, 

Sc = (5.11) 

Zn is height above the ground of negative charged region of 
thunderstorm; the ambient conductivity at the altitude Zn- 

The conductivity beneath the active thunderstorm at any time t is 
given by 

(5.12) 

Where, the conductivity, A(z,di,(pj) at any grid (ixj) has been 
calculated from 

A(z,0i,(Pi) = 2n(z,ei,^l)eii(z) (5.13) 

Where, z is the altitude; 9i the co-latitude; Oj, the longitude; /z(z), 
the mobility of the ions and e, the electronic change. 

The conductivity is determined by the steady state ion balance 
equation: 

q — an^ — (^NqU — 0 (5.14) 

_L _L 4- J- ' (5.15) 

% % T/? 

Where, the time constant associated with ion production (Xp), ion 
recombination (Tq.) and ion attachment (t^) are given by: 
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J_ = i 

Xp n 


1 

— = an. 





( 5 . 16 ) 


5.1 1 EESULTS AND DISCUSSION: 

A conductivity above the thundercloud has been calculated. Fig. 
5.6 (a) shows the vertical profile of ambient conductivity calculated at the 
equator, which has been compared with the observations of Gish and 
Wait (1950), Woessner et ah, (1958), Morita et ah, (1971), Widell et. ah, 
(1976) and Rosen and Hofmann, (1981). The altitude profile of 
conductivity at the time when thundercloud is active has been shown in 
fig. (5.6 b). Value of AdZ) has been calculated fro m equation (5.101 . 

The LT variation of conductivity beneath the thundercloud due to 
coronal ions has been shown in fig.( 5.6 c). 

The monitoring for Pb and Cd pollution using feral pigeons in 
rural, urban and industrial environments of Korea was studied (Dong-Ha, 
2006). 

Our calculations show that increase in AP concentration to five 



time the average concentration over ocean decrease the conductivity by 
48% as shown in fig. (5.7). 

In calculating the effective attachment coefficient we have used 
equations (5.8) and (5.9). Mason (1972) Chiu and Klett (1976) assumed 
the value of R. and N. as 10 ;im and 2 x lo' m’’, respectively. 
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fig. 5.6 : (a) Vertical profile of ambient conductivity at th e equator. 

(b) Effect of coronol ions on the conductivity. 

(c) Local time variation of conductivity beneath active thundercloud. 
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The calculated value as the function of are shown in Table 
5.5 the graph between (jum) and effective attachment coefficient 
(xIO''^) m^s"' are plotted in fig. (5.7). The variation of and D as the 
function of (Z-Zq) are given in Table 5.6. From these tables we note that 
the value of pc> and D vary exponentially with (Z-Zq) but pc varies 
linearly with D. The altitude profile of atmospheric electric conductivity, 
Ac(Z) has calculated using eq. (5.10) and eq. (5.11) at the time when 
thundercloud is active. These calculations have shown in Table 5.7. We 
have taken conductivity near earth’s surface A^CZ) to be one order of 
magnitude higher than its ambient value(/lciv = 4.4 x We 

find that A^CZ) decreases exponentially with increase in Zn. At about Zn 
= 4.5 km, Ac(Z) is minimum. At higher altitudes than Zn =4.5; Ac(Z) 
increases. Thus the values of A^CZ) depend on ambient conductivity Ag 
and(Z — Zq). Also at altitudes higher than 6 km both A^ and Ac(Z) 
increase and remain the same at a given height. The calculated values of 
Ac(z) are in good agreement with the experimental work of Gish and 
Wait (1950), Woessner et al., (1958), Morita et al. (1971), Rosen and 
Hofmann (1981), Widell et al. (1976) and Markson et al. (1981). 

CFD simulations of the wind environment around an airport 
terminal building was observed (P. Neofytou, 2006). 
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TABLE-5.5 

CALCliLA l ED VALUES OF EFFECTIVE ATTACHMENT 
COEFFICIENT p, AS THE FUNCTION OF MEAN RADIUS OF 
CLOUD PARTICLES (R^) AT (Z-Zg) = 1 KM. 


S.No. 

Rcifim) 

Pc(^ 10 ^ 


2 

1.157 

2 . 

6 

i 

3.471 

i 3. 

1 

10 

5.784 

4. 

14 

8.098 

5. 

18 

10.412 

6. 

22 

\2J25 

7 

26 

15.039 

8. 

30 

17.353 

9. 

34 

19.666 

wmgm 

38 

21.980 

IL 

42 

24.294 

12. 

46 

26.607 

13. 

50 

28.921 
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TABLE-S.6 

CALCULATED VALUES OF D, THE DIFFUSION COEFFICIENT OF 
IONS AND jSe THE EFFECTIVE ATTACHMENT COEFFICIENT AS 
THE FUNCTIONS OF (Z-Zc) AND R, = 20fim CONSTANT 


S.No. 

(Z-Z(;)kni 

D(xlO'^)mV' 


1. 

2 

4.2008 

1.0558 

2. 

6 

4.2023 

1.0562 

3. 

1 

10 

4.2038 

1.0565 

4. 

14 

4.2054 

1.0569 

5. 

18 

4.2069 

1.0573 

6. 

22 

4.2085 

1.0577 

7. 

26 

4.2100 

1.0581 

8. 

30 

4.2116 

1.0585 

1 

34 

4.2123 

1.0586 


38 

4.2146 

1.0593 

11 . 

42 

4.2162 

1.0596 

12. 

46 

4.2177 

1.0600 

13. 

50 

4.2193 

1.0604 


( 272 > 

^'-xrrrrrrrtl^ 



















5.12 CONCLUSION: 

From above study, we infer that there is a decrease in the electrical 
conductivity with altitude near the ground. It is minimum at about 4.5 km 
above the ground and it again increase above the height. The decrease in 
electrical conductivity may be attributed to the decrease in the 
concentration of ion because the ions get attached with aerosol particles 
or CCN and other particles. But at higher altitude, the concentration of 
ions increases due to ionization caused by galactic cosmic rays. This 

increase of ions increases the electric conductivity. 
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